CpG Density Research Articles

2 The landscape of accessible chromatin in bovine oocytes and early embryos

Chromatin reorganization governs gene expression regulation during pre-implantation development. However, the global chromatin landscape and its dynamics in this period remain unexplored in bovine. In this study, we constructed a genome-wide map of accessible chromatin in bovine oocytes and early embryos using an improved assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq). We analysed pools of 20 germinal vesicles or MII oocytes or 2-, 4-, 8-, 16-cell, morula, and blastocyst stage invitro-produced embryos. We conducted ATAC-seq on six pools for each stage and an additional four pools of invivo-derived morula and blastocysts and six replicates using individual Day 14 elongating embryos. We obtained ~110 million paired end reads uniquely mapped to the bovine reference genome for each stage. Hierarchical clustering, t-distributed stochastic neighbour embedding, and principal component analysis showed four distinct patterns for open chromatin status: (1) low accessibility in germinal vesicles and MII oocytes and in 2- and 4-cell embryos; (2) significantly elevated accessibility in 8-cell, 16-cell, and morula embryos; (3) less accessibility in blastocysts; and (4) extremely high accessibility in elongating embryos. This dynamic and sequential chromatin remodelling is consistent with transcription activation during the bovine minor embryonic genome activation from fertilization to 4-cell, major embryonic genome activation at 8-cell, first differentiation at blastocyst and drastic transcription initiation for embryo elongation. Genome-wide characteristics of accessible chromatin showed (1) accessible chromatin near the transcription start sites of active genes and CpG-rich promoters; (2) widespread accessible chromatin regions extensively overlapped with transposable elements; (3) distal peaks preferentially enriched for repeats including LINE, SINE, and LTR from 8-cell to morula embryos, especially for LTR, whereas enrichment in simple repeats were found from oocytes to 4-cell and in elongating embryos; and (4) highly stage-specific transcription factor motifs in distal peaks were unveiled. By integrating the maps of chromatin accessibility with bovine embryo transcriptomes and DNA methylomes, we found promoter accessibility and DNA methylation in bovine embryos correlated with both gene activities and CpG densities. Most importantly, we constructed the regulatory networks of stage-specific expressed genes and stage-specific activated genes with three omics datasets in bovine early embryos and revealed conserved and distinctive transcriptional regulatory networks between invivo- and invitro-derived embryos. This comprehensive analysis revealed critical features of the chromatin landscape and epigenetic reprogramming during bovine early embryo development.

A genomic predictor of lifespan in vertebrates

Biological ageing and its mechanistic underpinnings are of immense biomedical and ecological significance. Ageing involves the decline of diverse biological functions and places a limit on a species’ maximum lifespan. Ageing is associated with epigenetic changes involving DNA methylation. Furthermore, an analysis of mammals showed that the density of CpG sites in gene promoters, which are targets for DNA methylation, is correlated with lifespan. Using 252 whole genomes and databases of animal age and promotor sequences, we show a pattern across vertebrates. We also derive a predictive lifespan clock based on CpG density in a selected set of promoters. The lifespan clock accurately predicts maximum lifespan in vertebrates (R2 = 0.76) from the density of CpG sites within only 42 selected promoters. Our lifespan clock provides a wholly new method for accurately estimating lifespan using genome sequences alone and enables estimation of this challenging parameter for both poorly understood and extinct species.

GENE-06. INTEGRATIVE ANALYSIS OF DNA METHYLATION SUGGESTS DOWN-REGULATION OF ONCOGENIC PATHWAYS AND REDUCED SOMATIC MUTATION RATES IN SURVIVAL OUTLIERS OF GLIOBLASTOMA

