Search

Phyllodes tumours (PTs) are rare fibroepithelial lesions of the breast that are classified as benign, borderline, or malignant. As little is known about the molecular underpinnings of PTs, current diagnosis relies on histological examination. However, accurate classification is often difficult, particularly for distinguishing borderline from malignant PTs. Furthermore, PTs can be misdiagnosed as other tumour types with shared histological features, such as fibroadenoma and metaplastic breast cancers. As DNA methylation is a recognised hallmark of many cancers, we hypothesised that DNA methylation could provide novel biomarkers for diagnosis and tumour stratification in PTs, whilst also allowing insight into the molecular aetiology of this otherwise understudied tumour. We generated whole-genome methylation data using the Illumina EPIC microarray in a novel PT cohort (n = 33) and curated methylation microarray data from published datasets including PTs and other potentially histopathologically similar tumours (total n = 817 samples). Analyses revealed that PTs have a unique methylome compared to normal breast tissue and to potentially histopathologically similar tumours (metaplastic breast cancer, fibroadenoma and sarcomas), with PT-specific methylation changes enriched in gene sets involved in KRAS signalling and epithelial-mesenchymal transition. Next, we identified 53 differentially methylated regions (DMRs) (false discovery rate < 0.05) that specifically delineated malignant from non-malignant PTs. The top DMR in both discovery and validation cohorts was hypermethylation at the HSD17B8 CpG island promoter. Matched PT single-cell expression data showed that HSD17B8 had minimal expression in fibroblast (putative tumour) cells. Finally, we created a methylation classifier to distinguish PTs from metaplastic breast cancer samples, where we revealed a likely misdiagnosis for two TCGA metaplastic breast cancer samples. In conclusion, DNA methylation alterations are associated with PT histopathology and hold the potential to improve our understanding of PT molecular aetiology, diagnostics, and risk stratification. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland., (© 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.)

Background: The Illumina family of Infinium Methylation BeadChip microarrays has been widely used over the last 15 years for genome-wide DNA methylation profiling, including large-scale and population-based studies, due to their ease of use and cost effectiveness. Succeeding the popular HumanMethylationEPIC BeadChip (EPICv1), the recently released Infinium MethylationEPIC v2.0 BeadChip (EPICv2) claims to extend genomic coverage to more than 935,000 CpG sites. Here, we comprehensively characterise the reproducibility, reliability and annotation of the EPICv2 array, based on bioinformatic analysis of both manifest data and new EPICv2 data from diverse biological samples., Results: We find a high degree of reproducibility with EPICv1, evidenced by comparable sensitivity and precision from empirical cross-platform comparison incorporating whole genome bisulphite sequencing (WGBS), and high correlation between technical sample replicates, including between samples with DNA input levels below the manufacturer's recommendation. We provide a full assessment of probe content, evaluating genomic distribution and changes from previous array versions. We characterise EPICv2's new feature of replicated probes and provide recommendations as to the superior probes. In silico analysis of probe sequences demonstrates that probe cross-hybridisation remains a significant problem in EPICv2. By mapping the off-target sites at single nucleotide resolution and comparing with WGBS we show empirical evidence for preferential off-target binding., Conclusions: Overall, we find EPICv2 a worthy successor to the previous Infinium methylation microarrays, however some technical issues remain. To support optimal EPICv2 data analysis we provide an expanded version of the EPICv2 manifest to aid researchers in understanding probe design, data processing, choosing appropriate probes for analysis and for integration with methylation datasets from previous versions of the Infinium Methylation BeadChip., (© 2024. The Author(s).)

Three-dimensional (3D) epigenome remodeling is an important mechanism of gene deregulation in cancer. However, its potential as a target to counteract therapy resistance remains largely unaddressed. Here, we show that epigenetic therapy with decitabine (5-Aza-mC) suppresses tumor growth in xenograft models of pre-clinical metastatic estrogen receptor positive (ER+) breast tumor. Decitabine-induced genome-wide DNA hypomethylation results in large-scale 3D epigenome deregulation, including de-compaction of higher-order chromatin structure and loss of boundary insulation of topologically associated domains. Significant DNA hypomethylation associates with ectopic activation of ER-enhancers, gain in ER binding, creation of new 3D enhancer-promoter interactions and concordant up-regulation of ER-mediated transcription pathways. Importantly, long-term withdrawal of epigenetic therapy partially restores methylation at ER-enhancer elements, resulting in a loss of ectopic 3D enhancer-promoter interactions and associated gene repression. Our study illustrates the potential of epigenetic therapy to target ER+ endocrine-resistant breast cancer by DNA methylation-dependent rewiring of 3D chromatin interactions, which are associated with the suppression of tumor growth., (© 2024. The Author(s).)

