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.

The correct establishment of DNA methylation patterns is vital for mammalian development and is achieved largely by the de novo DNA methyltransferases DNMT3A and DNMT3B. Mutations in DNMT3B can cause immunodeficiency-centromeric instability-facial anomalies type 1 (ICF1) syndrome which is characterised by hypomethylated heterochromatin. However, in the genome, DNMT3B primarily localises to actively transcribing gene bodies through the interaction of its PWWP domain with the histone modification H3K36me3 and it is unclear how it is recruited to heterochromatin. Here we show that in DNMT3B knockout cells, loss of DNA methylation predominantly occurs in heterochromatic domains marked by H3K9me3. We also find that PWWP domain mutations which disrupt DNMT3B's interaction with H3K36me3 result in striking increases of DNA methylation in H3K9me3-marked heterochromatin. Gains of methylation are also observed when the PWWP domain of DNMT3B is deleted. In contrast, we find that the ICF1 syndrome-causing PWWP mutation, S270P, does not result in hypermethylation of heterochromatin and destabilises the protein. We also show that removal of the N-terminus region of DNMT3B affects its recruitment to chromatin and ability to methylate H3K9me3 marked regions. Our results suggest that DNMT3B is recruited to H3K9me3 marked heterochromatin in a PWWP-independent manner and that this recruitment is facilitated by the protein's N-terminus. More generally, we suggest that DNMT3B plays a role in DNA methylation homeostasis at heterochromatin, a process which is disrupted in ICF syndrome, cancer and aging.

We investigate an infinite, linear system of ordinary differential equations that models the evolution of fragmenting clusters. We assume that each cluster is composed of identical units (monomers), and we allow mass to be lost, gained or conserved during each fragmentation event. By formulating the initial-value problem for the system as an abstract Cauchy problem (ACP), posed in an appropriate weighted 1 space, and then applying perturbation results from the theory of operator semigroups, we prove the existence and uniqueness of physically relevant, classical solutions for a wide class of initial cluster distributions. Additionally, we establish that it is always possible to identify a weighted 1 space on which the fragmentation semigroup is analytic, which immediately implies that the corresponding ACP is well posed for any initial distribution belonging to this particular space. We also investigate the asymptotic behaviour of solutions and show that, under appropriate restrictions on the fragmentation coefficients, solutions display the expected long-term behaviour of converging to a purely monomeric steady state. Moreover, when the fragmentation semigroup is analytic, solutions are shown to decay to this steady state at an explicitly defined exponential rate.Mathematics Subject Classification: 47D06, 34G10, 80A30, 34D05