Abstract The study of survival outliers of glioblastoma (GBM) can provide important clues on gliomagenesis as well as on the ways to alter clinical course on this almost uniformly lethal cancer type. However, there has been little consensus on genetic and epigenetic signatures of the long-term survival outliers of GBM. In this study, we compared the two groups of survival outliers of glioblastoma with IDH wild-type, consisting of the GBM patients who lived longer than 3 years (n=17) and the patients who lived less than 1 year (n=12) in terms of genome-wide DNA methylation profile. Statistical analyses were performed to identify differentially methylated sites between the two groups. Functional implication of long-term survivors of GBM (LTS-GBM)-specific sites were investigated by comprehensive enrichment analyses with genomic and epigenomic features. We found that the genome of LTS-GBM is differentially methylated relative to short-term survivor patients depending on CpG density: hypermethylation near CpG islands (CGIs) and hypomethylation far from CGIs. Interestingly, these two patterns are associated with distinct oncogenic aspects in gliomagenesis. In the LTS-GBM-specific sites distant from CGI, somatic mutations of GBM are enriched with higher DNA methylation, suggesting that the hypomethylation in LTS-GBM can contribute to reduce the rate of somatic mutation. On the other hand, the hypermethylation near CGIs associates with transcriptional downregulation of genes involved in cancer progression pathways. Using independent cohorts of IDH1/2- wild type GBM, we also showed that these two patterns of DNA methylation can be used as molecular markers of LTS-GBM. Our results provide extended understanding of DNA methylation, especially of DNA hypomethylation, in cancer genome and reveal clinical importance of DNA methylation pattern as prognostic markers of GBM.

Radiation-Induced Chromosomal Breaks may be DNA Repair Fragile Sites with Larger-scale Correlations to Eight Double-Strand-Break Related Data Sets over the Human Genome.

In this work, we compared the genomic distribution of common radiation-induced chromosomal breaks to eight different data sets covering the whole human genome. Sites with a high probability of chromatid breakage after exposure to low and high ionization density radiations were often located inside common and rare fragile sites, indicating that they may be a new and more local type of DNA repair-related fragility. Breaks in specific chromosome bands after acute exposure to oil and benzene also showed strong correlation with these sites and fragile sites. In addition, close correlation was found with cytologically detected chiasma and MLH1 immunofluorescence sites and with the HapMap recombination density distributions. Also, of interest, copy number changes occurred predominantly at radiation-induced breaks and fragile sites, at least for breast cancers with poor prognosis, and they decreased weakly but significantly in regions with increasing recombination and CpG density. An increased CpG density is linked to regions of high gene density to secure high-fidelity reproduction and survival. To minimize cancer induction, cancer-related genes are often located in regions of decreased recombination density and/or higher-than-average CpG density. It is compelling that all these data sets were influenced by the cells' handling of double-strand breaks and, more generally, DNA damage on its genome. In fact, the DNA repair genes systematically avoid regions with a high recombination density, as they need to be intact to accurately handle repairable DNA lesions.

Abstract 5206: Role of super-enhancers in HPV+ head and neck squamous cell carcinoma

Abstract Head and neck squamous cell carcinomas (HNSCC) are the sixth leading cause of cancer worldwide with different incidences, mortalities, and prognosis for different subsites. Infection by Human Papilloma Viral (HPV) can cause the development of HPV+ HNSCC, the most rapidly growing population of cancer patients. The molecular biology of HPV+ HNSCC is related to abnormal transcriptional regulation. Super-enhancers (SEs) were recently identified as critical regulators of gene expression during cell differentiation and disease development. SEs are long clusters of enhancers that recruit master transcription factors (TFs) and coactivators to regulate the expression of target genes. We have recently developed a new whole-genome analytical pipeline to optimize ChIP-Seq procedures on primary surgical samples, and this protocol helped us to define the whole-genome distribution of H3K27ac, the hallmark of SEs, in HPV+ HNSCC samples. To detect SEs, we used a novel LILY algorithm that corrects the ChIP-Seq signal for copy number variation and CpG density. The analysis revealed that HPV+ HNSCC-specific SEs regulate genes critical for head and neck cancer development, such as EGFR, TP63, JAG1, IRF6, SMAD3, MAP4K2, SNAI1, TNFAIP3, ANO2, and TNFRSF1A. These genes are located in the vicinity of HNSCC-specific SEs, and the expression of those genes was altered by JQ1, the inhibitor of BRD4, the main component of SE machinery, supporting the role of SE in their regulation. The following Cistrome analysis allowed us to elucidate the TF composition of the SE machinery. Thus, we have found out that P63, P53, SOX2, FOSL1, SMAD3, and SNAI2 are TFs that are the most overrepresented in the HPV+ HNSCC-specific SEs. These comprehensive analyses have revealed novel insight into the HPV+ HNSCC biology and paved the basis for the implications for epigenetic therapeutics for this tumor types, and potentially other virus-related malignancies. Citation Format: Fernando T. Zamuner, Ilya Vorontsov, Emily Flam, Vera Mukhina, Ludmila Danilova, Tingting Ou, Elena Stavrovskaya, Theresa Guo, Eric Windsor, Dylan Z. Kelley, Michael Parfenov, Michael Considine, Elana J. Fertig, Alexander Favorov, Daria A. Gaykalova. Role of super-enhancers in HPV+ head and neck squamous cell carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5206.

