DNMT3B PWWP mutations cause hypermethylation of heterochromatin.

The correct establishment of DNA methylation patterns is vital for mammalian development and is achieved by the de novo DNA methyltransferases DNMT3A and DNMT3B. DNMT3B localises to H3K36me3 at actively transcribing gene bodies via its PWWP domain. It also functions at heterochromatin through an unknown recruitment mechanism. Here, we find that knockout of DNMT3B causes loss of methylation predominantly at H3K9me3-marked heterochromatin and that DNMT3B PWWP domain mutations or deletion result in striking increases of methylation in H3K9me3-marked heterochromatin. Removal of the N-terminal region of DNMT3B affects its ability to methylate H3K9me3-marked regions. This region of DNMT3B directly interacts with HP1α and facilitates the bridging of DNMT3B with H3K9me3-marked nucleosomes in vitro. Our results suggest that DNMT3B is recruited to H3K9me3-marked heterochromatin in a PWWP-independent manner that is facilitated by the protein's N-terminal region through an interaction with a key heterochromatin protein. More generally, we suggest that DNMT3B plays a role in DNA methylation homeostasis at heterochromatin, a process which is disrupted in cancer, aging and Immunodeficiency, Centromeric Instability and Facial Anomalies (ICF) syndrome.

Genome-wide single-molecule analysis of long-read DNA methylation reveals heterogeneous patterns at heterochromatin that reflect nucleosome organisation.

High-throughput sequencing technology is central to our current understanding of the human methylome. The vast majority of studies use chemical conversion to analyse bulk-level patterns of DNA methylation across the genome from a population of cells. While this technology has been used to probe single-molecule methylation patterns, such analyses are limited to short reads of a few hundred basepairs. DNA methylation can also be directly detected using Nanopore sequencing which can generate reads measuring megabases in length. However, thus far these analyses have largely focused on bulk-level assessment of DNA methylation. Here, we analyse DNA methylation in single Nanopore reads from human lymphoblastoid cells, to show that bulk-level metrics underestimate large-scale heterogeneity in the methylome. We use the correlation in methylation state between neighbouring sites to quantify single-molecule heterogeneity and find that heterogeneity varies significantly across the human genome, with some regions having heterogeneous methylation patterns at the single-molecule level and others possessing more homogeneous methylation patterns. By comparing the genomic distribution of the correlation to epigenomic annotations, we find that the greatest heterogeneity in single-molecule patterns is observed within heterochromatic partially methylated domains (PMDs). In contrast, reads originating from euchromatic regions and gene bodies have more ordered DNA methylation patterns. By analysing the patterns of single molecules in more detail, we show the existence of a nucleosome-scale periodicity in DNA methylation that accounts for some of the heterogeneity we uncover in long single-molecule DNA methylation patterns. We find that this periodic structure is partially masked in bulk data and correlates with DNA accessibility as measured by nanoNOMe-seq, suggesting that it could be generated by nucleosomes. Our findings demonstrate the power of single-molecule analysis of long-read data to understand the structure of the human methylome.

Germline de novo mutations in families with Mendelian cancer syndromes caused by defects in DNA repair.

DNA repair defects underlie many cancer syndromes. We tested whether de novo germline mutations (DNMs) are increased in families with germline defects in polymerase proofreading or base excision repair. A parent with a single germline POLE or POLD1 mutation, or biallelic MUTYH mutations, had 3-4 fold increased DNMs over sex-matched controls. POLE had the largest effect. The DNMs carried mutational signatures of the appropriate DNA repair deficiency. No DNM increase occurred in offspring of MUTYH heterozygous parents. Parental DNA repair defects caused about 20-150 DNMs per child, additional to the ~60 found in controls, but almost all extra DNMs occurred in non-coding regions. No increase in post-zygotic mutations was detected, excepting a child with bi-allelic MUTYH mutations who was excluded from the main analysis; she had received chemotherapy and may have undergone oligoclonal haematopoiesis. Inherited DNA repair defects associated with base pair-level mutations increase DNMs, but phenotypic consequences appear unlikely.

