Machine learning enables pan-cancer identification of mutational hotspots at persistent CTCF binding sites.

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

The H2A.Z-nuclesome code in mammals: emerging functions.

H2A.Z is a histone variant that provides specific structural and docking-side properties to the nucleosome, resulting in diverse and specialised molecular and cellular functions. In this review, we discuss the latest studies uncovering new functional aspects of mammalian H2A.Z in gene transcription, including pausing and elongation of RNA polymerase II (RNAPII) and enhancer activity; DNA repair; DNA replication; and 3D chromatin structure. We also review the recently described role of H2A.Z in embryonic development, cell differentiation, neurodevelopment, and brain function. In conclusion, our cumulative knowledge of H2A.Z over the past 40 years, in combination with the implementation of novel molecular technologies, is unravelling an unexpected and complex role of histone variants in gene regulation and disease.

Detailed DNA methylation characterisation of phyllodes tumours identifies a signature of malignancy and distinguishes phyllodes from metaplastic breast carcinoma.

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.

Characterisation and reproducibility of the HumanMethylationEPIC v2.0 BeadChip for DNA methylation profiling.

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.

Transient exposure to miR-203 enhances the differentiation capacity of established pluripotent stem cells.

Full differentiation potential along with self-renewal capacity is a major property of pluripotent stem cells (PSCs). However, the differentiation capacity frequently decreases during expansion of PSCs in vitro. We show here that transient exposure to a single microRNA, expressed at early stages during normal development, improves the differentiation capacity of already-established murine and human PSCs. Short exposure to miR-203 in PSCs (miPSCs) induces a transient expression of 2C markers that later results in expanded differentiation potency to multiple lineages, as well as improved efficiency in tetraploid complementation and human-mouse interspecies chimerism assays. Mechanistically, these effects are at least partially mediated by direct repression of de novo DNA methyltransferases Dnmt3a and Dnmt3b, leading to transient and reversible erasure of DNA methylation. These data support the use of transient exposure to miR-203 as a versatile method to reset the epigenetic memory in PSCs, and improve their effectiveness in regenerative medicine.

Calling differentially methylated regions from whole genome bisulphite sequencing with DMRcate.

Whole genome bisulphite sequencing (WGBS) permits the genome-wide study of single molecule methylation patterns. One of the key goals of mammalian cell-type identity studies, in both normal differentiation and disease, is to locate differential methylation patterns across the genome. We discuss the most desirable characteristics for DML (differentially methylated locus) and DMR (differentially methylated region) detection tools in a genome-wide context and choose a set of statistical methods that fully or partially satisfy these considerations to compare for benchmarking. Our data simulation strategy is both biologically informed-employing distribution parameters derived from large-scale consortium datasets-and thorough. We report DML detection ability with respect to coverage, group methylation difference, sample size, variability and covariate size, both marginally and jointly, and exhaustively with respect to parameter combination. We also benchmark these methods on FDR control and computational time. We use this result to backend and introduce an expanded version of DMRcate: an existing DMR detection tool for microarray data that we have extended to now call DMRs from WGBS data. We compare DMRcate to a set of alternative DMR callers using a similarly realistic simulation strategy. We find DMRcate and RADmeth are the best predictors of DMRs, and conclusively find DMRcate the fastest.

ELF5 modulates the estrogen receptor cistrome in breast cancer.

Acquired resistance to endocrine therapy is responsible for half of the therapeutic failures in the treatment of breast cancer. Recent findings have implicated increased expression of the ETS transcription factor ELF5 as a potential modulator of estrogen action and driver of endocrine resistance, and here we provide the first insight into the mechanisms by which ELF5 modulates estrogen sensitivity. Using chromatin immunoprecipitation sequencing we found that ELF5 binding overlapped with FOXA1 and ER at super enhancers, enhancers and promoters, and when elevated, caused FOXA1 and ER to bind to new regions of the genome, in a pattern that replicated the alterations to the ER/FOXA1 cistrome caused by the acquisition of resistance to endocrine therapy. RNA sequencing demonstrated that these changes altered estrogen-driven patterns of gene expression, the expression of ER transcription-complex members, and 6 genes known to be involved in driving the acquisition of endocrine resistance. Using rapid immunoprecipitation mass spectrometry of endogenous proteins, and proximity ligation assays, we found that ELF5 interacted physically with members of the ER transcription complex, such as DNA-PKcs. We found 2 cases of endocrine-resistant brain metastases where ELF5 levels were greatly increased and ELF5 patterns of gene expression were enriched, compared to the matched primary tumour. Thus ELF5 alters ER-driven gene expression by modulating the ER/FOXA1 cistrome, by interacting with it, and by modulating the expression of members of the ER transcriptional complex, providing multiple mechanisms by which ELF5 can drive endocrine resistance.

MethPanel: a parallel pipeline and interactive analysis tool for multiplex bisulphite PCR sequencing to assess DNA methylation biomarker panels for disease detection.

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.