Capture of a functionally active methyl-CpG binding domain by an arthropod retrotransposon family.

The repressive capacity of cytosine DNA methylation is mediated by recruitment of silencing complexes by methyl-CpG binding domain (MBD) proteins. Despite MBD proteins being associated with silencing, we discovered that a family of arthropod Copia retrotransposons have incorporated a host-derived MBD. We functionally show how retrotransposon-encoded MBDs preferentially bind to CpG-dense methylated regions, which correspond to transposable element regions of the host genome, in the myriapod Strigamia maritima. Consistently, young MBD-encoding Copia retrotransposons (CopiaMBD) accumulate in regions with higher CpG densities than other LTR-retrotransposons also present in the genome. This would suggest that retrotransposons use MBDs to integrate into heterochromatic regions in Strigamia, avoiding potentially harmful insertions into host genes. In contrast, CopiaMBD insertions in the spider Stegodyphus dumicola genome disproportionately accumulate in methylated gene bodies compared with other spider LTR-retrotransposons. Given that transposons are not actively targeted by DNA methylation in the spider genome, this distribution bias would also support a role for MBDs in the integration process. Together, these data show that retrotransposons can co-opt host-derived epigenome readers, potentially harnessing the host epigenome landscape to advantageously tune the retrotransposition process.

Intricate Regulation of Phosphoenolpyruvate Carboxykinase (PEPCK) Isoforms in Normal Physiology and Disease.

The phosphoenolpyruvate carboxykinase (PEPCK) isoforms are considered as rate-limiting enzymes for gluconeogenesis and glyceroneogenesis pathways. PEPCK exhibits several interesting features such as a) organelle-specific isoforms (cytosolic and a mitochondrial) in vertebrate clade, b) tissue-specific expression of isoforms and c) organism-specific requirement of ATP or GTP as a cofactor. In higher organisms, PEPCK isoforms are intricately regulated and activated through several physiological and pathological stimuli such as corticoids, hormones, nutrient starvation and hypoxia. Isoform-specific transcriptional/translational regulation and their interplay in maintaining glucose homeostasis remain to be fully understood. Mounting evidence indicates the significant involvement of PEPCK isoforms in physiological processes (development and longevity) and in the progression of a variety of diseases (metabolic disorders, cancer, Smith-Magenis syndrome). The present systematic review aimed to assimilate existing knowledge of transcriptional and translational regulation of PEPCK isoforms derived from cell, animal and clinical models. Based on current knowledge and extensive bioinformatics analysis, in this review we have provided a comparative (epi)genetic understanding of PCK1 and PCK2 genes encompassing regulatory elements, disease-associated polymorphisms, copy number variations, regulatory miRNAs and CpG densities. We have also discussed various exogenous and endogenous modulators of PEPCK isoforms and their signaling mechanisms. A comprehensive review of existing knowledge of PEPCK regulation and function may enable identification of the underlying gaps to design new pharmacological strategies and interventions for the diseases associated with gluconeogenesis.

Integrative analysis of DNA methylation suggests down-regulation of oncogenic pathways and reduced somatic mutation rates in survival outliers of glioblastoma