Local CpG density affects the trajectory and variance of age-associated DNA methylation changes.

BACKGROUND: DNA methylation is an epigenetic mark associated with the repression of gene promoters. Its pattern in the genome is disrupted with age and these changes can be used to statistically predict age with epigenetic clocks. Altered rates of aging inferred from these clocks are observed in human disease. However, the molecular mechanisms underpinning age-associated DNA methylation changes remain unknown. Local DNA sequence can program steady-state DNA methylation levels, but how it influences age-associated methylation changes is unknown. RESULTS: We analyze longitudinal human DNA methylation trajectories at 345,895 CpGs from 600 individuals aged between 67 and 80 to understand the factors responsible for age-associated epigenetic changes at individual CpGs. We show that changes in methylation with age occur at 182,760 loci largely independently of variation in cell type proportions. These changes are especially apparent at 8322 low CpG density loci. Using SNP data from the same individuals, we demonstrate that methylation trajectories are affected by local sequence polymorphisms at 1487 low CpG density loci. More generally, we find that low CpG density regions are particularly prone to change and do so variably between individuals in people aged over 65. This differs from the behavior of these regions in younger individuals where they predominantly lose methylation. CONCLUSIONS: Our results, which we reproduce in two independent groups of individuals, demonstrate that local DNA sequence influences age-associated DNA methylation changes in humans in vivo. We suggest that this occurs because interactions between CpGs reinforce maintenance of methylation patterns in CpG dense regions.

In Vivo Modeling of Patient Genetic Heterogeneity Identifies New Ways to Target Cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (ICC) is an aggressive malignancy of the bile ducts within the liver characterized by high levels of genetic heterogeneity. In the context of such genetic variability, determining which oncogenic mutations drive ICC growth has been difficult, and developing modes of patient stratification and targeted therapies remains challenging. Here we model the interactions between rare mutations with more common driver genes and combine in silico analysis of patient data with highly multiplexed in vivo CRISPR-spCas9 screens to perform a functional in vivo study into the role genetic heterogeneity plays in driving ICC. Novel tumor suppressors were uncovered, which, when lost, cooperate with the RAS oncoprotein to drive ICC growth. Focusing on a set of driver mutations that interact with KRAS to initiate aggressive, sarcomatoid-type ICC revealed that tumor growth relies on Wnt and PI3K signaling. Pharmacologic coinhibition of Wnt and PI3K in vivo impeded ICC growth regardless of mutational profile. Therefore, Wnt and PI3K activity should be considered as a signature by which patients can be stratified for treatment independent of tumor genotype, and inhibitors of these pathways should be levied to treat ICC. SIGNIFICANCE: This work shows that, despite significant genetic heterogeneity, intrahepatic cholangiocarcinoma relies on a limited number of signaling pathways to grow, suggesting common therapeutic vulnerabilities across patients.

Cluster mean-field theory accurately predicts statistical properties of large-scale DNA methylation patterns.

The accurate establishment and maintenance of DNA methylation patterns is vital for mammalian development and disruption to these processes causes human disease. Our understanding of DNA methylation mechanisms has been facilitated by mathematical modelling, particularly stochastic simulations. Megabase-scale variation in DNA methylation patterns is observed in development, cancer and ageing and the mechanisms generating these patterns are little understood. However, the computational cost of stochastic simulations prevents them from modelling such large genomic regions. Here, we test the utility of three different mean-field models to predict summary statistics associated with large-scale DNA methylation patterns. By comparison to stochastic simulations, we show that a cluster mean-field model accurately predicts the statistical properties of steady-state DNA methylation patterns, including the mean and variance of methylation levels calculated across a system of CpG sites, as well as the covariance and correlation of methylation levels between neighbouring sites. We also demonstrate that a cluster mean-field model can be used within an approximate Bayesian computation framework to accurately infer model parameters from data. As mean-field models can be solved numerically in a few seconds, our work demonstrates their utility for understanding the processes underpinning large-scale DNA methylation patterns.