Macrophage development and activation involve coordinated intron retention in key inflammatory regulators.

Monocytes and macrophages are essential components of the innate immune system. Herein, we report that intron retention (IR) plays an important role in the development and function of these cells. Using Illumina mRNA sequencing, Nanopore direct cDNA sequencing and proteomics analysis, we identify IR events that affect the expression of key genes/proteins involved in macrophage development and function. We demonstrate that decreased IR in nuclear-detained mRNA is coupled with increased expression of genes encoding regulators of macrophage transcription, phagocytosis and inflammatory signalling, including ID2, IRF7, ENG and LAT. We further show that this dynamic IR program persists during the polarisation of resting macrophages into activated macrophages. In the presence of proinflammatory stimuli, intron-retaining CXCL2 and NFKBIZ transcripts are rapidly spliced, enabling timely expression of these key inflammatory regulators by macrophages. Our study provides novel insights into the molecular factors controlling vital regulators of the innate immune response.

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.

Mother-child transmission of epigenetic information by tunable polymorphic imprinting.

Genomic imprinting mediated by DNA methylation restricts gene expression to a single allele determined by parental origin and is not generally considered to be under genetic or environmental influence. Here, we focused on a differentially methylated region (DMR) of approximately 1.9 kb that includes a 101-bp noncoding RNA gene (nc886/VTRNA2-1), which is maternally imprinted in ∼75% of humans. This is unlike other imprinted genes, which demonstrate monoallelic methylation in 100% of individuals. The DMR includes a CTCF binding site on the centromeric side defining the DMR boundary and is flanked by a CTCF binding site on the telomeric side. The centromeric CTCF binding site contains an A/C polymorphism (rs2346018); the C allele is associated with less imprinting. The frequency of imprinting of the nc886 DMR in infants was linked to at least two nongenetic factors, maternal age at delivery and season of conception. In a separate cohort, nc886 imprinting was associated with lower body mass index in children at 5 y of age. Thus, we propose that the imprinting status of the nc886 DMR is "tunable" in that it is associated with maternal haplotype and prenatal environment. This provides a potential mechanism for transmitting information, with phenotypic consequences, from mother to child.

Multiple links between 5-methylcytosine content of mRNA and translation.

BACKGROUND: 5-Methylcytosine (m5C) is a prevalent base modification in tRNA and rRNA but it also occurs more broadly in the transcriptome, including in mRNA, where it serves incompletely understood molecular functions. In pursuit of potential links of m5C with mRNA translation, we performed polysome profiling of human HeLa cell lysates and subjected RNA from resultant fractions to efficient bisulfite conversion followed by RNA sequencing (bsRNA-seq). Bioinformatic filters for rigorous site calling were devised to reduce technical noise. RESULTS: We obtained ~ 1000 candidate m5C sites in the wider transcriptome, most of which were found in mRNA. Multiple novel sites were validated by amplicon-specific bsRNA-seq in independent samples of either human HeLa, LNCaP and PrEC cells. Furthermore, RNAi-mediated depletion of either the NSUN2 or TRDMT1 m5C:RNA methyltransferases showed a clear dependence on NSUN2 for the majority of tested sites in both mRNAs and noncoding RNAs. Candidate m5C sites in mRNAs are enriched in 5'UTRs and near start codons and are embedded in a local context reminiscent of the NSUN2-dependent m5C sites found in the variable loop of tRNA. Analysing mRNA sites across the polysome profile revealed that modification levels, at bulk and for many individual sites, were inversely correlated with ribosome association. CONCLUSIONS: Our findings emphasise the major role of NSUN2 in placing the m5C mark transcriptome-wide. We further present evidence that substantiates a functional interdependence of cytosine methylation level with mRNA translation. Additionally, we identify several compelling candidate sites for future mechanistic analysis.

Evaluation of cross-platform and interlaboratory concordance via consensus modelling of genomic measurements.

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.

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.

Epigenomics of mammary gland development.

Differentiation of stem cells into highly specialised cells requires gene expression changes brought about by remodelling of the chromatin architecture. During this lineage-commitment process, the majority of DNA needs to be packaged into inactive heterochromatin, allowing only a subset of regulatory elements to remain open and functionally required genes to be expressed. Epigenetic mechanisms such as DNA methylation, post-translational modifications to histone tails, and nucleosome positioning all potentially contribute to the changes in higher order chromatin structure during differentiation. The mammary gland is a particularly useful model to study these complex epigenetic processes since the majority of its development is postnatal, the gland is easily accessible, and development occurs in a highly reproducible manner. Inappropriate epigenetic remodelling can also drive tumourigenesis; thus, insights into epigenetic remodelling during mammary gland development advance our understanding of breast cancer aetiology. We review the current literature surrounding DNA methylation and histone modifications in the developing mammary gland and its implications for breast cancer.

Long-range regulators of the lncRNA HOTAIR enhance its prognostic potential in breast cancer.