The study of survival outliers of glioblastoma can provide important clues on gliomagenesis as well as on the ways to alter clinical course of this almost uniformly lethal cancer type. However, there has been little consensus on genetic and epigenetic signatures of the long-term survival outliers of glioblastoma. Although the two classical molecular markers of glioblastoma including isocitrate dehydrogenase 1 or 2 (IDH1/2) mutation and O6-methylguanine DNA methyltransferase (MGMT) promoter methylation are associated with overall survival rate of glioblastoma patients, they are not specific to the survival outliers. In this study, we compared the two groups of survival outliers of glioblastoma with IDH wild-type, consisting of the glioblastoma patients who lived longer than 3 years (n = 17) and the patients who lived less than 1 year (n = 12) in terms of genome-wide DNA methylation profile. Statistical analyses were performed to identify differentially methylated sites between the two groups. Functional implication of DNA methylation patterns specific to long-term survivors of glioblastoma were investigated by comprehensive enrichment analyses with genomic and epigenomic features. We found that the genome of long-term survivors of glioblastoma is differentially methylated relative to short-term survivor patients depending on CpG density: hypermethylation near CpG islands (CGIs) and hypomethylation far from CGIs. Interestingly, these two patterns are associated with distinct oncogenic aspects in gliomagenesis. In the long-term survival glioblastoma-specific sites distant from CGI, somatic mutations of glioblastoma are enriched with higher DNA methylation, suggesting that the hypomethylation in long-term survival glioblastoma can contribute to reduce the rate of somatic mutation. On the other hand, the hypermethylation near CGIs associates with transcriptional downregulation of genes involved in cancer progression pathways. Using independent cohorts of IDH1/2- wild type glioblastoma, we also showed that these two patterns of DNA methylation can be used as molecular markers of long-term survival glioblastoma. Our results provide extended understanding of DNA methylation, especially of DNA hypomethylation, in cancer genome and reveal clinical importance of DNA methylation pattern as prognostic markers of glioblastoma.

SAT-551 Thyroid Hormone Effects on the Methylome During Post-Embryonic Development

Thyroid hormones (TH) have been shown to govern post-embryonic development in vertebrates inducing multiples effects. Recently, TH were shown to regulate the expression of genes coding for DNA methylation modifying enzymes. Acting on DNA methylation is an important means of regulating the expression of gene clusters and a central mechanism for coordinating developmental transitions. However, DNA methylation is dynamic and also integrate signals from the constantly changing environment. Certain variation of DNA methylation at these developmental transitions can confer on the genome a risk situation that can have lifelong consequences (cardiovascular and metabolism diseases, cancer or nervous system associated illness). Our objective is to locate TH induced variation of DNA methylation profiles to characterize potential determinants for certain pathologies and / or adaptation to a variation of the environment. We have thus initiated the mapping of DNA methylation following treatment with T3, using a method for capturing methylated DNA (MethylCap) combined with high-throughput sequencing (MethylCap-Seq). Magnetic beads allow for capture and fractionation of methylated DNA by CpG density. Here, we will focus on the highly methylated fraction of DNA. Our model of developmental transition is the amphibian metamorphosis, a post-embryonic process regulated by TH. Additionally, stress during metamorphosis might also induce potential phenotype propagation to adult forms, highlighting its importance in driving adaptation to fluctuating environments. We have compared TH induced methylation profiles in the hind limb and caudal epidermal, two organs with contrasted fate during metamorphosis (respectively morphogenesis and cell death). The results shown that TH affect the methylome with a strong tendency for DNA methylation marks to be removed. The majority of differential methylated regions (DMR) are tissue specific. However, few DMRs are found in both tissues and show either identical or opposite variations. The DMRs are generally not localize near genes regulated or not by the hormones. In conclusion, at post-embryonic transition in amphibians, TH induce mainly tissues-specific DNA demethylation not directly related to gene transcription.

CG dinucleotides enhance promoter activity independent of DNA methylation.

