Region-based epigenetic clock design improves RRBS-based age prediction

Recent studies suggest that epigenetic rejuvenation can be achieved using drugs that mimic calorie restriction and techniques such as reprogramming-induced rejuvenation. To effectively test rejuvenation in vivo, mouse models are the safest alternative. However, we have found that the recent epigenetic clocks developed for mouse reduced-representation bisulphite sequencing (RRBS) data have significantly poor performance when applied to external datasets. We show that the sites captured and the coverage of key CpGs required for age prediction vary greatly between datasets, which likely contributes to the lack of transferability in RRBS clocks. To mitigate these coverage issues in RRBS-based age prediction, we present two novel design strategies that use average methylation over large regions rather than individual CpGs, whereby regions are defined by sliding windows (e.g. 5 kb), or density-based clustering of CpGs. We observe improved correlation and error in our regional blood clocks (RegBCs) compared to published individual-CpG-based techniques when applied to external datasets. The RegBCs are also more robust when applied to low coverage data and detect a negative age acceleration in mice undergoing calorie restriction. Our RegBCs offer a proof of principle that age prediction of RRBS datasets can be improved by accounting for multiple CpGs over a region, which negates the lack of read depth currently hindering individual-CpG-based approaches.     

CLOCK,PER2 and BMAL1 DNA Methylation: Association with Obesity and Metabolic Syndrome Characteristics and Monounsaturated Fat Intake

The circadian clock system instructs 24-h rhythmicity on gene expression in essentially all cells, including adipocytes, and epigenetic mechanisms may participate in this regulation. The aim of this research was to investigate the influence of obesity and metabolic syndrome (MetS) features in clock gene methylation and the involvement of these epigenetic modifications in the outcomes. Sixty normal-weight, overweight and obese women followed a 16-weeks weight reduction program. DNA methylation levels at different CpG sites of CLOCK, BMAL1 and PER2 genes were analyzed by Sequenom's MassARRAY in white blood cells obtained before the treatment. Statistical differences between normal-weight and overweight?+?obese subjects were found in the methylation status of different CpG sites of CLOCK (CpGs 1, 5-6, 8 and 11-14) and, with lower statistical significance, in BMAL1 (CpGs 6-7, 8, 15 and 16-17). The methylation pattern of different CpG sites of the three genes showed significant associations with anthropometric parameters such as body mass index and adiposity, and with a MetS score. Moreover, the baseline methylation levels of CLOCK CpG 1 and PER2 CpGs 2-3 and 25 correlated with the magnitude of weight loss. Interestingly, the percentage of methylation of CLOCK CpGs 1 and 8 showed associations with the intake of monounsaturated and polyunsaturated fatty acids. This study demonstrates for the first time an association between methylation status of CpG sites located in clock genes (CLOCK, BMAL1 and PER2) with obesity, MetS and weight loss. Moreover, the methylation status of different CpG sites in CLOCK and PER2 could be used as biomarkers of weight-loss success, particularly CLOCK CPGs 5-6. (Author correspondence: garaulet@um.es)     

Epigenetic clock and methylation studies in elephants

Age‐associated DNA‐methylation profiles have been used successfully to develop highly accurate biomarkers of age ("epigenetic clocks") in humans, mice, dogs, and other species. Here we present epigenetic clocks for African and Asian elephants. These clocks were developed using novel DNA methylation profiles of 140 elephant blood samples of known age, at loci that are highly conserved between mammalian species, using a custom Infinium array (HorvathMammalMethylChip40). We present epigenetic clocks for Asian elephants (Elephas maximus), African elephants (Loxodonta africana), and both elephant species combined. Two additional human‐elephant clocks were constructed by combining human and elephant samples. Epigenome‐wide association studies identified elephant age‐related CpGs and their proximal genes. The products of these genes play important roles in cellular differentiation, organismal development, metabolism, and circadian rhythms. Intracellular events observed to change with age included the methylation of bivalent chromatin domains, and targets of polycomb repressive complexes. These readily available epigenetic clocks can be used for elephant conservation efforts where accurate estimates of age are needed to predict demographic trends.     

DNA Methylation Analysis of Chromosome 21 Gene Promoters at Single Base Pair and Single Allele Resolution

Differential DNA methylation is an essential epigenetic signal for gene regulation, development, and disease processes. We mapped DNA methylation patterns of 190 gene promoter regions on chromosome 21 using bisulfite conversion and subclone sequencing in five human cell types. A total of 28,626 subclones were sequenced at high accuracy using (long-read) Sanger sequencing resulting in the measurement of the DNA methylation state of 580427 CpG sites. Our results show that average DNA methylation levels are distributed bimodally with enrichment of highly methylated and unmethylated sequences, both for amplicons and individual subclones, which represent single alleles from individual cells. Within CpG-rich sequences, DNA methylation was found to be anti-correlated with CpG dinucleotide density and GC content, and methylated CpGs are more likely to be flanked by AT-rich sequences. We observed over-representation of CpG sites in distances of 9, 18, and 27 bps in highly methylated amplicons. However, DNA sequence alone is not sufficient to predict an amplicon's DNA methylation status, since 43% of all amplicons are differentially methylated between the cell types studied here. DNA methylation in promoter regions is strongly correlated with the absence of gene expression and low levels of activating epigenetic marks like H3K4 methylation and H3K9 and K14 acetylation. Utilizing the single base pair and single allele resolution of our data, we found that i) amplicons from different parts of a CpG island frequently differ in their DNA methylation level, ii) methylation levels of individual cells in one tissue are very similar, and iii) methylation patterns follow a relaxed site-specific distribution. Furthermore, iv) we identified three cases of allele-specific DNA methylation on chromosome 21. Our data shed new light on the nature of methylation patterns in human cells, the sequence dependence of DNA methylation, and its function as epigenetic signal in gene regulation. Further, we illustrate genotype–epigenotype interactions by showing novel examples of allele-specific methylation.     