Glioblastomas acquire myeloid-affiliated transcriptional programs via epigenetic immunoediting to elicit immune evasion.

Glioblastoma multiforme (GBM) is an aggressive brain tumor for which current immunotherapy approaches have been unsuccessful. Here, we explore the mechanisms underlying immune evasion in GBM. By serially transplanting GBM stem cells (GSCs) into immunocompetent hosts, we uncover an acquired capability of GSCs to escape immune clearance by establishing an enhanced immunosuppressive tumor microenvironment. Mechanistically, this is not elicited via genetic selection of tumor subclones, but through an epigenetic immunoediting process wherein stable transcriptional and epigenetic changes in GSCs are enforced following immune attack. These changes launch a myeloid-affiliated transcriptional program, which leads to increased recruitment of tumor-associated macrophages. Furthermore, we identify similar epigenetic and transcriptional signatures in human mesenchymal subtype GSCs. We conclude that epigenetic immunoediting may drive an acquired immune evasion program in the most aggressive mesenchymal GBM subtype by reshaping the tumor immune microenvironment.

De novo DNA methyltransferase activity in colorectal cancer is directed towards H3K36me3 marked CpG islands.

The aberrant gain of DNA methylation at CpG islands is frequently observed in colorectal tumours and may silence the expression of tumour suppressors such as MLH1. Current models propose that these CpG islands are targeted by de novo DNA methyltransferases in a sequence-specific manner, but this has not been tested. Using ectopically integrated CpG islands, here we find that aberrantly methylated CpG islands are subject to low levels of de novo DNA methylation activity in colorectal cancer cells. By delineating DNA methyltransferase targets, we find that instead de novo DNA methylation activity is targeted primarily to CpG islands marked by the histone modification H3K36me3, a mark associated with transcriptional elongation. These H3K36me3 marked CpG islands are heavily methylated in colorectal tumours and the normal colon suggesting that de novo DNA methyltransferase activity at CpG islands in colorectal cancer is focused on similar targets to normal tissues and not greatly remodelled by tumourigenesis.

An epigenetic predictor of death captures multi-modal measures of brain health.

Individuals of the same chronological age exhibit disparate rates of biological ageing. Consequently, a number of methodologies have been proposed to determine biological age and primarily exploit variation at the level of DNA methylation (DNAm). A novel epigenetic clock, termed 'DNAm GrimAge' has outperformed its predecessors in predicting the risk of mortality as well as many age-related morbidities. However, the association between DNAm GrimAge and cognitive or neuroimaging phenotypes remains unknown. We explore these associations in the Lothian Birth Cohort 1936 (n = 709, mean age 73 years). Higher DNAm GrimAge was strongly associated with all-cause mortality over the eighth decade (Hazard Ratio per standard deviation increase in GrimAge: 1.81, P < 2.0 × 10-16). Higher DNAm GrimAge was associated with lower age 11 IQ (ß = -0.11), lower age 73 general cognitive ability (ß = -0.18), decreased brain volume (ß = -0.25) and increased brain white matter hyperintensities (ß = 0.17). There was tentative evidence for a longitudinal association between DNAm GrimAge and cognitive decline from age 70 to 79. Sixty-nine of 137 health- and brain-related phenotypes tested were significantly associated with GrimAge. Adjusting all models for childhood intelligence attenuated to non-significance a small number of associations (12/69 associations; 6 of which were cognitive traits), but not the association with general cognitive ability (33.9% attenuation). Higher DNAm GrimAge associates with lower cognitive ability and brain vascular lesions in older age, independently of early-life cognitive ability. This epigenetic predictor of mortality associates with different measures of brain health and may aid in the prediction of age-related cognitive decline.