Most mammalian RNA polymerase II initiation events occur at CpG islands, which are rich in CpGs and devoid of DNA methylation. Despite their relevance for gene regulation, it is unknown to what extent the CpG dinucleotide itself actually contributes to promoter activity. To address this question, we determined the transcriptional activity of a large number of chromosomally integrated promoter constructs and monitored binding of transcription factors assumed to play a role in CpG island activity. This revealed that CpG density significantly improves motif-based prediction of transcription factor binding. Our experiments also show that high CpG density alone is insufficient for transcriptional activity, yet results in increased transcriptional output when combined with particular transcription factor motifs. However, this CpG contribution to promoter activity is independent of DNA methyltransferase activity. Together, this refines our understanding of mammalian promoter regulation as it shows that high CpG density within CpG islands directly contributes to an environment permissive for full transcriptional activity.

MeDEStrand: an improved method to infer genome-wide absolute methylation levels from DNA enrichment data

BackgroundDNA methylation of CpG dinucleotides is an essential epigenetic modification that plays a key role in transcription. Widely used DNA enrichment-based methods offer high coverage for measuring methylated CpG dinucleotides, with the lowest cost per CpG covered genome-wide. However, these methods measure the DNA enrichment of methyl-CpG binding, and thus do not provide information on absolute methylation levels. Further, the enrichment is influenced by various confounding factors in addition to methylation status, for example, CpG density. Computational models that can accurately derive absolute methylation levels from DNA enrichment data are needed.ResultsWe developed “MeDEStrand,” a method that uses a sigmoid function to estimate and correct the CpG bias from enrichment results to infer absolute DNA methylation levels. Unlike previous methods, which estimate CpG bias based on reads mapped at the same genomic loci, MeDEStrand processes the reads for the positive and negative DNA strands separately. We compared the performance of MeDEStrand to that of three other state-of-the-art methods “MEDIPS,” “BayMeth,” and “QSEA” on four independent datasets generated using immortalized cell lines (GM12878 and K562) and human primary cells (foreskin fibroblasts and mammary epithelial cells). Based on the comparison of the inferred absolute methylation levels from MeDIP-seq data and the corresponding reduced-representation bisulfite sequencing data from each method, MeDEStrand showed the best performance at high resolution of 25, 50, and 100 base pairs.ConclusionsThe MeDEStrand tool can be used to infer whole-genome absolute DNA methylation levels at the same cost of enrichment-based methods with adequate accuracy and resolution. R package MeDEStrand and its tutorial is freely available for download at https://github.com/jxu1234/MeDEStrand.git.

DNA methylation subpatterns at distinct regulatory regions in human early embryos.

DNA methylation has been investigated for many years, but recent technologies have allowed for single-cell- and single-base-resolution DNA methylation datasets and more accurate assessment of DNA methylation dynamics at the key genomic regions that regulate gene expression in human early embryonic development. In this study, the region from upstream 20 kb to downstream 20 kb of RefSeq gene was selected and divided into 12 distinct regions (up20, up10, up5, up2, 5'UTR, exon, intron, 3'UTR, down2, down5, down10 and down20). The candidate promoter region (TSS ± 2 kb) was further divided into 20 consecutive subregions, which were termed ‘bins’. The DNA methylation dynamics of these regions were systematically analysed along with their effects on gene expression in human early embryos. The dynamic DNA methylation subpatterns at the distinct genomic regions with a focus on promoter regions were mapped. For the 12 distinct genomic regions, up2 and 5'UTR had the lowest DNA methylation levels, and their methylation dynamics were different with other regions. The region 3'UTR had the highest DNA methylation levels, and the correlation analysis with gene expression proved that it was a feature of transcribed genes. For the 20 bins in promoter region, the CpG densities showed a normal distribution pattern, and the trend of the methylated CpG counts was inverse with the DNA methylation levels, especially for the bin 1 (downstream 200 bp of the TSS). Through the correlation analysis between DNA methylation and gene expression, the current study finally revealed that the region bin −4 to 6 (800 bp upstream to 1200 bp downstream of the TSS) was the best candidate for the promoter region in human early embryos, and bin 1 was the putative key regulator of gene activity. This study provided a global and high-resolution view of DNA methylation subpatterns at the distinct genomic regions in human early embryos.

