non-CpG Sites Research Articles

Mapping the epigenomic landscape of post-traumatic stress disorder in human cortical neurons.

The study conducted a comprehensive genome-wide analysis of differential 5mC and 5hmC modifications at both CpG and non-CpG sites in postmortem orbitofrontal neurons from 25 PTSD cases and 13 healthy controls. It was observed that PTSD patients exhibit a greater number of differential 5hmC sites compared to 5mC sites. Specifically, individuals with PTSD tend to show hyper-5mC/5hmC at CpG sites, particularly within CpG islands and promoter regions, and hypo-5mC/5hmC at non-CpG sites, especially within intragenic regions. Functional enrichment analysis indicated distinct yet interconnected roles for 5mC and 5hmC in PTSD. The 5mC marks primarily regulate cell-cell adhesion processes, whereas 5hmC marks are involved in embryonic morphogenesis and cell fate commitment. By integrating published PTSD findings from central and peripheral tissues through multi-omics approaches, several biological mechanisms were prioritized, including developmental processes, HPA axis regulation, and immune responses. Based on the consistent enrichment in developmental processes, we hypothesize that if epigenetic changes occur during early developmental stages, they may increase the risk of developing PTSD following trauma exposure. Conversely, if these epigenetic changes occur in adulthood, they may influence neuronal apoptosis and survival mechanisms.

Maternal e-cigarette exposure alters DNA methylome, site-specific CpG and CH methylation, and transcriptomic signatures in the neonatal brain

Maternal use of e-cigarette (e-cig) aerosols poses significant risks to fetal brain development, potentially increasing susceptibility to neurodevelopmental disorders in later life. However, the underlying mechanisms remain incompletely understood. This study aimed to understand the effects of fetal e-cig exposure on DNA methylome and transcriptomic changes in the neonatal brain. Pregnant rats were exposed to e-cig aerosols, and neonatal brains (5 males and 5 females/group) from both control and e-cig-exposed groups were used for experimental analysis. Results indicated that prenatal e-cig exposure altered site-specific DNA methylation patterns at both CpG and CH (non-CpG) sites, predominantly in intergenic and intronic regions, with sex dimorphism in methylation and gene expression changes. Gene ontology analysis revealed that e-cig exposure not only affected neuron projection development and axonogenesis but also altered pathways related to neurodegeneration and long-term depression. These findings provide novel insights into the dynamic changes of CpG and CH methylation induced by e-cig exposure, underscoring the susceptibility of the developing brain to maternal e-cig exposure and its potential implications for developmental disorders and neurodegenerative diseases later in life.

Identifying the “stripe” transcription factors and cooperative binding related to DNA methylation

DNA methylation plays a critical role in gene regulation by modulating the DNA binding of transcription factors (TFs). This study integrates TFs’ ChIP-seq profiles with WGBS profiles to investigate how DNA methylation affects protein interactions. Statistical methods and a 5-letter DNA motif calling model have been developed to characterize DNA sequences bound by proteins, while considering the effects of DNA modifications. By employing these methods, 79 significant universal “stripe” TFs and cofactors (USFs), 2360 co-binding protein pairs, and distinct protein modules associated with various DNA methylation states have been identified. The USFs hint a regulatory hierarchy within these protein interactions. Proteins preferentially bind to non-CpG sites in methylated regions, indicating binding affinity is not solely CpG-dependent. Proteins involved in methylation-specific USFs and cobinding pairs play essential roles in promoting and sustaining DNA methylation through interacting with DNMTs or inhibiting TET binding. These findings underscore the interplay between protein binding and methylation, offering insights into epigenetic regulation in cellular biology.

Structural basis of water-mediated cis Watson-Crick/Hoogsteen base-pair formation in non-CpG methylation.