Epigenetic profile of developmentally important genes in bovine oocytes

Assisted reproductive technologies are associated with an increased incidence of epigenetic aberrations, specifically in imprinted genes. Here, we used the bovine oocyte as a model to determine putative epigenetic mutations at three imprinted gene loci caused by the type of maturation, either in vitro maturation (IVM) in Tissue Culture Medium 199 (TCM) or modified synthetic oviduct fluid (mSOF) medium, or in vivo maturation. We applied a limiting dilution approach and direct bisulfite sequencing to analyze the methylation profiles of individual alleles (DNA molecules) for H19/IGF2, PEG3, and SNRPN, which are each associated with imprinting defects in humans and/or the mouse model, and are known to be differentially methylated in bovine embryos. Altogether, we obtained the methylation patterns of 203 alleles containing 4,512 CpG sites from immature oocytes, 213 alleles with 4,779 CpG sites from TCM-matured oocytes, 215 alleles/4,725 CpGs in mSOF-matured oocytes, and 78 alleles/1,672 CpGs from in vivo-matured oocytes. The total rate of individual CpGs and entire allele methylation errors did not differ significantly between the two IVM and the in vivo group, indicating that current IVM protocols have no or only marginal effects on these critical epigenetic marks. Furthermore, the mRNA expression profiles of the three imprinted genes and a panel of eight other genes indicative of oocyte competence were determined by quantitative real-time PCR. We found different mRNA expression profiles between in vivo-matured oocytes versus their in vitro-matured counterparts, suggesting an influence on regulatory mechanisms other than DNA methylation.     

DNA methylation pattern of CALCA in preterm neonates with bacterial sepsis as a putative epigenetic biomarker

Diagnosis of bacterial sepsis in preterm neonates can be difficult when using serum markers that rely on physiological changes because these changes may not necessarily be the result of bacterial infections alone. This retrospective investigation explores the potential use of the DNA methylation pattern of CpG sites in the promoter region of the calcitonin-related polypeptide α (CALCA) gene as an epigenetic biomarker for bacterial sepsis in preterm newborns. Four novel changes in the DNA methylation of eight CpG sites were detected in this gene and are present only in neonates with bacterial sepsis: (1) partial methylation at -769 CpG in gram-negative or gram-positive early onset sepsis (EOS) and late onset sepsis (LOS) episodes; (2) demethylation of 8 CpGs in gram-negative EOS followed by LOS (ELS) and in gram-negative EOS; (3) demethylation of 7 CpGs in gram-positive ELS and gram-positive EOS; (4) -771 C:G > T:A; 5′ de novo -778 CpG mutation on both alleles in EOS. These changes were not detected in birth weight and gestational age matched controls or in newborns with isolated infections. Our results indicate that the DNA methylation pattern of the promoter region of the CALCA gene varies in different types of bacterial preterm sepsis, thus suggesting a potential use as an epigenetic biomarker. A prospective confirmation of these results is essential.     

Factors affecting the persistence of drug-induced reprogramming of the cancer methylome

Aberrant DNA methylation is a critical feature of cancer. Epigenetic therapy seeks to reverse these changes to restore normal gene expression. DNA demethylating agents, including 5-aza-2′-deoxycytidine (DAC), are currently used to treat certain leukemias, and can sensitize solid tumors to chemotherapy and immunotherapy. However, it has been difficult to pin the clinical efficacy of these agents to specific demethylation events, and the factors that contribute to the durability of response remain largely unknown. Here we examined the genome-wide kinetics of DAC-induced DNA demethylation and subsequent remethylation after drug withdrawal in breast cancer cells. We find that CpGs differ in both their susceptibility to demethylation and propensity for remethylation after drug removal. DAC-induced demethylation was most apparent at CpGs with higher initial methylation levels and further from CpG islands. Once demethylated, such sites exhibited varied remethylation potentials. The most rapidly remethylating CpGs regained >75% of their starting methylation within a month of drug withdrawal. These sites had higher pretreatment methylation levels, were enriched in gene bodies, marked by H3K36me3, and tended to be methylated in normal breast cells. In contrast, a more resistant class of CpG sites failed to regain even 20% of their initial methylation after 3 months. These sites had lower pretreatment methylation levels, were within or near CpG islands, marked by H3K79me2 or H3K4me2/3, and were overrepresented in sites that become aberrantly hypermethylated in breast cancers. Thus, whereas DAC-induced demethylation affects both endogenous and aberrantly methylated sites, tumor-specific hypermethylation is more slowly regained, even as normal methylation promptly recovers. Taken together, these data suggest that the durability of DAC response is linked to its selective ability to stably reset at least a portion of the cancer methylome.     