DNA methylation reprogramming of functional elements during mammalian embryonic development

DNA methylation plays important roles during development. However, the DNA methylation reprogramming of functional elements has not been fully investigated during mammalian embryonic development. Herein, using our modified MethylC-Seq library generation method and published post-bisulphite adapter-tagging (PBAT) method, we generated genome-wide DNA methylomes of human gametes and early embryos at single-base resolution and compared them with mouse methylomes. We showed that the dynamics of DNA methylation in functional elements are conserved between humans and mice during early embryogenesis, except for satellite repeats. We further found that oocyte-specific hypomethylated promoters usually exhibit low CpG densities. Genes with oocyte-specific hypomethylated promoters generally show oocyte-specific hypomethylated genic and intergenic regions, and these hypomethylated regions contribute to the hypomethylation pattern of mammalian oocytes. Furthermore, hypomethylated genic regions with low CG densities correlate with gene silencing in oocytes, whereas hypomethylated genic regions with high CG densities correspond to high gene expression. We further show that methylation reprogramming of enhancers during early embryogenesis is highly associated with the development of almost all human organs. Our data support the hypothesis that DNA methylation plays important roles during mammalian development.

DNMT3A and TET1 cooperate to regulate promoter epigenetic landscapes in mouse embryonic stem cells

BackgroundDNA methylation is a heritable epigenetic mark, enabling stable but reversible gene repression. In mammalian cells, DNA methyltransferases (DNMTs) are responsible for modifying cytosine to 5-methylcytosine (5mC), which can be further oxidized by the TET dioxygenases to ultimately cause DNA demethylation. However, the genome-wide cooperation and functions of these two families of proteins, especially at large under-methylated regions, called canyons, remain largely unknown.ResultsHere we demonstrate that DNMT3A and TET1 function in a complementary and competitive manner in mouse embryonic stem cells to mediate proper epigenetic landscapes and gene expression. The longer isoform of DNMT3A, DNMT3A1, exhibits significant enrichment at distal promoters and canyon edges, but is excluded from proximal promoters and canyons where TET1 shows prominent binding. Deletion of Tet1 increases DNMT3A1 binding capacity at and around genes with wild-type TET1 binding. However, deletion of Dnmt3a has a minor effect on TET1 binding on chromatin, indicating that TET1 may limit DNA methylation partially by protecting its targets from DNMT3A and establishing boundaries for DNA methylation. Local CpG density may determine their complementary binding patterns and therefore that the methylation landscape is encoded in the DNA sequence. Furthermore, DNMT3A and TET1 impact histone modifications which in turn regulate gene expression. In particular, they regulate Polycomb Repressive Complex 2 (PRC2)-mediated H3K27me3 enrichment to constrain gene expression from bivalent promoters.ConclusionsWe conclude that DNMT3A and TET1 regulate the epigenome and gene expression at specific targets via their functional interplay.

Conflicts of CpG density and DNA methylation are proximally and distally involved in gene regulation in human and mouse tissues

ABSTRACTThe relationship between CpG content and DNA methylation has attracted considerable interest in recent years. Direct or indirect methods have been developed to investigate their regulatory functions based on various hypotheses, large cohort studies, and meta-analyses. However, all of these analyses were performed at units of CpG blocks and, thus, the influence of finer genome structure has been neglected. Herein, we present a novel algorithm of base-pair resolution to systematically investigate the relationship between CpG contents and DNA methylation. By introducing the concept of ‘complementary index’ we examined the methylomes of 34 adult and 7 embryonic tissues and successfully fitted the relationship of DNA methylation and CpG density into a nonlinear mathematical model. A further algorithm was developed to locate the regions where CpG density does not match expectations from the model, termed ‘conflict of gap’ (COG) regions. Interestingly, COGs are highly concordant in human and mouse and their distributions display a tissue-specific pattern. Based on COG methylation patterns we correctly classified tissues according to their function or origin. We demonstrate that COGs based on our method can reveal more and deeper information than traditional differential methylation region (DMR) approaches. We also found that when COGs are located near to transcription start site (TSS), these regions can determine which promoters will be utilized for initiating gene transcription. Furthermore, COGs located far from the TSS perform as enhancers in terms of histone modification, sequence conservation, transcription factor binding, and DNase I-hypersensitivity.