cGAS-mediated induction of type I interferon due to inborn errors of histone pre-mRNA processing.

Inappropriate stimulation or defective negative regulation of the type I interferon response can lead to autoinflammation. In genetically uncharacterized cases of the type I interferonopathy Aicardi-Goutières syndrome, we identified biallelic mutations in LSM11 and RNU7-1, which encode components of the replication-dependent histone pre-mRNA-processing complex. Mutations were associated with the misprocessing of canonical histone transcripts and a disturbance of linker histone stoichiometry. Additionally, we observed an altered distribution of nuclear cyclic guanosine monophosphate-adenosine monophosphate synthase (cGAS) and enhanced interferon signaling mediated by the cGAS-stimulator of interferon genes (STING) pathway in patient-derived fibroblasts. Finally, we established that chromatin without linker histone stimulates cyclic guanosine monophosphate-adenosine monophosphate (cGAMP) production in vitro more efficiently. We conclude that nuclear histones, as key constituents of chromatin, are essential in suppressing the immunogenicity of self-DNA.

Transition to naïve human pluripotency mirrors pan-cancer DNA hypermethylation.

Epigenetic reprogramming is a cancer hallmark, but how it unfolds during early neoplastic events and its role in carcinogenesis and cancer progression is not fully understood. Here we show that resetting from primed to naïve human pluripotency results in acquisition of a DNA methylation landscape mirroring the cancer DNA methylome, with gradual hypermethylation of bivalent developmental genes. We identify a dichotomy between bivalent genes that do and do not become hypermethylated, which is also mirrored in cancer. We find that loss of H3K4me3 at bivalent regions is associated with gain of methylation. Additionally, we observe that promoter CpG island hypermethylation is not restricted solely to emerging naïve cells, suggesting that it is a feature of a heterogeneous intermediate population during resetting. These results indicate that transition to naïve pluripotency and oncogenic transformation share common epigenetic trajectories, which implicates reprogramming and the pluripotency network as a central hub in cancer formation.

Activation of transcription factor circuity in 2i-induced ground state pluripotency is independent of repressive global epigenetic landscapes.

Mouse embryonic stem cells (mESCs) cultured with MEK/ERK and GSK3ß (2i) inhibitors transition to ground state pluripotency. Gene expression changes, redistribution of histone H3K27me3 profiles and global DNA hypomethylation are hallmarks of 2i exposure, but it is unclear whether epigenetic alterations are required to achieve and maintain ground state or occur as an outcome of 2i signal induced changes. Here we show that ESCs with three epitypes, WT, constitutively methylated, or hypomethylated, all undergo comparable morphological, protein expression and transcriptome changes independently of global alterations of DNA methylation levels or changes in H3K27me3 profiles. Dazl and Fkbp6 expression are induced by 2i in all three epitypes, despite exhibiting hypermethylated promoters in constitutively methylated ESCs. We identify a number of activated gene promoters that undergo 2i dependent loss of H3K27me3 in all three epitypes, however genetic and pharmaceutical inhibition experiments show that H3K27me3 is not required for their silencing in non-2i conditions. By separating and defining their contributions, our data suggest that repressive epigenetic systems play minor roles in mESC self-renewal and naïve ground state establishment by core sets of dominant pluripotency associated transcription factor networks, which operate independently from these epigenetic processes.

Age-related clonal haemopoiesis is associated with increased epigenetic age.