DNA methylation and chromatin accessibility predict age in the domestic dog

Across mammals, the epigenome is highly predictive of chronological age. These “epigenetic clocks,” most of which have been built using DNA methylation (DNAm) profiles, have gained traction as biomarkers of aging and organismal health. While the ability of DNAm to predict chronological age has been repeatedly demonstrated, the ability of other epigenetic features to predict age remains unclear. Here, we use two types of epigenetic information—DNAm, and chromatin accessibility as measured by ATAC-seq—to develop age predictors in peripheral blood mononuclear cells sampled from a population of domesticated dogs. We measured DNAm and ATAC-seq profiles for 71 dogs, building separate predictive clocks from each, as well as the combined dataset. We also use fluorescence-assisted cell sorting to quantify major lymphoid populations for each sample. We found that chromatin accessibility can accurately predict chronological age (R²_ATAC = 26%), though less accurately than the DNAm clock (R²_DNAm = 33%), and the clock built from the combined datasets was comparable to both (R²_combined = 29%). We also observed various populations of CD62L+ T cells significantly correlated with dog age. Finally, we found that all three clocks selected features that were in or near at least two protein-coding genes: BAIAP2 and SCARF2, both previously implicated in processes related to cognitive or neurological impairment. Taken together, these results highlight the potential of chromatin accessibility as a complementary epigenetic resource for modeling and investigating biologic age.     

Human-specific CpG “beacons” identify loci associated with human-specific traits and disease

Regulatory change has long been hypothesized to drive the delineation of the human phenotype from other closely related primates. Here we provide evidence that CpG dinucleotides play a special role in this process. CpGs enable epigenome variability via DNA methylation, and this epigenetic mark functions as a regulatory mechanism. Therefore, species-specific CpGs may influence species-specific regulation. We report non-polymorphic species-specific CpG dinucleotides (termed “CpG beacons”) as a distinct genomic feature associated with CpG island (CGI) evolution, human traits and disease. Using an inter-primate comparison, we identified 21 extreme CpG beacon clusters (≥ 20/kb peaks, empirical p < 1.0 × 10^?3) in humans, which include associations with four monogenic developmental and neurological disease related genes (Benjamini-Hochberg corrected p = 6.03 × 10^?3). We also demonstrate that beacon-mediated CpG density gain in CGIs correlates with reduced methylation in these species in orthologous CGIs over time, via human, chimpanzee and macaque MeDIP-seq. Therefore mapping into both the genomic and epigenomic space the identified CpG beacon clusters define points of intersection where a substantial two-way interaction between genetic sequence and epigenetic state has occurred. Taken together, our data support a model for CpG beacons to contribute to CGI evolution from genesis to tissue-specific to constitutively active CGIs.     

The defining DNA methylation signature of Kabuki syndrome enables functional assessment of genetic variants of unknown clinical significance

Kabuki syndrome (KS) is caused by mutations in KMT2D, which is a histone methyltransferase involved in methylation of H3K4, a histone marker associated with DNA methylation. Analysis of >450,000 CpGs in 24 KS patients with pathogenic mutations in KMT2D and 216 controls, identified 24 genomic regions, along with 1,504 CpG sites with significant DNA methylation changes including a number of Hox genes and the MYO1F gene. Using the most differentiating and significant probes and regions we developed a “methylation variant pathogenicity (MVP) score,” which enables 100% sensitive and specific identification of individuals with KS, which was confirmed using multiple public and internal patient DNA methylation databases. We also demonstrated the ability of the MVP score to accurately reclassify variants of unknown significance in subjects with apparent clinical features of KS, enabling its potential use in molecular diagnostics. These findings provide novel insights into the molecular etiology of KS and illustrate that DNA methylation patterns can be interpreted as ‘epigenetic echoes’ in certain clinical disorders.     

Nucleotide distance influences co-methylation between nearby CpG sites

The tendency of individual CpG sites to be methylated is distinctive, non-random and well-regulated throughout the genome. We investigated the structural and spatial factors influencing CpGs methylation by performing an ultra-deep targeted methylation analysis on human, mouse and zebrafish genes. We found that methylation is not a random process and that closer neighboring CpG sites are more likely to share the same methylation status. Moreover, if the distance between CpGs increases, the degree of co-methylation decreases. We set up a simulation model to analyze the contribution of both the intrinsic susceptibility and the distance effect on the probability of a CpG to be methylated. Our finding suggests that the establishment of a specific methylation pattern follows a universal rule that must take into account of the synergistic and dynamic interplay of these two main factors: the intrinsic methylation susceptibility of specific CpG and the nucleotide distance between two CpG sites.