Switching roles for DNA and histone methylation depend on evolutionary ages of human endogenous retroviruses.

We provide a comprehensive genomic and epigenomic map of the more than 500,000 endogenous retroviruses (ERVs) and fragments that populate the intergenic regions of the human genome. The repressive epigenetic marks associated with the ERVs, particularly long terminal repeats (LTRs), show a remarkable switch in silencing mechanisms, depending on the evolutionary age of the LTRs. Young LTRs tend to be CpG rich and are mainly suppressed by DNA methylation, whereas intermediate age LTRs are associated predominantly with histone modifications, particularly histone H3 lysine 9 (H3K9) methylation. Young LTRs can be reactivated by treatment with the DNA methylation inhibitor 5-aza-2′-deoxycytidine (5-aza-CdR) alone, but their level of expression is much increased by 5-aza-CdR treatment plus knockdown of one of several H3K9 methyltransferases or of the H3K27 methyltransferase EZH2. The removal of cytosine methylation led to rapid, widespread increases in H3K9me3 in the LTRs. Intermediate age LTRs had lower CpG densities and were not up-regulated by 5-aza-CdR treatment, but they were sensitive to knockdown of H3K9 methyltransferases. Unlike the situation in embryonic stem cells, the polycomb repressive complex (PRC2) has a minor role in LTR suppression by itself and is only a player after removal of cytosine methylation in the analyzed cancer cell line. Up-regulation of LTRs and induction of “viral mimicry” is rapidly becoming of interest for predicting cancer patient response to epigenetic therapies. Understanding the mechanism for LTR suppression is of major importance in order to improve patient treatment strategies.

Genome-Scale Oscillations in DNA Methylation during Exit from Pluripotency.

SummaryPluripotency is accompanied by the erasure of parental epigenetic memory, with naïve pluripotent cells exhibiting global DNA hypomethylation both in vitro and in vivo. Exit from pluripotency and priming for differentiation into somatic lineages is associated with genome-wide de novo DNA methylation. We show that during this phase, co-expression of enzymes required for DNA methylation turnover, DNMT3s and TETs, promotes cell-to-cell variability in this epigenetic mark. Using a combination of single-cell sequencing and quantitative biophysical modeling, we show that this variability is associated with coherent, genome-scale oscillations in DNA methylation with an amplitude dependent on CpG density. Analysis of parallel single-cell transcriptional and epigenetic profiling provides evidence for oscillatory dynamics both in vitro and in vivo. These observations provide insights into the emergence of epigenetic heterogeneity during early embryo development, indicating that dynamic changes in DNA methylation might influence early cell fate decisions.

The evolution of CpG density and lifespan in conserved primate and mammalian promoters.

Gene promoters are evolutionarily conserved across holozoans and enriched in CpG sites, the target for DNA methylation. As animals age, the epigenetic pattern of DNA methylation degrades, with highly methylated CpG sites gradually becoming demethylated while CpG islands increase in methylation. Across vertebrates, aging is a trait that varies among species. We used this variation to determine whether promoter CpG density correlates with species’ maximum lifespan. Human promoter sequences were used to identify conserved regions in 131 mammals and a subset of 28 primate genomes. We identified approximately 1000 gene promoters (5% of the total), that significantly correlated CpG density with lifespan. The correlations were performed via the phylogenetic least squares method to account for trait similarity by common descent using phylogenetic branch lengths. Gene set enrichment analysis revealed no significantly enriched pathways or processes, consistent with the hypothesis that aging is not under positive selection. However, within both mammals and primates, 95% of the promoters showed a positive correlation between increasing CpG density and species lifespan, and two thirds were shared between the primate subset and mammalian datasets. Thus, these genes may require greater buffering capacity against age-related dysregulation of DNA methylation in longer-lived species.