Predicting response to endocrine therapy and survival in oestrogen receptor positive breast cancer is a significant clinical challenge and novel prognostic biomarkers are needed. Long-range regulators of gene expression are emerging as promising biomarkers and therapeutic targets for human diseases, so we have explored the potential of distal enhancer elements of non-coding RNAs in the prognostication of breast cancer survival. HOTAIR is a long non-coding RNA that is overexpressed, promotes metastasis and is predictive of decreased survival. Here, we describe a long-range transcriptional enhancer of the HOTAIR gene that binds several hormone receptors and associated transcription factors, interacts with the HOTAIR promoter and augments transcription. This enhancer is dependent on Forkhead-Box transcription factors and functionally interacts with a novel alternate HOTAIR promoter. HOTAIR expression is negatively regulated by oestrogen, positively regulated by FOXA1 and FOXM1, and is inversely correlated with oestrogen receptor and directly correlated with FOXM1 in breast tumours. The combination of HOTAIR and FOXM1 enables greater discrimination of endocrine therapy responders and non-responders in patients with oestrogen receptor positive breast cancer. Consistent with this, HOTAIR expression is increased in cell-line models of endocrine resistance. Analysis of breast cancer gene expression data indicates that HOTAIR is co-expressed with FOXA1 and FOXM1 in HER2-enriched tumours, and these factors enhance the prognostic power of HOTAIR in aggressive HER2+ breast tumours. Our study elucidates the transcriptional regulation of HOTAIR, identifies HOTAIR and its regulators as novel biomarkers of patient response to endocrine therapy and corroborates the importance of transcriptional enhancers in cancer.

ELF5 Drives Lung Metastasis in Luminal Breast Cancer through Recruitment of Gr1+ CD11b+ Myeloid-Derived Suppressor Cells.

During pregnancy, the ETS transcription factor ELF5 establishes the milk-secreting alveolar cell lineage by driving a cell fate decision of the mammary luminal progenitor cell. In breast cancer, ELF5 is a key transcriptional determinant of tumor subtype and has been implicated in the development of insensitivity to anti-estrogen therapy. In the mouse mammary tumor virus-Polyoma Middle T (MMTV-PyMT) model of luminal breast cancer, induction of ELF5 levels increased leukocyte infiltration, angiogenesis, and blood vessel permeability in primary tumors and greatly increased the size and number of lung metastasis. Myeloid-derived suppressor cells, a group of immature neutrophils recently identified as mediators of vasculogenesis and metastasis, were recruited to the tumor in response to ELF5. Depletion of these cells using specific Ly6G antibodies prevented ELF5 from driving vasculogenesis and metastasis. Expression signatures in luminal A breast cancers indicated that increased myeloid cell invasion and inflammation were correlated with ELF5 expression, and increased ELF5 immunohistochemical staining predicted much shorter metastasis-free and overall survival of luminal A patients, defining a group who experienced unexpectedly early disease progression. Thus, in the MMTV-PyMT mouse mammary model, increased ELF5 levels drive metastasis by co-opting the innate immune system. As ELF5 has been previously implicated in the development of antiestrogen resistance, this finding implicates ELF5 as a defining factor in the acquisition of the key aspects of the lethal phenotype in luminal A breast cancer.

Mining cancer methylomes: prospects and challenges.

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.

Reconfiguration of nucleosome-depleted regions at distal regulatory elements accompanies DNA methylation of enhancers and insulators in cancer.

It is well established that cancer-associated epigenetic repression occurs concomitant with CpG island hypermethylation and loss of nucleosomes at promoters, but the role of nucleosome occupancy and epigenetic reprogramming at distal regulatory elements in cancer is still poorly understood. Here, we evaluate the scope of global epigenetic alterations at enhancers and insulator elements in prostate and breast cancer cells using simultaneous genome-wide mapping of DNA methylation and nucleosome occupancy (NOMe-seq). We find that the genomic location of nucleosome-depleted regions (NDRs) is mostly cell type specific and preferentially found at enhancers in normal cells. In cancer cells, however, we observe a global reconfiguration of NDRs at distal regulatory elements coupled with a substantial reorganization of the cancer methylome. Aberrant acquisition of nucleosomes at enhancer-associated NDRs is associated with hypermethylation and epigenetic silencing marks, and conversely, loss of nucleosomes with demethylation and epigenetic activation. Remarkably, we show that nucleosomes remain strongly organized and phased at many facultative distal regulatory elements, even in the absence of a NDR as an anchor. Finally, we find that key transcription factor (TF) binding sites also show extensive peripheral nucleosome phasing, suggesting the potential for TFs to organize NDRs genome-wide and contribute to deregulation of cancer epigenomes. Together, our findings suggest that "decommissioning" of NDRs and TFs at distal regulatory elements in cancer cells is accompanied by DNA hypermethylation susceptibility of enhancers and insulator elements, which in turn may contribute to an altered genome-wide architecture and epigenetic deregulation in malignancy.