Gene expression noise is known to promote stochastic drug resistance through the elevated expression of individual genes in rare cancer cells. However, we now demonstrate that chemoresistant neuroblastoma cells emerge at a much higher frequency when the influence of noise is integrated across multiple components of an apoptotic signaling network. Using a JNK activity biosensor with longitudinal high-content and in vivo intravital imaging, we identify a population of stochastic, JNK-impaired, chemoresistant cells that exist because of noise within this signaling network. Furthermore, we reveal that the memory of this initially random state is retained following chemotherapy treatment across a series of in vitro, in vivo, and patient models. Using matched PDX models established at diagnosis and relapse from individual patients, we show that HDAC inhibitor priming cannot erase the memory of this resistant state within relapsed neuroblastomas but improves response in the first-line setting by restoring drug-induced JNK activity within the chemoresistant population of treatment-naïve tumors.

Background: Prostate cancer is a clinically heterogeneous disease with a subset of patients rapidly progressing to lethal-metastatic prostate cancer. Current clinicopathological measures are imperfect predictors of disease progression. Epigenetic changes are amongst the earliest molecular changes in tumourigenesis. To find new prognostic biomarkers to enable earlier intervention and improved outcomes, we performed methylome sequencing of DNA from patients with localised prostate cancer and long-term clinical follow-up., Methods: We used whole-genome bisulphite sequencing (WGBS) to comprehensively map and compare DNA methylation of radical prostatectomy tissue between patients with lethal disease (n = 7) and non-lethal (n = 8) disease (median follow-up 19.5 years). Validation of differentially methylated regions (DMRs) was performed in an independent cohort (n = 185, median follow-up 15 years) using targeted multiplex bisulphite sequencing of candidate regions. Survival was assessed via univariable and multivariable analyses including clinicopathological measures (log-rank and Cox regression models)., Results: WGBS data analysis identified cancer-specific methylation patterns including CpG island hypermethylation, and hypomethylation of repetitive elements, with increasing disease risk. We identified 1420 DMRs associated with prostate cancer-specific mortality (PCSM), which showed enrichment for gene sets downregulated in prostate cancer and de novo methylated in cancer. Through comparison with public prostate cancer datasets, we refined the DMRs to develop an 18-gene prognostic panel. Applying this panel to an independent cohort, we found significant associations between PCSM and hypermethylation at EPHB3, PARP6, TBX1, MARCH6 and a regulatory element within CACNA2D4. Strikingly in a multivariable model, inclusion of CACNA2D4 methylation was a better predictor of PCSM versus grade alone (Harrell's C-index: 0.779 vs. 0.684)., Conclusions: Our study provides detailed methylome maps of non-lethal and lethal prostate cancer and identifies novel genic regions that distinguish these patient groups. Inclusion of our DNA methylation biomarkers with existing clinicopathological measures improves prognostic models of prostate cancer mortality, and holds promise for clinical application., (© 2022 The Authors. Clinical and Translational Medicine published by John Wiley & Sons Australia, Ltd on behalf of Shanghai Institute of Clinical Bioinformatics.)

Epigenetic therapies remain a promising, but still not widely used, approach in the management of patients with cancer. To date, the efficacy and use of epigenetic therapies has been demonstrated primarily in the management of haematological malignancies, with limited supportive data in solid malignancies. The most studied epigenetic therapies in breast cancer are those that target DNA methylation and histone modification; however, none have been approved for routine clinical use. The majority of pre-clinical and clinical studies have focused on triple negative breast cancer (TNBC) and hormone-receptor positive breast cancer. Even though the use of epigenetic therapies alone in the treatment of breast cancer has not shown significant clinical benefit, these therapies show most promise in use in combinations with other treatments. With improving technologies available to study the epigenetic landscape in cancer, novel epigenetic alterations are increasingly being identified as potential biomarkers of response to conventional and epigenetic therapies. In this review, we describe epigenetic targets and potential epigenetic biomarkers in breast cancer, with a focus on clinical trials of epigenetic therapies. We describe alterations to the epigenetic landscape in breast cancer and in treatment resistance, highlighting mechanisms and potential targets for epigenetic therapies. We provide an updated review on epigenetic therapies in the pre-clinical and clinical setting in breast cancer, with a focus on potential real-world applications. Finally, we report on the potential value of epigenetic biomarkers in diagnosis, prognosis and prediction of response to therapy, to guide and inform the clinical management of breast cancer patients.

Neoadjuvant chemotherapy (NAC) is used to treat triple-negative breast cancer (TNBC) prior to resection. Biomarkers that accurately predict a patient's response to NAC are needed to individualise therapy and avoid chemotoxicity from unnecessary chemotherapy. We performed whole-genome DNA methylation profiling on diagnostic TNBC biopsy samples from the Sequential Evaluation of Tumours Undergoing Preoperative (SETUP) NAC study. We found 9 significantly differentially methylated regions (DMRs) at diagnosis which were associated with response to NAC. We show that 4 of these DMRs are associated with TNBC overall survival (P < 0.05). Our results highlight the potential of DNA methylation biomarkers for predicting NAC response in TNBC., (© 2021. The Author(s).)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Summary: DNA methylation patterns in a cell are associated with gene expression and the phenotype of a cell, including disease states. Bisulphite PCR sequencing is commonly used to assess the methylation profile of genomic regions between different cells. Here we have developed MethPanel, a computational pipeline with an interactive graphical interface to rapidly analyse multiplex bisulphite PCR sequencing data. MethPanel comprises a complete analysis workflow from genomic alignment to DNA methylation calling and supports an unlimited number of PCR amplicons and input samples. MethPanel offers important and unique features, such as calculation of an epipolymorphism score and bisulphite PCR bias correction capabilities, and is designed so that the methylation data from all samples can be processed in parallel. The outputs are automatically forwarded to a shinyApp for convenient display, visualization and remotely sharing data with collaborators and clinicians., Availabilityand Implementation: MethPanel is freely available at https://github.com/thinhong/MethPanel., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press.)

Modification of the cancer-associated chromatin landscape in response to therapeutic DNA damage influences gene expression and contributes to cell fate. The central histone mark H2Bub1 results from addition of a single ubiquitin on lysine 120 of histone H2B and is an important regulator of gene expression. Following treatment with a platinum-based chemotherapeutic, there is a reduction in global levels of H2Bub1 accompanied by an increase in levels of the tumor suppressor p53. Although total H2Bub1 decreases following DNA damage, H2Bub1 is enriched downstream of transcription start sites of specific genes. Gene-specific H2Bub1 enrichment was observed at a defined group of genes that clustered into cancer-related pathways and correlated with increased gene expression. H2Bub1-enriched genes encompassed fifteen p53 target genes including PPM1D, BTG2, PLK2, MDM2, CDKN1A and BBC3, genes related to ERK/MAPK signalling, those participating in nucleotide excision repair including XPC, and genes involved in the immune response and platinum drug resistance including POLH. Enrichment of H2Bub1 at key cancer-related genes may function to regulate gene expression and influence the cellular response to therapeutic DNA damage.

There is a major clinical need for accurate biomarkers for prostate cancer prognosis, to better inform treatment strategies and disease monitoring. Current clinically recognised prognostic factors, including prostate-specific antigen (PSA) levels, lack sensitivity and specificity in distinguishing aggressive from indolent disease, particularly in patients with localised intermediate grade prostate cancer. There has therefore been a major focus on identifying molecular biomarkers that can add prognostic value to existing markers, including investigation of DNA methylation, which has a known role in tumorigenesis. In this review, we will provide a comprehensive overview of the current state of DNA methylation biomarker studies in prostate cancer prognosis, and highlight the advances that have been made in this field. We cover the numerous studies into well-established candidate genes, and explore the technological transition that has enabled hypothesis-free genome-wide studies and the subsequent discovery of novel prognostic genes.

Background: DNA methylation is a well-studied epigenetic mark that is frequently altered in diseases such as cancer, where specific changes are known to reflect the type and severity of the disease. Therefore, there is a growing interest in assessing the clinical utility of DNA methylation as a biomarker for diagnosing disease and guiding treatment. The development of an accurate loci-specific methylation assay, suitable for use on low-input clinical material, is crucial for advancing DNA methylation biomarkers into a clinical setting. A targeted multiplex bisulphite PCR sequencing approach meets these needs by allowing multiple DNA methylated regions to be interrogated simultaneously in one experiment on limited clinical material., Results: Here, we provide an updated protocol and recommendations for multiplex bisulphite PCR sequencing (MBPS) assays for target DNA methylation analysis. We describe additional steps to improve performance and reliability: (1) pre-sequencing PCR optimisation which includes assessing the optimal PCR cycling temperature and primer concentration and (2) post-sequencing PCR optimisation to achieve uniform coverage of each amplicon. We use a gradient of methylated controls to demonstrate how PCR bias can be assessed and corrected. Methylated controls also allow assessment of the sensitivity of methylation detection for each amplicon. Here, we show that the MBPS assay can amplify as little as 0.625 ng starting DNA and can detect methylation differences of 1% with a sequencing coverage of 1000 reads. Furthermore, the multiplex bisulphite PCR assay can comprehensively interrogate multiple regions on 1-5 ng of formalin-fixed paraffin-embedded DNA or circulating cell-free DNA., Conclusions: The MBPS assay is a valuable approach for assessing methylated DNA regions in clinical samples with limited material. The optimisation and additional quality control steps described here improve the performance and reliability of this method, advancing it towards potential clinical applications in biomarker studies.

Background: Prostate cancer changes the phenotype of cells within the stromal microenvironment, including fibroblasts, which in turn promote tumour progression. Functional changes in prostate cancer-associated fibroblasts (CAFs) coincide with alterations in DNA methylation levels at loci-specific regulatory regions. Yet, it is not clear how these methylation changes compare across CAFs from different patients. Therefore, we examined the consistency and prognostic significance of genome-wide DNA methylation profiles between CAFs from patients with different grades of primary prostate cancer., Results: We used Infinium MethylationEPIC BeadChips to evaluate genome-wide DNA methylation profiles from 18 matched CAFs and non-malignant prostate tissue fibroblasts (NPFs) from men with moderate to high grade prostate cancer, as well as five unmatched benign prostate tissue fibroblasts (BPFs) from men with benign prostatic hyperplasia. We identified two sets of differentially methylated regions (DMRs) in patient CAFs. One set of DMRs reproducibly differed between CAFs and fibroblasts from non-malignant tissue (NPFs and BPFs). Indeed, more than 1200 DMRs consistently changed in CAFs from every patient, regardless of tumour grade. The second set of DMRs varied between CAFs according to the severity of the tumour. Notably, hypomethylation of the EDARADD promoter occurred specifically in CAFs from high-grade tumours and correlated with increased transcript abundance and increased EDARADD staining in patient tissue. Across multiple cohorts, tumours with low EDARADD DNA methylation and high EDARADD mRNA expression were consistently associated with adverse clinical features and shorter recurrence free survival., Conclusions: We identified a large set of DMRs that are commonly shared across CAFs regardless of tumour grade and outcome, demonstrating highly consistent epigenome changes in the prostate tumour microenvironment. Additionally, we found that CAFs from aggressive prostate cancers have discrete methylation differences compared to CAFs from moderate risk prostate cancer. Together, our data demonstrates that the methylome of the tumour microenvironment reflects both the presence and the severity of the prostate cancer and, therefore, may provide diagnostic and prognostic potential.

Endocrine therapy resistance frequently develops in estrogen receptor positive (ER+) breast cancer, but the underlying molecular mechanisms are largely unknown. Here, we show that 3-dimensional (3D) chromatin interactions both within and between topologically associating domains (TADs) frequently change in ER+ endocrine-resistant breast cancer cells and that the differential interactions are enriched for resistance-associated genetic variants at CTCF-bound anchors. Ectopic chromatin interactions are preferentially enriched at active enhancers and promoters and ER binding sites, and are associated with altered expression of ER-regulated genes, consistent with dynamic remodelling of ER pathways accompanying the development of endocrine resistance. We observe that loss of 3D chromatin interactions often occurs coincidently with hypermethylation and loss of ER binding. Alterations in active A and inactive B chromosomal compartments are also associated with decreased ER binding and atypical interactions and gene expression. Together, our results suggest that 3D epigenome remodelling is a key mechanism underlying endocrine resistance in ER+ breast cancer.

The architectural protein CTCF is a mediator of chromatin conformation, but how CTCF binding to DNA is orchestrated to maintain long-range gene expression is poorly understood. Here we perform RNAi knockdown to reduce CTCF levels and reveal a shared subset of CTCF-bound sites are robustly resistant to protein depletion. The 'persistent' CTCF sites are enriched at domain boundaries and chromatin loops constitutive to all cell types. CRISPR-Cas9 deletion of 2 persistent CTCF sites at the boundary between a long-range epigenetically active (LREA) and silenced (LRES) region, within the Kallikrein (KLK) locus, results in concordant activation of all 8 KLK genes within the LRES region. CTCF genome-wide depletion results in alteration in Topologically Associating Domain (TAD) structure, including the merging of TADs, whereas TAD boundaries are not altered where persistent sites are maintained. We propose that the subset of essential CTCF sites are involved in cell-type constitutive, higher order chromatin architecture.

As one of the most abundant and well-studied epigenetic modifications, DNA methylation plays an essential role in normal development and cellular biology. Global alterations to the DNA methylation landscape contribute to alterations in the transcriptome and deregulation of cellular pathways. Indeed, improved methods to study DNA methylation patterning and dynamics at base pair resolution and across individual DNA molecules on a genome-wide scale has highlighted the scope of change to the DNA methylation landscape in disease states, particularly during tumorigenesis. More recently has been the development of DNA hydroxymethylation profiling techniques, which allows differentiation between 5mC and 5hmC profiles and provides further insights into DNA methylation dynamics and remodeling in tumorigenesis. In this review, we describe the distribution of DNA methylation and DNA hydroxymethylation in different genomic contexts, first in normal cells, and how this is altered in cancer. Finally, we discuss DNA methylation profiling technologies and the most recent advances in single-cell methods, bisulfite-free approaches and ultra-long read sequencing techniques., (© 2019 The Author(s).)

Adipocytes support key metabolic and endocrine functions of adipose tissue. Lipid is stored in two major classes of depots, namely visceral adipose (VA) and subcutaneous adipose (SA) depots. Increased visceral adiposity is associated with adverse health outcomes, whereas the impact of SA tissue is relatively metabolically benign. The precise molecular features associated with the functional differences between the adipose depots are still not well understood. Here, we characterised transcriptomes and methylomes of isolated adipocytes from matched SA and VA tissues of individuals with normal BMI to identify epigenetic differences and their contribution to cell type and depot-specific function. We found that DNA methylomes were notably distinct between different adipocyte depots and were associated with differential gene expression within pathways fundamental to adipocyte function. Most striking differential methylation was found at transcription factor and developmental genes. Our findings highlight the importance of developmental origins in the function of different fat depots.

Motivation: A synoptic view of the human genome benefits chiefly from the application of nucleic acid sequencing and microarray technologies. These platforms allow interrogation of patterns such as gene expression and DNA methylation at the vast majority of canonical loci, allowing granular insights and opportunities for validation of original findings. However, problems arise when validating against a "gold standard" measurement, since this immediately biases all subsequent measurements towards that particular technology or protocol. Since all genomic measurements are estimates, in the absence of a "gold standard" we instead empirically assess the measurement precision and sensitivity of a large suite of genomic technologies via a consensus modelling method called the row-linear model. This method is an application of the American Society for Testing and Materials Standard E691 for assessing interlaboratory precision and sources of variability across multiple testing sites. Both cross-platform and cross-locus comparisons can be made across all common loci, allowing identification of technology- and locus-specific tendencies., Results: We assess technologies including the Infinium MethylationEPIC BeadChip, whole genome bisulfite sequencing (WGBS), two different RNA-Seq protocols (PolyA+ and Ribo-Zero) and five different gene expression array platforms. Each technology thus is characterised herein, relative to the consensus. We showcase a number of applications of the row-linear model, including correlation with known interfering traits. We demonstrate a clear effect of cross-hybridisation on the sensitivity of Infinium methylation arrays. Additionally, we perform a true interlaboratory test on a set of samples interrogated on the same platform across twenty-one separate testing laboratories., Availability and Implementation: A full implementation of the row-linear model, plus extra functions for visualisation, are found in the R package consensus at https://github.com/timpeters82/consensus., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2018. Published by Oxford University Press.)

Background: ATP-dependent chromatin remodelling complexes are responsible for establishing and maintaining the positions of nucleosomes. Chromatin remodellers are targeted to chromatin by transcription factors and non-coding RNA to remodel the chromatin into functional states. However, the influence of chromatin remodelling on shaping the functional epigenome is not well understood. Moreover, chromatin remodellers have not been extensively explored as a collective group across two-dimensional and three-dimensional epigenomic layers., Results: Here, we have integrated the genome-wide binding profiles of eight chromatin remodellers together with DNA methylation, nucleosome positioning, histone modification and Hi-C chromosomal contacts to reveal that chromatin remodellers can be stratified into two functional groups. Group 1 (BRG1, SNF2H, CHD3 and CHD4) has a clear preference for binding at 'actively marked' chromatin and Group 2 (BRM, INO80, SNF2L and CHD1) for 'repressively marked' chromatin. We find that histone modifications and chromatin architectural features, but not DNA methylation, stratify the remodellers into these functional groups., Conclusions: Our findings suggest that chromatin remodelling events are synchronous and that chromatin remodellers themselves should be considered simultaneously and not as individual entities in isolation or necessarily by structural similarity, as they are traditionally classified. Their coordinated function should be considered by preference for chromatin features in order to gain a more accurate and comprehensive picture of chromatin regulation.

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., (Copyright © 2019 Elsevier Inc. All rights reserved.)

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.

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., (© 2018 Pidsley et al.; Published by Cold Spring Harbor Laboratory Press.)

Epigenetic regulation plays a critical role in gene expression, cellular differentiation, and disease. There is a complex interplay between the different layers of epigenetic information, including DNA methylation, nucleosome positions, histone modifications, histone variants, and other important epigenetic regulators. The different modifications do not act independently of each other and their relationship plays an important role in governing the regulation of the epigenome. Of these, DNA methylation is the best-studied epigenetic modification in mammals. However, the direct relationship between DNA methylation and chromatin modifications has been difficult to unravel with existing technologies, with epigenome-wide integration studies still based on "overlaying" independent chromatin modification and DNA methylation maps. Bisulphite sequencing enables the methylation state of every cytosine residue to be analyzed across a given molecule in a strand-specific context. Here, we describe a direct approach to interrogating the DNA methylation status of specific chromatin-marked DNA, using high-throughput sequencing of bisulphite-treated chromatin immunoprecipitated DNA (BisChIP-seq). This combined approach enables the exquisite relationship between chromatin-modified DNA or transcription factor-associated DNA and the methylation state of each targeted allele to be directly interrogated. BisChIP-Seq can now be widely applied genome-wide to further understand the molecular relationship between DNA methylation and other important epigenetic regulators.

Acetylation of the histone variant H2A.Z (H2A.Zac) occurs at active promoters and is associated with oncogene activation in prostate cancer, but its role in enhancer function is still poorly understood. Here we show that H2A.Zac containing nucleosomes are commonly redistributed to neo-enhancers in cancer resulting in a concomitant gain of chromatin accessibility and ectopic gene expression. Notably incorporation of acetylated H2A.Z nucleosomes is a pre-requisite for activation of Androgen receptor (AR) associated enhancers. H2A.Zac nucleosome occupancy is rapidly remodeled to flank the AR sites to initiate the formation of nucleosome-free regions and the production of AR-enhancer RNAs upon androgen treatment. Remarkably higher levels of global H2A.Zac correlate with poorer prognosis. Altogether these data demonstrate the novel contribution of H2A.Zac in activation of newly formed enhancers in prostate cancer.

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.

Cancer is characterised by DNA hypermethylation and gene silencing of CpG island-associated promoters, including tumour-suppressor genes. The methyl-CpG-binding domain (MBD) family of proteins bind to methylated DNA and can aid in the mediation of gene silencing through interaction with histone deacetylases and histone methyltransferases. However, the mechanisms responsible for eliciting CpG island hypermethylation in cancer, and the potential role that MBD proteins play in modulation of the methylome remain unclear. Our previous work demonstrated that MBD2 preferentially binds to the hypermethylated GSTP1 promoter CpG island in prostate cancer cells. Here, we use functional genetic approaches to investigate if MBD2 plays an active role in reshaping the DNA methylation landscape at this locus and genome-wide. First, we show that loss of MBD2 results in inhibition of both maintenance and spread of de novo methylation of a transfected construct containing the GSTP1 promoter CpG island in prostate cancer cells and Mbd2-/- mouse fibroblasts. De novo methylation was rescued by transient expression of Mbd2 in Mbd2-/- cells. Second, we show that MBD2 depletion triggers significant hypomethylation genome-wide in prostate cancer cells with concomitant loss of MBD2 binding at promoter and enhancer regulatory regions. Finally, CpG islands and shores that become hypomethylated after MBD2 depletion in LNCaP cancer cells show significant hypermethylation in clinical prostate cancer samples, highlighting a potential active role of MBD2 in promoting cancer-specific hypermethylation. Importantly, co-immunoprecipiation of MBD2 shows that MBD2 associates with DNA methyltransferase enzymes 1 and 3A. Together our results demonstrate that MBD2 has a critical role in 'rewriting' the cancer methylome at specific regulatory regions.

Background: In recent years the Illumina HumanMethylation450 (HM450) BeadChip has provided a user-friendly platform to profile DNA methylation in human samples. However, HM450 lacked coverage of distal regulatory elements. Illumina have now released the MethylationEPIC (EPIC) BeadChip, with new content specifically designed to target these regions. We have used HM450 and whole-genome bisulphite sequencing (WGBS) to perform a critical evaluation of the new EPIC array platform., Results: EPIC covers over 850,000 CpG sites, including >90 % of the CpGs from the HM450 and an additional 413,743 CpGs. Even though the additional probes improve the coverage of regulatory elements, including 58 % of FANTOM5 enhancers, only 7 % distal and 27 % proximal ENCODE regulatory elements are represented. Detailed comparisons of regulatory elements from EPIC and WGBS show that a single EPIC probe is not always informative for those distal regulatory elements showing variable methylation across the region. However, overall data from the EPIC array at single loci are highly reproducible across technical and biological replicates and demonstrate high correlation with HM450 and WGBS data. We show that the HM450 and EPIC arrays distinguish differentially methylated probes, but the absolute agreement depends on the threshold set for each platform. Finally, we provide an annotated list of probes whose signal could be affected by cross-hybridisation or underlying genetic variation., Conclusion: The EPIC array is a significant improvement over the HM450 array, with increased genome coverage of regulatory regions and high reproducibility and reliability, providing a valuable tool for high-throughput human methylome analyses from diverse clinical samples.

Triple negative breast cancer (TNBC) is the most aggressive breast cancer subtype with no targeted treatment available. Our previous study identified 38 TNBC-specific genes with altered expression comparing tumour to normal samples. This study aimed to establish whether DNA methylation contributed to these expression changes in the same cohort as well as disease progression from primary breast tumour to lymph node metastasis associated with changes in the epigenome. We obtained DNA from 23 primary TNBC samples, 12 matched lymph node metastases, and 11 matched normal adjacent tissues and assayed for differential methylation profiles using Illumina HumanMethylation450 BeadChips. The results were validated in an independent cohort of 70 primary TNBC samples. The expression of 16/38 TNBC-specific genes was associated with alteration in DNA methylation. Novel methylation changes between primary tumours and lymph node metastases, as well as those associated with survival were identified. Altered methylation of 18 genes associated with lymph node metastasis were identified and validated. This study reveals the important role DNA methylation plays in altered gene expression of TNBC-specific genes and lymph node metastases. The novel insights into progression of TNBC to secondary disease may provide potential prognostic indicators for this hard-to-treat breast cancer subtype.

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., (© 2016 Taberlay et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

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.

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.

Background: The clinical utility of DNA methylation as a predictive or prognostic biomarker requires scalable resequencing protocols for bisulfite-converted DNA. Key features of any validation method should be adaptability for low- or high-throughput needs and high reproducibility, and should only require minimal amounts of precious clinical sample as input material. Critically, this method should also deliver robust results when working with bisulfite-converted DNA extracted from formalin-fixed, paraffin-embedded (FFPE) blocks., Results: We report here for the first time on comparison studies between the Fluidigm Access Array system and multiplex assays for multiplex bisulfite PCR resequencing. The requirement of the Fluidigm Access Array system for high template amounts and its sensitivity to variations in template quality rendered it unsuitable for bisulfite PCR applications utilizing FFPE DNA. In response to this limitation, we established a multiplex bisulfite PCR assay capable of delivering robust methylation data using minimal amounts of FFPE clinical DNA. To evaluate the parameters and reproducibility of this assay, 57 amplicons were used to prepare sequencing libraries in triplicate for 13 FFPE tumour samples, as well as a series of 5 methylated controls (0%, 25%, 50%, 75%, and 100%). Analysis of this data demonstrated that this multiplex assay had high reproducibility (mean standard deviation of 1.4% for methylation values), was low cost, required low sample input (50 ng of DNA or less), and could be scaled for both low- and high-throughput needs. Notably, ExoSAP-IT (exonuclease I) treatment to remove residual primers in bisulfite resequencing libraries appeared to degrade the library and generate a high-molecular weight smear which may impact on the degree of methylation assessed., Conclusions: Multiplex bisulfite PCR assays represent a convenient and scalable method for validation and screening of methylated DNA regions from archival FFPE DNA. Moreover, the library construction process detailed here can be rapidly optimized and implemented with a minimal amount of work, can be performed using the standard equipment found in any molecular biology laboratory, and can be easily adapted for use on both genomic DNA and bisulfite DNA applications. However, in preparing bisulfite libraries for sequencing, the use of ExoSAP-IT is not recommended due to potential off-target nuclease effects which may impact downstream methylation analysis.

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.

How DNA methylation is interpreted and influences genome regulation remains largely unknown. Proteins of the methyl-CpG-binding domain (MBD) family are primary candidates for the readout of DNA methylation as they recruit chromatin remodelers, histone deacetylases and methylases to methylated DNA associated with gene repression. MBD protein binding requires both functional MBD domains and methyl-CpGs; however, some MBD proteins also bind unmethylated DNA and active regulatory regions via alternative regulatory domains or interaction with the nucleosome remodeling deacetylase (NuRD/Mi-2) complex members. Mutations within MBD domains occur in many diseases, including neurological disorders and cancers, leading to loss of MBD binding specificity to methylated sites and gene deregulation. Here, we summarize the current state of knowledge about MBD proteins and their role as readers of the epigenome.

Affinity capture of DNA methylation combined with high-throughput sequencing strikes a good balance between the high cost of whole genome bisulfite sequencing and the low coverage of methylation arrays. We present BayMeth, an empirical Bayes approach that uses a fully methylated control sample to transform observed read counts into regional methylation levels. In our model, inefficient capture can readily be distinguished from low methylation levels. BayMeth improves on existing methods, allows explicit modeling of copy number variation, and offers computationally efficient analytical mean and variance estimators. BayMeth is available in the Repitools Bioconductor package.

There are over 28 million CpG sites in the human genome. Assessing the methylation status of each of these sites will be required to understand fully the role of DNA methylation in health and disease. Genome-wide analysis, using arrays and high-throughput sequencing, has enabled assessment of large fractions of the methylome, but each protocol comes with unique advantages and disadvantages. Notably, except for whole-genome bisulfite sequencing, most commonly used genome-wide methods detect <5% of all CpG sites. Here, we discuss approaches for methylome studies and compare genome coverage of promoters, genes, and intergenic regions, and capacity to quantitate individual CpG methylation states. Finally, we examine the extent of published cancer methylomes that have been generated using genome-wide approaches., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

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., (Crown Copyright © 2013. Published by Elsevier Ltd. All rights reserved.)

Deregulation of microRNA (miRNA) expression can have a critical role in carcinogenesis. Here we show in prostate cancer that miRNA-205 (miR-205) transcription is commonly repressed and the MIR-205 locus is hypermethylated. LOC642587, the MIR-205 host gene of unknown function, is also concordantly inactivated. We show that miR-205 targets mediator 1 (MED1, also called TRAP220 and PPARBP) for transcriptional silencing in normal prostate cells, leading to reduction in MED1 mRNA levels, and in total and active phospho-MED1 protein. Overexpression of miR-205 in prostate cancer cells negatively affects cell viability, consistent with a tumor suppressor function. We found that hypermethylation of the MIR-205 locus was strongly related with a decrease in miR-205 expression and an increase in MED1 expression in primary tumor samples (n=14), when compared with matched normal prostate (n=7). An expanded patient cohort (tumor n=149, matched normal n=30) also showed significant MIR-205 DNA methylation in tumors compared with normal, and MIR-205 hypermethylation is significantly associated with biochemical recurrence (hazard ratio=2.005, 95% confidence interval (1.109, 3.625), P=0.02), in patients with low preoperative prostate specific antigen. In summary, these results suggest that miR-205 is an epigenetically regulated tumor suppressor that targets MED1 and may provide a potential biomarker in prostate cancer management.

Epigenetic gene deregulation in cancer commonly occurs through chromatin repression and promoter hypermethylation of tumor-associated genes. However, the mechanism underpinning epigenetic-based gene activation in carcinogenesis is still poorly understood. Here, we identify a mechanism of domain gene deregulation through coordinated long-range epigenetic activation (LREA) of regions that typically span 1 Mb and harbor key oncogenes, microRNAs, and cancer biomarker genes. Gene promoters within LREA domains are characterized by a gain of active chromatin marks and a loss of repressive marks. Notably, although promoter hypomethylation is uncommon, we show that extensive DNA hypermethylation of CpG islands or "CpG-island borders" is strongly related to cancer-specific gene activation or differential promoter usage. These findings have wide ramifications for cancer diagnosis, progression, and epigenetic-based gene therapies., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Developments in microarray and high-throughput sequencing (HTS) technologies have resulted in a rapid expansion of research into epigenomic changes that occur in normal development and in the progression of disease, such as cancer. Not surprisingly, copy number variation (CNV) has a direct effect on HTS read densities and can therefore bias differential detection results. We have developed a flexible approach called ABCD-DNA (affinity-based copy-number-aware differential quantitative DNA sequencing analyses) that integrates CNV and other systematic factors directly into the differential enrichment engine.

The complex relationship between DNA methylation, chromatin modification, and underlying DNA sequence is often difficult to unravel with existing technologies. Here, we describe a novel technique based on high-throughput sequencing of bisulfite-treated chromatin immunoprecipitated DNA (BisChIP-seq), which can directly interrogate genetic and epigenetic processes that occur in normal and diseased cells. Unlike most previous reports based on correlative techniques, we found using direct bisulfite sequencing of Polycomb H3K27me3-enriched DNA from normal and prostate cancer cells that DNA methylation and H3K27me3-marked histones are not always mutually exclusive, but can co-occur in a genomic region-dependent manner. Notably, in cancer, the co-dependency of marks is largely redistributed with an increase of the dual repressive marks at CpG islands and transcription start sites of silent genes. In contrast, there is a loss of DNA methylation in intergenic H3K27me3-marked regions. Allele-specific methylation status derived from the BisChIP-seq data clearly showed that both methylated and unmethylated alleles can simultaneously be associated with H3K27me3 histones, highlighting that DNA methylation status in these regions is not dependent on Polycomb chromatin status. BisChIP-seq is a novel approach that can be widely applied to directly interrogate the genomic relationship between allele-specific DNA methylation, histone modification, or other important epigenetic regulators.