Age-related clonal haemopoiesis (ARCH) in healthy individuals was initially observed through an increased skewing in X-chromosome inactivation [1]. More recently, several groups reported that ARCH is driven by somatic mutations [2], with the most prevalent ARCH mutations being in the DNMT3A and TET2 genes, previously described as drivers of myeloid malignancies. ARCH is associated with an increased risk for haematological cancers [2]. ARCH also confers an increased risk for non-haematological diseases, such as cardiovascular disease, atherosclerosis, and chronic ischemic heart failure, for which age is a main risk factor [3,4]. Whether ARCH is linked to accelerated ageing has remained unexplored. The most accurate and commonly used tools to measure age acceleration are epigenetic clocks: they are based on age-related methylation differences at specific CpG sites [5]. Deviations from chronological age towards an increased epigenetic age have been associated with increased risk of earlier mortality and age-related morbidities [5,6]. Here we present evidence of accelerated epigenetic age in individuals with ARCH.

Publisher Correction: Parent of origin genetic effects on methylation in humans are common and influence complex trait variation.

In the original version of this Article, the legend in the upper panel of Figure 2 incorrectly read 'paternal imprinting' and should have read 'maternal imprinting'. This has been corrected in both the PDF and HTML versions of the Article.

Parent of origin genetic effects on methylation in humans are common and influence complex trait variation.

Parent-of-origin effects (POE) exist when there is differential expression of alleles inherited from the two parents. A genome-wide scan for POE on DNA methylation at 639,238 CpGs in 5,101 individuals identifies 733 independent methylation CpGs potentially influenced by POE at a false discovery rate ≤ 0.05 of which 331 had not previously been identified. Cis and trans methylation quantitative trait loci (mQTL) regulate methylation variation through POE at 54% (399/733) of the identified POE-influenced CpGs. The combined results provide strong evidence for previously unidentified POE-influenced CpGs at 171 independent loci. Methylation variation at 14 of the POE-influenced CpGs is associated with multiple metabolic traits. A phenome-wide association analysis using the POE mQTL SNPs identifies a previously unidentified imprinted locus associated with waist circumference. These results provide a high resolution population-level map for POE on DNA methylation sites, their local and distant regulators and potential consequences for complex traits.

Gain-of-function DNMT3A mutations cause microcephalic dwarfism and hypermethylation of Polycomb-regulated regions.

DNA methylation and Polycomb are key factors in the establishment of vertebrate cellular identity and fate. Here we report de novo missense mutations in DNMT3A, which encodes the DNA methyltransferase DNMT3A. These mutations cause microcephalic dwarfism, a hypocellular disorder of extreme global growth failure. Substitutions in the PWWP domain abrogate binding to the histone modifications H3K36me2 and H3K36me3, and alter DNA methylation in patient cells. Polycomb-associated DNA methylation valleys, hypomethylated domains encompassing developmental genes, become methylated with concomitant depletion of H3K27me3 and H3K4me3 bivalent marks. Such de novo DNA methylation occurs during differentiation of Dnmt3aW326R pluripotent cells in vitro, and is also evident in Dnmt3aW326R/+ dwarf mice. We therefore propose that the interaction of the DNMT3A PWWP domain with H3K36me2 and H3K36me3 normally limits DNA methylation of Polycomb-marked regions. Our findings implicate the interplay between DNA methylation and Polycomb at key developmental regulators as a determinant of organism size in mammals.

Genotype effects contribute to variation in longitudinal methylome patterns in older people.

BACKGROUND: DNA methylation levels change along with age, but few studies have examined the variation in the rate of such changes between individuals. METHODS: We performed a longitudinal analysis to quantify the variation in the rate of change of DNA methylation between individuals using whole blood DNA methylation array profiles collected at 2-4 time points (N = 2894) in 954 individuals (67-90 years). RESULTS: After stringent quality control, we identified 1507 DNA methylation CpG sites (rsCpGs) with statistically significant variation in the rate of change (random slope) of DNA methylation among individuals in a mixed linear model analysis. Genes in the vicinity of these rsCpGs were found to be enriched in Homeobox transcription factors and the Wnt signalling pathway, both of which are related to ageing processes. Furthermore, we investigated the SNP effect on the random slope. We found that 4 out of 1507 rsCpGs had one significant (P < 5 × 10-8/1507) SNP effect and 343 rsCpGs had at least one SNP effect (436 SNP-probe pairs) reaching genome-wide significance (P < 5 × 10-8). Ninety-five percent of the significant (P < 5 × 10-8) SNPs are on different chromosomes from their corresponding probes. CONCLUSIONS: We identified CpG sites that have variability in the rate of change of DNA methylation between individuals, and our results suggest a genetic basis of this variation. Genes around these CpG sites have been reported to be involved in the ageing process.