No evidence for association of MTHFR 677C>T and 1298A>C variants with placental DNA methylation

Background5,10-Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme in one-carbon metabolism that ensures the availability of methyl groups for methylation reactions. Two single-nucleotide polymorphisms (SNPs) in the MTHFR gene, 677C>T and 1298A>C, result in a thermolabile enzyme with reduced function. These variants, in both the maternal and/or fetal genes, have been associated with pregnancy complications including miscarriage, neural tube defects (NTDs), and preeclampsia (PE), perhaps due to altered capacity for DNA methylation (DNAm). In this study, we assessed the association between MTHFR 677TT and 1298CC genotypes and risk of NTDs, PE, or normotensive intrauterine growth restriction (nIUGR). Additionally, we assessed whether these high-risk genotypes are associated with altered DNAm in the placenta.ResultsIn 303 placentas screened for this study, we observed no significant association between the occurrence of NTDs (N = 55), PE (early-onset: N = 28, late-onset: N = 20), or nIUGR (N = 21) and placental (fetal) MTHFR 677TT or 1298CC genotypes compared to healthy pregnancies (N = 179), though a trend of increased 677TT genotype in PE/IUGR together was observed (OR 2.53, p = 0.048). DNAm was profiled in 10 high-risk 677 (677TT + 1298AA), 10 high-risk 1298 (677CC + 1298CC), and 10 reference (677CC + 1298AA) genotype placentas. Linear modeling identified no significantly differentially methylated sites between high-risk 677 or 1298 and reference placentas at a false discovery rate < 0.05 and Δβ ≥ 0.05 using the Illumina Infinium HumanMethylation450 BeadChip. Using a differentially methylated region analysis or separating cytosine-guanine dinucleotides (CpGs) by CpG density to reduce multiple comparisons also did not identify differential methylation. Additionally, there was no consistent evidence for altered methylation of repetitive DNA between high-risk and reference placentas.ConclusionsWe conclude that large-scale, genome-wide disruption in DNAm does not occur in placentas with the high-risk MTHFR 677TT or 1298CC genotypes. Furthermore, there was no evidence for an association of the 1298CC genotype and only a tendency to higher 677TT in pregnancy complications of PE/IUGR. This may be due to small sample sizes or folate repletion in our Canadian population attenuating effects of the high-risk MTHFR variants. However, given our results and the conflicting results in the literature, investigations into alternative mechanisms that may explain the link between MTHFR variants and pregnancy complications, or in populations at risk of folate deficiencies, are warranted.

Exploring the Link between Nucleosome Occupancy and DNA Methylation.

Near promoters, both nucleosomes and CpG sites form characteristic spatial patterns. Previously, nucleosome depleted regions were observed upstream of transcription start sites and nucleosome occupancy was reported to correlate both with CpG density and the level of CpG methylation. Several studies imply a causal link where CpG methylation might induce nucleosome formation, whereas others argue the opposite, i.e., that nucleosome occupancy might influence CpG methylation. Correlations are indeed evident between nucleosomes, CpG density and CpG methylation—at least near promoter sites. It is however less established whether there is an immediate causal relation between nucleosome occupancy and the presence of CpG sites—or if nucleosome occupancy could be influenced by other factors. In this work, we test for such causality in human genomes by analyzing the three quantities both near and away from promoter sites. For data from the human genome we compare promoter regions with given CpG densities with genomic regions without promoters but of similar CpG densities. We find the observed correlation between nucleosome occupancy and CpG density, respectively CpG methylation, to be specific to promoter regions. In other regions along the genome nucleosome occupancy is statistically independent of the positioning of CpGs or their methylation levels. Anti-correlation between CpG density and methylation level is however similarly strong in both regions. On promoters, nucleosome occupancy is more strongly affected by the level of gene expression than CpG density or CpG methylation—calling into question any direct causal relation between nucleosome occupancy and CpG organization. Rather, our results suggest that for organisms with cytosine methylation nucleosome occupancy might be primarily linked to gene expression, with no strong impact on methylation.