Non-CpG methylation is associated with several cellular processes, especially neuronal development and cancer, while its effect on DNA structure remains unclear. We have determined the crystal structures of DNA duplexes containing -CGCCG- regions as CCG repeat motifs that comprise a non-CpG site with or without cytosine methylation. Crystal structure analyses have revealed that the mC:G base-pair can simultaneously form two alternative conformations arising from non-CpG methylation, including a unique water-mediated cis Watson-Crick/Hoogsteen, (w)cWH, and Watson-Crick (WC) geometries, with partial occupancies of 0.1 and 0.9, respectively. NMR studies showed that an alternative conformation of methylated mC:G base-pair at non-CpG step exhibits characteristics of cWH with a syn-guanosine conformation in solution. DNA duplexes complexed with the DNA binding drug echinomycin result in increased occupancy of the (w)cWH geometry in the methylated base-pair (from 0.1 to 0.3). Our structural results demonstrated that cytosine methylation at a non-CpG step leads to an anti→syntransition of its complementary guanosine residue toward the (w)cWH geometry as a partial population of WC, in both drug-bound and naked mC:G base pairs. This particular geometry is specific to non-CpG methylated dinucleotide sites in B-form DNA. Overall, the current study provides new insights into DNA conformation during epigenetic regulation.

Unravelling reference bias in ancient DNA datasets.