Variable outcome and methylation status according to CEBPA mutant type in double-mutated acute myeloid leukemia patients and the possible implications for treatment.

Although CEBPA double-mutated (CEBPADM) acute myeloid leukemia is considered to be a favorable-risk disease, relapse remains a major cause of treatment failure. Most CEBPADM patients have a classic biallelic mutant combination with an N-terminal mutation leading to production of p30 protein plus a C-terminal loss-of-function in-frame indel mutation (CEBPAClassic-DM), but approximately one-third of cases have one or more non-classic mutations, with diverse combinations reported, and there is little information on the consequences of such mutants. We evaluated outcome in a cohort of 104 CEBPADM patients, 79 CEBPAClassic-DM and 25 with non-classic mutants, and found that the latter may have poorer survival (5-year overall survival 64% vs. 46%; P=0.05), particularly post relapse (41% vs. 0%; P=0.02). However, for this analysis, all non-classic cases were grouped together, irrespective of mutant combination. As CEBPADM cases have been reported to be hypermethylated, we used methylation profiling to assess whether this could segregate the different mutants. We developed a CEBPAClassic-DM methylation signature from a preliminary cohort of 10 CEBPADM (including 8 CEBPAClassic-DM) and 30 CEBPA wild-type (CEBPAWT) samples, and independently validated the signature in 17 CEBPAClassic-DM cases. Assessment of the signature in 16 CEBPADM cases with different non-classic mutant combinations showed that only 31% had a methylation profile equivalent to CEBPAClassic-DM whereas for 69% the profile was either intermediate between CEBPAClassic-DM and CEBPAWT or equivalent to CEBPAWT These results suggest that CEBPADM cases with non-classic mutants may be functionally different from those with CEBPAClassic-DM mutants, and should not automatically be included in the same prognostic group. (AML12 is registered under ISRCTN17833622 and AML15 under ISRCTN17161961).

Elevated FOXG1 and SOX2 in glioblastoma enforces neural stem cell identity through transcriptional control of cell cycle and epigenetic regulators.

Glioblastoma multiforme (GBM) is an aggressive brain tumor driven by cells with hallmarks of neural stem (NS) cells. GBM stem cells frequently express high levels of the transcription factors FOXG1 and SOX2. Here we show that increased expression of these factors restricts astrocyte differentiation and can trigger dedifferentiation to a proliferative NS cell state. Transcriptional targets include cell cycle and epigenetic regulators (e.g., Foxo3, Plk1, Mycn, Dnmt1, Dnmt3b, and Tet3). Foxo3 is a critical repressed downstream effector that is controlled via a conserved FOXG1/SOX2-bound cis-regulatory element. Foxo3 loss, combined with exposure to the DNA methylation inhibitor 5-azacytidine, enforces astrocyte dedifferentiation. DNA methylation profiling in differentiating astrocytes identifies changes at multiple polycomb targets, including the promoter of Foxo3 In patient-derived GBM stem cells, CRISPR/Cas9 deletion of FOXG1 does not impact proliferation in vitro; however, upon transplantation in vivo, FOXG1-null cells display increased astrocyte differentiation and up-regulate FOXO3. In contrast, SOX2 ablation attenuates proliferation, and mutant cells cannot be expanded in vitro. Thus, FOXG1 and SOX2 operate in complementary but distinct roles to fuel unconstrained self-renewal in GBM stem cells via transcriptional control of core cell cycle and epigenetic regulators.