The alignment of sequencing reads is a critical step in the characterization of ancient genomes. However, reference bias and spurious mappings pose a significant challenge, particularly as cutting-edge wet lab methods generate datasets that push the boundaries of alignment tools. Reference bias occurs when reference alleles are favoured over alternative alleles during mapping, whereas spurious mappings stem from either contamination or when endogenous reads fail to align to their correct position. Previous work has shown that these phenomena are correlated with read length but a more thorough investigation of reference bias and spurious mappings for ancient DNA has been lacking. Here, we use a range of empirical and simulated palaeogenomic datasets to investigate the impacts of mapping tools, quality thresholds, and reference genome on mismatch rates across read lengths. For these analyses, we introduce AMBER, a new bioinformatics tool for assessing the quality of ancient DNA mapping directly from BAM-files and informing on reference bias, read length cut-offs and reference selection. AMBER rapidly and simultaneously computes the sequence read mapping bias in the form of the mismatch rates per read length, cytosine deamination profiles at both CpG and non-CpG sites, fragment length distributions, and genomic breadth and depth of coverage. Using AMBER, we find that mapping algorithms and quality threshold choices dictate reference bias and rates of spurious alignment at different read lengths in a predictable manner, suggesting that optimized mapping parameters for each read length will be a key step in alleviating reference bias and spurious mappings. AMBER is available for noncommercial use on GitHub (https://github.com/tvandervalk/AMBER.git). Scripts used to generate and analyse simulated datasets are available on Github (https://github.com/sdolenz/refbias_scripts).

Comprehensive DNA methylation profiling by MeDIP-NGS identifies potential genes and pathways for epithelial ovarian cancer.

Ovarian cancer, among all gynecologic malignancies, exhibits the highest incidence and mortality rate, primarily because it is often presents with non-specific or no symptoms during its early stages. For the advancement of Ovarian Cancer Diagnosis, it is crucial to identify the potential molecular signatures that could significantly differentiate between healthy and ovarian cancerous tissues and can be used further as a diagnostic biomarker for detecting ovarian cancer. In this study, we investigated the genome-wide methylation patterns in ovarian cancer patients using Methylated DNA Immunoprecipitation (MeDIP-Seq) followed by NGS. Identified differentially methylated regions (DMRs) were further validated by targeted bisulfite sequencing for CpG site-specific methylation profiles. Furthermore, expression validation of six genes by Quantitative Reverse Transcriptase-PCR was also performed. Out of total 120 differentially methylated genes (DMGs), 68 genes were hypermethylated, and 52 were hypomethylated in their promoter region. After analysis, we identified the top 6 hub genes, namely POLR3B, PLXND1, GIGYF2, STK4, BMP2 and CRKL. Interestingly we observed Non-CpG site methylation in the case of POLR3B and CRKL which was statistically significant in discriminating ovarian cancer samples from normal controls. The most significant pathways identified were focal adhesion, the MAPK signaling pathway, and the Ras signaling pathway. Expression analysis of hypermethylated genes was correlated with the downregulation of the genes. POLR3B and GIGYF2 turned out to be the novel genes associated with the carcinogenesis of EOC. Our study demonstrated that methylation profiling through MeDIP-sequencing has effectively identified six potential hub genes and pathways that might exacerbate our understanding of underlying molecular mechanisms of ovarian carcinogenesis.

TET3 is a positive regulator of mitochondrial respiration in Neuro2A cells.

Ten-Eleven-Translocase (TET) enzymes contribute to the regulation of the methylome via successive oxidation of 5-methyl cytosine (5mC) to derivatives which can be actively removed by base-excision-repair (BER) mechanisms in the absence of cell division. This is particularly important in post-mitotic neurons where changes in DNA methylation are known to associate with changes in neural function. TET3, specifically, is a critical regulator of both neuronal differentiation in development and mediates dynamic changes in the methylome of adult neurons associated with cognitive function. While DNA methylation is understood to regulate transcription, little is known of the specific targets of TET3-dependent catalytic activity in neurons. We report the results of an unbiased transcriptome analysis of the neuroblastoma-derived cell line; Neuro2A, in which Tet3 was silenced. Oxidative phosphorylation (OxPhos) was identified as the most significantly down-regulated functional canonical pathway, and these findings were confirmed by measurements of oxygen consumption rate in the Seahorse bioenergetics analyser. The mRNA levels of both nuclear- and mitochondrial-encoded OxPhos genes were reduced by Tet3-silencing, but we found no evidence for differential (hydroxy)methylation deposition at these gene loci. However, the mRNA expression of genes known to be involved in mitochondrial quality control were also shown to be significantly downregulated in the absence of TET3. One of these genes; EndoG, was identified as a direct target of TET3-catalytic activity at non-CpG methylated sites within its gene body. Accordingly, we propose that aberrant mitochondrial homeostasis may contribute to the decrease in OxPhos, observed upon Tet3-downregulation in Neuro2A cells.

Biochemical and structural characterization of the first-discovered metazoan DNA cytosine-N4 methyltransferase from the bdelloid rotifer Adineta vaga

Much is known about the generation, removal, and roles of 5-methylcytosine (5mC) in eukaryote DNA, and there is a growing body of evidence regarding N6-methyladenine, but very little is known about N4-methylcytosine (4mC) in the DNA of eukaryotes. The gene for the first metazoan DNA methyltransferase generating 4mC (N4CMT) was reported and characterized recently by others, in tiny freshwater invertebrates called bdelloid rotifers. Bdelloid rotifers are ancient, apparently asexual animals, and lack canonical 5mC DNA methyltransferases. Here, we characterize the kinetic properties and structural features of the catalytic domain of the N4CMT protein from the bdelloid rotifer Adineta vaga. We find that N4CMT generates high-level methylation at preferred sites, (a/c)CG(t/c/a), and low-level methylation at disfavored sites, exemplified by ACGG. Like the mammalian de novo 5mC DNA methyltransferase 3A/3B (DNMT3A/3B), N4CMT methylates CpG dinucleotides on both DNA strands, generating hemimethylated intermediates and eventually fully methylated CpG sites, particularly in the context of favored symmetric sites. In addition, like DNMT3A/3B, N4CMT methylates non-CpG sites, mainly CpA/TpG, though at a lower rate. Both N4CMT and DNMT3A/3B even prefer similar CpG-flanking sequences. Structurally, the catalytic domain of N4CMT closely resembles the Caulobacter crescentus cell cycle-regulated DNA methyltransferase. The symmetric methylation of CpG, and similarity to a cell cycle-regulated DNA methyltransferase, together suggest that N4CMT might also carry out DNA synthesis-dependent methylation following DNA replication.

Synonymous alterations of cancer-associated Trp53 CpG mutational hotspots cause fatal developmental jaw malocclusions but no tumors in knock-in mice.

Intragenic CpG dinucleotides are tightly conserved in evolution yet are also vulnerable to methylation-dependent mutation, raising the question as to why these functionally critical sites have not been deselected by more stable coding sequences. We previously showed in cell lines that altered exonic CpG methylation can modify promoter start sites, and hence protein isoform expression, for the human TP53 tumor suppressor gene. Here we extend this work to the in vivo setting by testing whether synonymous germline modifications of exonic CpG sites affect murine development, fertility, longevity, or cancer incidence. We substituted the DNA-binding exons 5-8 of Trp53, the mouse ortholog of human TP53, with variant-CpG (either CpG-depleted or -enriched) sequences predicted to encode the normal p53 amino acid sequence; a control construct was also created in which all non-CpG sites were synonymously substituted. Homozygous Trp53-null mice were the only genotype to develop tumors. Mice with variant-CpG Trp53 sequences remained tumor-free, but were uniquely prone to dental anomalies causing jaw malocclusion (p < .0001). Since the latter phenotype also characterises murine Rett syndrome due to dysfunction of the trans-repressive MeCP2 methyl-CpG-binding protein, we hypothesise that CpG sites may exert non-coding phenotypic effects via pre-translational cis-interactions of 5-methylcytosine with methyl-binding proteins which regulate mRNA transcript initiation, expression or splicing, although direct effects on mRNA structure or translation are also possible.

Abstract 2609: Non-CpG sites preference in G:C&amp;gt;A:T transition of TP53 in gastric cancer of Eastern Europe (Poland, Romania and Hungary) compared to East Asian countries (China and Japan) : gastric cancer in the world, mutation spectrum revisited

Abstract Aim: Mutation spectrum of TP53 in gastric cancer (GC) has been investigated world-widely, but a comparison of mutation spectrum among GCs from various regions in the world are still sparsely documented. In order to identify the difference of TP53 mutation spectrum in GCs in Eastern Europe and in East Asia, we sequenced TP53 in GCs from Eastern Europe, Lujiang (China), and Yokohama, Kanagawa (Japan) and identified the feature of TP53 mutations of GC in these regions. Subjects and Method: In total, 689 tissue samples of GC were analyzed: 288 samples from East European populations (25 from Hungary, 71 from Poland and 192 from Romania), 268 from Yokohama, Kanagawa, Japan and 133 from Lujiang, Anhui province, China. DNA was extracted from FFPE tissue of Chinese, East European cases; and from frozen tissue of Japanese GCs. PCR products were direct-sequenced by Sanger method, and in ambiguous cases, PCR product was cloned and up to 8 clones were sequenced. We used No. NC_000017.11(hg38) as the reference sequence of TP53. Mutation patterns were categorized into nine groups: six base substitutions, insertion, deletion and deletion-insertion. Within G:C&amp;gt;A:T mutations the mutations in CpG and non-CpG sites were divided. Furthermore, 40 cases of Rwanda GC were compared in the world for mutation spectrum. The Cancer Genome Atlas data (TCGA, ver.R20, July, 2019) having somatic mutation list of GCs from Whites, Asians, and other ethnicities were used as a reference for our data. Results: The most frequent base substitutions were G:C&amp;gt;A:T transition in all the areas investigated. The G:C&amp;gt;A:T transition in non-CpG sites were prominent in East European GCs, compared with Asian ones. Mutation pattern from TCGA data revealed the same trend between GCs from White (TCGA category) vs Asian countries. Chinese and Japanese GCs showed higher ratio of G:C&amp;gt;A:T transition in CpG sites and A:T&amp;gt;G:C mutation was more prevalent in Asian countries. Conclusion: The divergence in mutation spectrum of GC in different areas in the world may reflect various pathogeneses and etiologies of GC, region to region. Diversified mutation spectrum in GC in Eastern Europe may suggest GC in Europe has different carcinogenic pathway of those from Asia. Whether the Rwanda pattern belong to Asian or East European would be interesting. Citation Format: Hiroko Natsume, Kinga Szczepaniak, Hidetaka Yamada, Yuji Iwashita, Marta Gędek, Jelena Šuto, Rika Kasajima, Felix Manirakiza, Augustin Nzitakera, Rei Ishikawa, Eric Rutaganda, François Ngabonziza, Benoit Seminega, Placide Kamali, Vincent Dusabejambo, Belson Rugwizangoga, Nobuhito Nobuhito, Yaping Wang, Yohei Miyagi, Simona Gurzu, Haruhiko Sugimura. Non-CpG sites preference in G:C&amp;gt;A:T transition of TP53 in gastric cancer of Eastern Europe (Poland, Romania and Hungary) compared to East Asian countries (China and Japan) : gastric cancer in the world, mutation spectrum revisited [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2609.

SNCA Gene Methylation in Parkinson's Disease and Multiple System Atrophy.

In recent years, epigenetic mechanisms have been implicated in the development of multifactorial diseases including neurodegenerative disorders. In Parkinson's disease (PD), as a synucleinopathy, most studies focused on DNA methylation of SNCA gene coding alpha-synuclein but obtained results were rather contradictory. In another neurodegenerative synucleinopathy, multiple system atrophy (MSA), very few studies investigated the epigenetic regulation. This study included patients with PD (n = 82), patients with MSA (n = 24), and a control group (n = 50). In three groups, methylation levels of CpG and non-CpG sites in regulatory regions of the SNCA gene were analyzed. We revealed hypomethylation of CpG sites in the SNCA intron 1 in PD and hypermethylation of predominantly non-CpG sites in the SNCA promoter region in MSA. In PD patients, hypomethylation in the intron 1 was associated with earlier age at the disease onset. In MSA patients, hypermethylation in the promotor was associated with shorter disease duration (before examination). These results showed different patterns of the epigenetic regulation in two synucleinopathies-PD and MSA.

CpH methylome analysis in human cortical neurons identifies novel gene pathways and drug targets for opioid use disorder.

DNA methylation (DNAm), an epigenetic mechanism, has been associated with opioid use disorder (OUD) in preclinical and human studies. However, most of the studies have focused on DNAm at CpG sites. DNAm at non-CpG sites (mCpHs, where H indicates A, T, or C) has been recently shown to have a role in gene regulation and to be highly abundant in neurons. However, its role in OUD is unknown. This work aims to evaluate mCpHs in the human postmortem orbital frontal cortex (OFC) in the context of OUD. A total of 38 Postmortem OFC samples were obtained from the VA Brain Bank (OUD = 12; Control = 26). mCpHs were assessed using reduced representation oxidative bisulfite sequencing in neuronal nuclei. Differential analysis was performed using the "methylkit" R package. Age, ancestry, postmortem interval, PTSD, and smoking status were included as covariates. Significant mCpHs were set at q-value < 0.05. Gene Ontology (GO) and KEGG enrichment analyses were performed for the annotated genes of all differential mCpH loci using String, ShinyGO, and amiGO software. Further, all annotated genes were analyzed using the Drug gene interaction database (DGIdb). A total of 2,352 differentially methylated genome-wide significant mCpHs were identified in OUD, mapping to 2,081 genes. GO analysis of genes with differential mCpH loci showed enrichment for nervous system development (p-value = 2.32E-19). KEGG enrichment analysis identified axon guidance and glutamatergic synapse (FDR 9E-4-2.1E-2). Drug interaction analysis found 3,420 interactions between the annotated genes and drugs, identifying interactions with 15 opioid-related drugs, including lofexidine and tizanidine, both previously used for the treatment of OUD-related symptoms. Our findings suggest a role of mCpHs for OUD in cortical neurons and reveal important biological pathways and drug targets associated with the disorder.

Non-CpG sites preference in G:C > A:T transition of TP53 in gastric cancer of Eastern Europe (Poland, Romania and Hungary) compared to East Asian countries (China and Japan)

AimMutation spectrum of TP53 in gastric cancer (GC) has been investigated world-widely, but a comparison of mutation spectrum among GCs from various regions in the world are still sparsely documented. In order to identify the difference of TP53 mutation spectrum in GCs in Eastern Europe and in East Asia, we sequenced TP53 in GCs from Eastern Europe, Lujiang (China), and Yokohama, Kanagawa (Japan) and identified the feature of TP53 mutations of GC in these regions.Subjects and methodIn total, 689 tissue samples of GC were analyzed: 288 samples from East European populations (25 from Hungary, 71 from Poland and 192 from Romania), 268 from Yokohama, Kanagawa, Japan and 133 from Lujiang, Anhui province, China. DNA was extracted from FFPE tissue of Chinese, East European cases; and from frozen tissue of Japanese GCs. PCR products were direct-sequenced by Sanger method, and in ambiguous cases, PCR product was cloned and up to 8 clones were sequenced. We used No. NC_000017.11(hg38) as the reference sequence of TP53. Mutation patterns were categorized into nine groups: six base substitutions, insertion, deletion and deletion-insertion. Within G:C > A:T mutations the mutations in CpG and non-CpG sites were divided. The Cancer Genome Atlas data (TCGA, ver.R20, July, 2019) having somatic mutation list of GCs from Whites, Asians, and other ethnicities were used as a reference for our data.ResultsThe most frequent base substitutions were G:C > A:T transition in all the areas investigated. The G:C > A:T transition in non-CpG sites were prominent in East European GCs, compared with Asian ones. Mutation pattern from TCGA data revealed the same trend between GCs from White (TCGA category) vs Asian countries. Chinese and Japanese GCs showed higher ratio of G:C > A:T transition in CpG sites and A:T > G:C mutation was more prevalent in Asian countries.ConclusionThe divergence in mutation spectrum of GC in different areas in the world may reflect various pathogeneses and etiologies of GC, region to region. Diversified mutation spectrum in GC in Eastern Europe may suggest GC in Europe has different carcinogenic pathway of those from Asia.

DNA methyltransferase DNMT3A forms interaction networks with the CpG site and flanking sequence elements for efficient methylation

Specific DNA methylation at CpG and non-CpG sites is essential for chromatin regulation. The DNA methyltransferase DNMT3A interacts with target sites surrounded by variable DNA sequences with its TRD and RD loops, but the functional necessity of these interactions is unclear. We investigated CpG and non-CpG methylation in a randomized sequence context using WT DNMT3A and several DNMT3A variants containing mutations at DNA-interacting residues. Our data revealed that the flanking sequence of target sites between the −2 and up to the +8 position modulates methylation rates >100-fold. Non-CpG methylation flanking preferences were even stronger and favor C(+1). R836 and N838 in concert mediate recognition of the CpG guanine. R836 changes its conformation in a flanking sequence-dependent manner and either contacts the CpG guanine or the +1/+2 flank, thereby coupling the interaction with both sequence elements. R836 suppresses activity at CNT sites but supports methylation of CAC substrates, the preferred target for non-CpG methylation of DNMT3A in cells. N838 helps to balance this effect and prevent the preference for C(+1) from becoming too strong. Surprisingly, we found L883 reduces DNMT3A activity despite being highly conserved in evolution. However, mutations at L883 disrupt the DNMT3A-specific DNA interactions of the RD loop, leading to altered flanking sequence preferences. Similar effects occur after the R882H mutation in cancer cells. Our data reveal that DNMT3A forms flexible and interdependent interaction networks with the CpG guanine and flanking residues that ensure recognition of the CpG and efficient methylation of the cytosine in contexts of variable flanking sequences.

High-resolution analysis for urinary DNA jagged ends

Single-stranded ends of double-stranded DNA (jagged ends) are more abundant in urinary DNA than in plasma DNA. However, the lengths of jagged ends in urinary DNA remained undetermined, as a previous method used for urinary DNA jagged end sequencing analysis (Jag-seq) relied on unmethylation at CpG sites, limiting the resolution. Here, we performed high-resolution Jag-seq analysis using methylation at non-CpG cytosine sites, allowing determination of exact length of jagged ends. The urinary DNA bore longer jagged ends (~26-nt) than plasma DNA (~17-nt). The jagged end length distribution displayed 10-nt periodicities in urinary DNA, which were much less observable in plasma DNA. Amplitude of the 10-nt periodicities increased in patients with renal cell carcinoma. Heparin treatment of urine diminished the 10-nt periodicities. The urinary DNA jagged ends often extended into nucleosomal cores, suggesting potential interactions with histones. This study has thus advanced our knowledge of jagged ends in urine DNA.

Flanking sequences influence the activity of TET1 and TET2 methylcytosine dioxygenases and affect genomic 5hmC patterns

TET dioxygenases convert 5-methylcytosine (5mC) preferentially in a CpG context into 5-hydroxymethylcytosine (5hmC) and higher oxidized forms, thereby initiating DNA demethylation, but details regarding the effects of the DNA sequences flanking the target 5mC site on TET activity are unknown. We investigated oxidation of libraries of DNA substrates containing one 5mC or 5hmC residue in randomized sequence context using single molecule readout of oxidation activity and sequence and show pronounced 20 and 70-fold flanking sequence effects on the catalytic activities of TET1 and TET2, respectively. Flanking sequence preferences were similar for TET1 and TET2 and also for 5mC and 5hmC substrates. Enhanced flanking sequence preferences were observed at non-CpG sites together with profound effects of flanking sequences on the specificity of TET2. TET flanking sequence preferences are reflected in genome-wide and local patterns of 5hmC and DNA demethylation in human and mouse cells indicating that they influence genomic DNA modification patterns in combination with locus specific targeting of TET enzymes.

Comparison of EM-seq and PBAT methylome library methods for low-input DNA

ABSTRACT DNA methylation is the most studied epigenetic mark involved in regulation of gene expression. For low input samples, a limited number of methods for quantifying DNA methylation genome-wide has been evaluated. Here, we compared a series of input DNA amounts (1–10ng) from two methylome library preparation protocols, enzymatic methyl-seq (EM-seq) and post-bisulfite adaptor tagging (PBAT) adapted from single-cell PBAT. EM-seq takes advantage of enzymatic activity while PBAT relies on conventional bisulfite conversion for detection of DNA methylation. We found that both methods accurately quantified DNA methylation genome-wide. They produced expected distribution patterns around genomic features, high C-T transition efficiency at non-CpG sites and high correlation between input amounts. However, EM-seq performed better in regard to library and sequencing quality, i.e. EM-seq produced larger insert sizes, higher alignment rates and higher library complexity with lower duplication rate compared to PBAT. Moreover, EM-seq demonstrated higher CpG coverage, better CpG site overlap and higher consistency between input series. In summary, our data suggests that EM-seq overall performed better than PBAT in whole-genome methylation quantification of low input samples.

Demethylation of Non-CpG Sites in DNA Is Initiated by TET2 5-Methylcytosine Dioxygenase

In the mammalian genome, cytosine methylation predominantly occurs at CpG sites. In addition, a number of recent studies have uncovered extensive C5 cytosine methylation (5mC) at non-CpG (5mCpH, where H = A/C/T) sites. Little is known about the enzyme responsible for active demethylation of 5mCpH sites. Using a very sensitive and quantitative LC–MS/MS method, we demonstrate that the human TET2, an iron (II)- and 2OG-dependent dioxygenase, which is a frequently mutated gene in several myeloid malignancies, as well as in a number of other types of cancers, can oxidize 5mCpH sites in double-stranded DNA in vitro. Similar to oxidation of 5mCpG, oxidation of 5mC at CpH sites produces 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxycytosine (5caC) bases in DNA. After 5mCpG, which is the most preferred substrate, TET2 prefers 5mCpC as a substrate, followed by 5mCpA and then 5mCpT. Since the TDG/BER pathway and deformylation or decarboxylation of 5fC or 5caC, respectively, can convert 5fCpH and 5caCpH to an unmodified cytosine base in DNA, our results suggest a novel demethylation pathway of 5mCpH sites initiated by TET2 dioxygenase.

DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation

DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, reduced representation bisulfite sequencing (RRBS) profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, this study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

Human non-CpG methylation patterns display both tissue-specific and inter-individual differences suggestive of underlying function

ABSTRACT DNA methylation (DNAm) in mammals is mostly examined within the context of CpG dinucleotides. Non-CpG DNAm is also widespread across the human genome, but the functional relevance, tissue-specific disposition, and inter-individual variability has not been widely studied. Our aim was to examine non-CpG DNAm in the wider methylome across multiple tissues from the same individuals to better understand non-CpG DNAm distribution within different tissues and individuals and in relation to known genomic regulatory features. DNA methylation in umbilical cord and cord blood at birth, and peripheral venous blood at age 12–13 y from 20 individuals from the Southampton Women’s Survey cohort was assessed by Agilent SureSelect methyl-seq. Hierarchical cluster analysis (HCA) was performed on CpG and non-CpG sites and stratified by specific cytosine environment. Analysis of tissue and inter-individual variation was then conducted in a second dataset of 12 samples: eight muscle tissues, and four aliquots of cord blood pooled from two individuals. HCA using methylated non-CpG sites showed different clustering patterns specific to the three base-pair triplicate (CNN) sequence. Analysis of CAC sites with non-zero methylation showed that samples clustered first by tissue type, then by individual (as observed for CpG methylation), while analysis using non-zero methylation at CAT sites showed samples grouped predominantly by individual. These clustering patterns were validated in an independent dataset using cord blood and muscle tissue. This research suggests that CAC methylation can have tissue-specific patterns, and that individual effects, either genetic or unmeasured environmental factors, can influence CAT methylation.