Search

BackgroundBread wheat has a large complex genome that makes whole genome resequencing costly. Therefore, genome complexity reduction techniques such as sequence capture make re-sequencing cost effective. With a high-quality draft wheat genome now available it is possible to design capture probe sets and to use them to accurately genotype and anchor SNPs to the genome. Furthermore, in addition to genetic variation, epigenetic variation provides a source of natural variation contributing to changes in gene expression and phenotype that can be profiled at the base pair level using sequence capture coupled with bisulphite treatment. Here, we present a new 12 Mbp wheat capture probe set, that allows both the profiling of genotype and methylation from the same DNA sample. Furthermore, we present a method, based on Agilent SureSelect Methyl-Seq, that will use a single capture assay as a starting point to allow both DNA sequencing and methyl-seq.ResultsOur method uses a single capture assay that is sequentially split and used for both DNA sequencing and methyl-seq. The resultant genotype and epi-type data is highly comparable in terms of coverage and SNP/methylation site identification to that generated from separate captures for DNA sequencing and methyl-seq. Furthermore, by defining SNP frequencies in a diverse landrace from the Watkins collection we highlight the importance of having genotype data to prevent false positive methylation calls. Finally, we present the design of a new 12 Mbp wheat capture and demonstrate its successful application to re-sequence wheat.ConclusionWe present a cost-effective method for performing both DNA sequencing and methyl-seq from a single capture reaction thus reducing reagent costs, sample preparation time and DNA requirements for these complementary analyses.

DNA methylation was the first epigenetic modification to be detected in human cancers with specific relation to aberrant gene expression. Herein, DNA methylation analysis explains how epigenetic patterns affect gene expression level. Hypermethylation at tumor suppressor gene loci leads to increased tumorigenesis due to tumor suppressor gene silencing, whereas global hypomethylation of CpG islands (CGIs) is followed by genomic instability and aberrant activation of multiple oncogenes. Therefore, characterization of the genes which silenced or activated epigenetically in human tumor cells can improve our understanding of cancer biology. Different genome-wide methodologies are applied to evaluate methylation status. Various commonly conducted techniques for this evaluation are reviewed in this paper. We provided comparative description of the procedures, advantages, and drawbacks of genome-wide DNA methylation analysis methods and biological applications, to give information on selecting the appropriate method for different methylation studies.

Hypermethylation of CpG islands in the promoter region of many tumor suppressor genes downregulates their expression and in a result promotes tumorigenesis. Therefore, detection of DNA methylation status is a convenient diagnostic tool for cancer detection. Here, we reported a novel method for the integrative detection of methylation by the microfluidic chip-based digital PCR. This method relies on methylation-sensitive restriction enzyme HpaII, which cleaves the unmethylated DNA strands while keeping the methylated ones intact. After HpaII treatment, the DNA methylation level is determined quantitatively by the microfluidic chip-based digital PCR with the lower limit of detection equal to 0.52%. To validate the applicability of this method, promoter methylation of two tumor suppressor genes (PCDHGB6 and HOXA9) was tested in 10 samples of early stage lung adenocarcinoma and their adjacent non-tumorous tissues. The consistency was observed in the analysis of these samples using our method and a conventional bisulfite pyrosequencing. Combining high sensitivity and low cost, the microfluidic chip-based digital PCR method might provide a promising alternative for the detection of DNA methylation and early diagnosis of epigenetics-related diseases.

Epigenetic modifications, in particular DNA methylation, are gaining increasing interest as complementary information to DNA mutations for cancer diagnostics and prognostics. We introduce a method to simultaneously profile DNA mutation and methylation events for an array of sites with single site specificity. Genomic (mutation) or bisulphite-treated (methylation) DNA is amplified using nondiscriminatory primers, and the amplicons are then hybridized to a giant magnetoresistive (GMR) biosensor array followed by melting curve measurements. The GMR biosensor platform offers scalable multiplexed detection of DNA hybridization, which is insensitive to temperature variation. The melting curve approach further enhances the assay specificity and tolerance to variations in probe length. We demonstrate the utility of this method by simultaneously profiling five mutation and four methylation sites in human melanoma cell lines. The method correctly identified all mutation and methylation events and further provided quantitative assessment of methylation density validated by bisulphite pyrosequencing.

Background Formalin fixed paraffin embedded (FFPE) tumor samples are a major source of DNA from patients in cancer research. However, FFPE is a challenging material to work with due to macromolecular fragmentation and nucleic acid crosslinking. FFPE tissue particularly possesses challenges for methylation analysis and for preparing sequencing-based libraries relying on bisulfite conversion. Successful bisulfite conversion is a key requirement for sequencing-based methylation analysis. Methods Here we describe a complete and streamlined workflow for preparing next generation sequencing libraries for methylation analysis from FFPE tissues. This includes, counting cells from FFPE blocks and extracting DNA from FFPE slides, testing bisulfite conversion efficiency with a polymerase chain reaction (PCR) based test, preparing reduced representation bisulfite sequencing libraries and massively parallel sequencing. Results The main features and advantages of this protocol are:An optimized method for extracting good quality DNA from FFPE tissues.An efficient bisulfite conversion and next generation sequencing library preparation protocol that uses 50 ng DNA from FFPE tissue.Incorporation of a PCR-based test to assess bisulfite conversion efficiency prior to sequencing. Conclusions We provide a complete workflow and an integrated protocol for performing DNA methylation analysis at the genome-scale and we believe this will facilitate clinical epigenetic research that involves the use of FFPE tissue. Electronic supplementary material The online version of this article (10.1186/s12920-017-0290-1) contains supplementary material, which is available to authorized users.

DNA methylation is one of the most extensively studied epigenetic modifications of genomic DNA. In recent years, sequencing of bisulfite-converted DNA, particularly via next-generation sequencing technologies, has become a widely popular method to study DNA methylation. This method can be readily applied to a variety of species, dramatically expanding the scope of DNA methylation studies beyond the traditionally studied human and mouse systems. In parallel to the increasing wealth of genomic methylation profiles, many statistical tools have been developed to detect differentially methylated loci (DMLs) or differentially methylated regions (DMRs) between biological conditions. We discuss and summarize several key properties of currently available tools to detect DMLs and DMRs from sequencing of bisulfite-converted DNA. However, the majority of the statistical tools developed for DML/DMR analyses have been validated using only mammalian data sets, and less priority has been placed on the analyses of invertebrate or plant DNA methylation data. We demonstrate that genomic methylation profiles of non-mammalian species are often highly distinct from those of mammalian species using examples of honey bees and humans. We then discuss how such differences in data properties may affect statistical analyses. Based on these differences, we provide three specific recommendations to improve the power and accuracy of DML and DMR analyses of invertebrate data when using currently available statistical tools. These considerations should facilitate systematic and robust analyses of DNA methylation from diverse species, thus advancing our understanding of DNA methylation.

Background There is increasing evidence for inter-individual methylation differences at CpG dinucleotides in the human genome, but the regional extent and function of these differences have not yet been studied in detail. For identifying regions of common methylation differences, we used whole genome bisulfite sequencing data of monocytes from five donors and a novel bioinformatic strategy. Results We identified 157 differentially methylated regions (DMRs) with four or more CpGs, almost none of which has been described before. The DMRs fall into different chromatin states, where methylation is inversely correlated with active, but not repressive histone marks. However, methylation is not correlated with the expression of associated genes. High-resolution single nucleotide polymorphism (SNP) genotyping of the five donors revealed evidence for a role of cis-acting genetic variation in establishing methylation patterns. To validate this finding in a larger cohort, we performed genome-wide association studies (GWAS) using SNP genotypes and 450k array methylation data from blood samples of 1128 individuals. Only 30/157 (19%) DMRs include at least one 450k CpG, which shows that these arrays miss a large proportion of DNA methylation variation. In most cases, the GWAS peak overlapped the CpG position, and these regions are enriched for CREB group, NF-1, Sp100 and CTCF binding motifs. In two cases, there was tentative evidence for a trans-effect by KRAB zinc finger proteins. Conclusions Allele-specific DNA methylation occurs in discrete chromosomal regions and is driven by genetic variation in cis and trans, but in general has little effect on gene expression. Electronic supplementary material The online version of this article (doi:10.1186/s13072-017-0144-2) contains supplementary material, which is available to authorized users.

Mammalian genomes exhibit complex patterns of gene expression regulated, in part, by DNA methylation. The advent of engineered DNA methyltransferases (MTases) to target DNA methylation to specific sites in the genome will accelerate many areas of biological research. However, targeted MTases require clear design rules to direct site-specific DNA methylation and minimize the unintended effects of off-target DNA methylation. Here we report a targeted MTase composed of an artificially split CpG MTase (sMTase) with one fragment fused to a catalytically-inactive Cas9 (dCas9) that directs the functional assembly of sMTase fragments at the targeted CpG site. We precisely map RNA-programmed DNA methylation to targeted CpG sites as a function of distance and orientation from the protospacer adjacent motif (PAM). Expression of the dCas9-sMTase in mammalian cells led to predictable and efficient (up to ~70%) DNA methylation at targeted sites. Multiplexing sgRNAs enabled targeting methylation to multiple sites in a single promoter and to multiple sites in multiple promoters. This programmable de novo MTase tool might be used for studying mechanisms of initiation, spreading and inheritance of DNA methylation, and for therapeutic gene silencing.

Comprehensive studies have shown that DNA methylation plays vital roles in both loss of pluripotency and governance of the transcriptome during embryogenesis and subsequent developmental processes. Aberrant DNA methylation patterns have been widely observed in tumorigenesis, ageing and neurodegenerative diseases, highlighting the importance of a systematic understanding of DNA methylation and the dynamic changes of methylomes during disease onset and progression. Here we describe a facile and convenient approach for efficient targeted DNA methylation by fusing inactive Cas9 (dCas9) with an engineered prokaryotic DNA methyltransferase MQ1. Our study presents a rapid and efficient strategy to achieve locus-specific cytosine modifications in the genome without obvious impact on global methylation in 24 h. Finally, we demonstrate our tool can induce targeted CpG methylation in mice by zygote microinjection, thereby demonstrating its potential utility in early development., Understanding how DNA methylation regulates gene expression requires the capacity to deploy it to regions of interest. The authors generate a highly rapid and locus-specific CpG methylation tool by fusing dCas9 to MQ1 DNA methyltransferase and show efficacy at multiple sites in vitro and in vivo.

Aberrant DNA methylation is associated with various diseases. Quantitative analysis of regional DNA methylation levels of some specific genes would aid in diseases diagnosis and risk stratification. In this study, we developed a robust method for detecting DNA methylation level over targeted genomic regions using nucleobases quantification in bisulfite amplicons by isotope dilution liquid chromatography tandem mass spectrometry coupled with a simple equation. This method had wide detection range (from 0% to 100% methylation) and high accuracy while more time-saving compared to clonal bisulfite sequencing method. The application for clinical tissue samples showed good applicability and cost effectiveness. This analytical method is suitable for quantifying average DNA methylation level over targeted genomic regions and expected to be a useful tool for detecting DNA methylation biomarkers.

Background Cytosine methylation in plant genomes is important for the regulation of gene transcription and transposon activity. Genome-wide methylomes are studied upon mutation of the DNA methyltransferases, adaptation to environmental stresses or during development. However, from basic biology to breeding programs, there is a need to monitor multiple samples to determine transgenerational methylation inheritance or differential cytosine methylation. Methylome data obtained by sodium hydrogen sulfite (bisulfite)-conversion and next-generation sequencing (NGS) provide genome-wide information on cytosine methylation. However, a profiling method that detects cytosine methylation state dispersed over the genome would allow high-throughput analysis of multiple plant samples with distinct epigenetic signatures. We use specific restriction endonucleases to enrich for cytosine coverage in a bisulfite and NGS-based profiling method, which was compared to whole-genome bisulfite sequencing of the same plant material. Methods We established an effective methylome profiling method in plants, termed plant-reduced representation bisulfite sequencing (plant-RRBS), using optimized double restriction endonuclease digestion, fragment end repair, adapter ligation, followed by bisulfite conversion, PCR amplification and NGS. We report a performant laboratory protocol and a straightforward bioinformatics data analysis pipeline for plant-RRBS, applicable for any reference-sequenced plant species. Results As a proof of concept, methylome profiling was performed using an Oryza sativa ssp. indica pure breeding line and a derived epigenetically altered line (epiline). Plant-RRBS detects methylation levels at tens of millions of cytosine positions deduced from bisulfite conversion in multiple samples. To evaluate the method, the coverage of cytosine positions, the intra-line similarity and the differential cytosine methylation levels between the pure breeding line and the epiline were determined. Plant-RRBS reproducibly covers commonly up to one fourth of the cytosine positions in the rice genome when using MspI-DpnII within a group of five biological replicates of a line. The method predominantly detects cytosine methylation in putative promoter regions and not-annotated regions in rice. Conclusions Plant-RRBS offers high-throughput and broad, genome-dispersed methylation detection by effective read number generation obtained from reproducibly covered genome fractions using optimized endonuclease combinations, facilitating comparative analyses of multi-sample studies for cytosine methylation and transgenerational stability in experimental material and plant breeding populations. Electronic supplementary material The online version of this article (doi:10.1186/s12870-017-1070-y) contains supplementary material, which is available to authorized users.

Despite numerous studies done on understanding the role of DNA methylation, limited work has focused on systems integration of cell type-specific interplay between DNA methylation and gene transcription. Through a genome-wide analysis of DNA methylation across 19 cell types with T-47D as reference, we identified 106,252 cell type-specific differentially-methylated CpGs categorized into 7,537 differentially (46.6% hyper- and 53.4% hypo-) methylated regions. We found 44% promoter regions and 75% CpG islands were T-47D cell type-specific methylated. Pyrosequencing experiments validated the cell type-specific methylation across three benchmark cell lines. Interestingly, these DMRs overlapped with 1,145 known tumor suppressor genes. We then developed a Bayesian Gaussian Regression model to measure the relationship among DNA methylation, genomic segment distribution, differential gene expression and tumor suppressor gene status. The model uncovered that 3′UTR methylation has much less impact on transcriptional activity than other regions. Integration of DNA methylation and 82 transcription factor binding information across the 19 cell types suggested diverse interplay patterns between the two regulators. Our integrative analysis reveals cell type-specific and genomic region-dependent regulatory patterns and provides a perspective for integrating hundreds of various omics-seq data together.

Background Differential DNA methylation plays a critical role in the regulation of imprinted genes. The differentially methylated state of the imprinting control region is inherited via the gametes at fertilization, and is stably maintained in somatic cells throughout development, influencing the expression of genes across the imprinting cluster. In contrast, DNA methylation patterns are more labile at secondary differentially methylated regions which are established at imprinted loci during post-implantation development. To investigate the nature of these more variably methylated secondary differentially methylated regions, we adopted a hairpin linker bisulfite mutagenesis approach to examine CpG dyad methylation at differentially methylated regions associated with the murine Dlk1/Gtl2 imprinting cluster on both complementary strands. Results We observed homomethylation at greater than 90% of the methylated CpG dyads at the IG-DMR, which serves as the imprinting control element. In contrast, homomethylation was only observed at 67–78% of the methylated CpG dyads at the secondary differentially methylated regions; the remaining 22–33% of methylated CpG dyads exhibited hemimethylation. Conclusions We propose that this high degree of hemimethylation could explain the variability in DNA methylation patterns at secondary differentially methylated regions associated with imprinted loci. We further suggest that the presence of 5-hydroxymethylation at secondary differentially methylated regions may result in hemimethylation and methylation variability as a result of passive and/or active demethylation mechanisms. Electronic supplementary material The online version of this article (doi:10.1186/s13072-017-0138-0) contains supplementary material, which is available to authorized users.

Aim: Validation of sequencing-based DNA methylation data is an important step for meaningful translation of findings. However, there has been limited assessment of different platforms to validate methylation data from next generation sequencing. Methods: We performed a comparative methylation analysis between the genome-wide platform of reduced representation bisulfite sequencing with a targeted, Sequenom EpiTyper platform (four genes were analyzed in 15 cell lines covering 52 CpG sites). Results: We show that the accuracy of validation substantially improves if results from multiple and adjacent CpG sites are combined rather than at single CpG sites. We demonstrate increased read number improves accuracy of reduced representation bisulfite sequencing results. Further, by using series of replicates, we document variation in samples analyzed by Sequenom EpiTyper, indicating the importance of including replicates to increase precision. Conclusion: The results reveal potential sources of bias and provide a guideline for refining study design for DNA methylation analysis.

Methylated cytosines (5mCs) are frequently mutated in the genome. However, no studies have yet comprehensively analysed mutation–methylation associations across cancer types. Here we analyse 916 cancer genomes, together with tissue type-specific methylation and replication timing data. We describe a strong mutation–methylation association across colorectal cancer subtypes, most interestingly in samples with microsatellite instability (MSI) or Polymerase epsilon (POLE) exonuclease domain mutations. By analysing genomic regions with differential mismatch repair (MMR) efficiency, we suggest a possible role for MMR in the correction of 5mC deamination events, potentially accounting for the high rate of 5mC mutation accumulation in MSI tumours. Additionally, we propose that mutant POLE asserts a mutator phenotype specifically at 5mCs, and we find coding mutation hotspots in POLE-mutant cancers at highly-methylated CpGs in the tumour-suppressor genes APC and TP53. Finally, using multivariable regression models, we demonstrate that different cancers exhibit distinct mutation–methylation associations, with DNA repair influencing such associations in certain cancer genomes. Taken together, we find differential associations with methylation that are vital for accurately predicting expected mutation loads across cancer types. Our findings reveal links between methylation and common mutation and repair processes, with these mechanisms defining a key part of the mutational landscape of cancer genomes.

Background DNA methylation is an important epigenetic mechanism involved in gene regulation, with alterations in DNA methylation in the nuclear genome being linked to numerous complex diseases. Mitochondrial DNA methylation is a phenomenon that is receiving ever-increasing interest, particularly in diseases characterized by mitochondrial dysfunction; however, most studies have been limited to the investigation of specific target regions. Analyses spanning the entire mitochondrial genome have been limited, potentially due to the amount of input DNA required. Further, mitochondrial genetic studies have been previously confounded by nuclear-mitochondrial pseudogenes. Methylated DNA Immunoprecipitation Sequencing is a technique widely used to profile DNA methylation across the nuclear genome; however, reads mapped to mitochondrial DNA are often discarded. Here, we have developed an approach to control for nuclear-mitochondrial pseudogenes within Methylated DNA Immunoprecipitation Sequencing data. We highlight the utility of this approach in identifying differences in mitochondrial DNA methylation across regions of the human brain and pre-mortem blood. Results We were able to correlate mitochondrial DNA methylation patterns between the cortex, cerebellum and blood. We identified 74 nominally significant differentially methylated regions (p

The genome-wide investigation of DNA methylation levels has been limited to reference transposable element positions. The methylation analysis of non-reference and mobile transposable elements has only recently been performed, but required both genome resequencing and MethylC-seq datasets. We have created epiTEome, a program that detects both new transposable element insertion sites and their methylation states from a single MethylC-seq dataset. EpiTEome outperforms other split-read insertion site detection programs, even while functioning on bisulfite-converted reads. EpiTEome characterizes the previously discarded fraction of DNA methylation at sites of new insertions, enabling future investigation into the epigenetic regulation of non-reference and transposed elements.

DNA methylation and demethylation is part of the essential biological processes regulating gene expression in normal cell development. Abnormal methylation status of specific genes and their irregularly translated proteins are normally associated with certain kinds of diseases or cancer. The rapid development of innovative DNA methylation mapping techniques provides a better understanding of DNA methylation pattern and its mechanisms in the human genome and its correlation with numerous diseases. These new techniques can lead us to develop new epigenetic medications, such as DNA methyltransferase inhibitors. As part of the approaches to probe DNA methylation and evaluate the effects of epigenetic therapy, mass spectrometry has been taking an important role in the identification, validation, and quantification of DNA methylation and demethylation. In this review, we will briefly summarize the current breadth of knowledge on the topic of DNA methylation and its occurrence in diseases, DNA methylation drugs, and mass spectrometry based approaches used to study DNA methylation.

We empirically examined the strengths and weaknesses of two human genome-wide DNA methylation platforms: rapid multiplexed reduced representation bisulfite sequencing and Illumina’s Infinium BeadChip. Rapid multiplexed reduced representation bisulfite sequencing required less input DNA, offered more flexibility in coverage, and interrogated more CpG loci at a higher regional density. The Infinium covered slightly more protein coding, cancer-associated and mitochondrial-related genes, both platforms covered all known imprinting clusters, and rapid multiplexed reduced representation bisulfite sequencing covered more microRNA genes than the HumanMethylation450, but fewer than the MethylationEPIC. Rapid multiplexed reduced representation bisulfite sequencing did not always interrogate exactly the same CpG loci, but genomic tiling improved overlap between different libraries. Reproducibility of rapid multiplexed reduced representation bisulfite sequencing and concordance between the platforms increased with CpG density. Only rapid multiplexed reduced representation bisulfite sequencing could genotype samples and measure allele-specific methylation, and we confirmed that Infinium measurements are influenced by nearby single-nucleotide polymorphisms. The respective strengths and weaknesses of these two genome-wide DNA methylation platforms need to be considered when conducting human epigenetic studies., Epigenetics: choose your DNA methylation probing tool wisely! Researchers who study human epigenetics need to carefully consider the platform used to measure genome-wide patterns of DNA methylation. A team led by Karin Michels and Andrea Baccarelli from Harvard University in Boston, Massachusetts, USA, empirically examined the strengths and weaknesses of two methylation profiling tools: Illumina’s Infinium BeadChip, which uses a microarray system to interrogate hundreds of thousands of methylation sites across the genome at single-nucleotide resolution; and a high-throughput sequencing-based approach known as rapid multiplexed reduced representation bisulfite sequencing, or rmRRBS. The former did a better job at reading methylation in protein-coding and mitochondrial-related genes, while the latter required less input DNA and covered more methylation sites across the genome. The authors conclude that a scientist’s platform preference should depend on the nature of his or her investigation.

There is a clear need in molecular markers of prostate cancer (PC) that would specify stratification of patients into risk categories and serve as additional parameters to clinical prognostic factors. Altered DNA methylation status is one of molecular processes that occur during PC development. The gene methylation status may be evaluated in laboratory diagnostics by a PCR approach known as means of methylation-sensitive high resolution melting (MS-HRM) curve analysis. This method is particularly promising due to minimal amounts of DNA required for analysis (10 ng). To date, numerous data based on the use of Infinium HumanMethylation450 BeadChip (HM450) technology have been accumulated on DNA methylation in PC samples. However, it remains unclear whether the MS-HRM results are consistent with chip hybridization data. The aim of this study was to perform the comparative analysis of the diagnostic informativeness of MS-HRM and HM450 chip hybridization methods in determination of the methylation status of three genes GSTP1, APC, and RASSF1 in samples obtained from PC patients. The study has shown that using the MS-HRM method it is possible to discriminate PC tumor samples from samples of histologically unaltered prostate tissue by analyzing the methylation status of each gene. Analysis of chip hybridization data confirmed the results obtained by MS-HRM. Differences between tumor tissue and histologically intact tissue obtained from each PC patient in the methylation levels of GSTP1, APC, and RASSF1, determined by MS-HRM or chip hybridization were consistent with each other. Thus, the assessment of the methylation level of gene GSTP1, APC, and RASSF1 using MS-HRM can become the basis for the development of a laboratory test, the specifying diagnosis of prostate cancer.

Aberrant methylation is one of the driving forces of cancer genome development. Although the rate of methylation appears massively variable across the genome, it is mainly observed in histone modification, chromatin organization, DNA accessibility, or promoter sequence. Methylation of promoter sequence occurs mostly to cytosine nucleotides, which can affect transcription factors' binding affinities. In this study, we demonstrated that cytosine repeats (C types density), consisting of CC, CCC, CCCC, CCCCC, CCCCCC, CCCCCCC motifs and CpG islands density in 25 proto-oncogenes, tumor suppressor genes and control genes may play a role in the pathogenesis of acute myeloid leukemia. The promoter sequences were divided into a 100 nucleotide window from -500 to +100 nucleotides and 20 nucleotide window from -100 to +100. Each window is analyzed to find the higher C type and CpG islands density, which may cause the increased methylation in the promoter sequence. Our novel findings show that promoter sequence cytosine repeats and CpG density increase closer to transcription sites, especially just before and after the transcription start site (TSS). The results demonstrate that cytosine density increases while proto-oncogenes and TSG promoter sequences are closer to TSS 50.8% and 41.0% respectively, if (-500 to -200) and (-100 to +100) windows of the nucleotide sequences are compared. This proves that around TSS location has special nucleotide motifs and could be an important implication for our understanding of potential methylating locations in promoters.

The discovery of the enormous role methylated cytosine plays in regulating gene expression has led to the development of a variety of techniques for detecting cytosine modification. A majority of these techniques are geared towards analyzing genomic DNA, which is typically available in large quantities. The concentration of cell-free DNAs (cfDNA) extracted from biological fluids including plasma, saliva, tears, or urine is relatively low and their degree of the fragmentation is high. Moreover, for noninvasive diagnostics of cancer, methylation patterns must be studied in minor cancer-specific fractions of DNA molecules substantially diluted by excess unmethylated molecules. The above limitations complicate the application of traditional techniques for cfDNA methylation analysis. In this manuscript, we review the state-of-art analysis of cfDNA methylation, hydroxymethylation, and noncanonical methylation (outside of CpG islands). The review covers methodological approaches to studying individual CpGs and genomic loci, as well as techniques for the large-scale analysis of methylation.

In humans, DNA methylation at specific CpG sites can be used to estimate the ‘epigenetic clock’, a biomarker of aging and health. The mechanisms that regulate the aging epigenome and level of conservation are not entirely clear. We performed affinity-based enrichment with methyl-CpG binding domain protein followed by high-throughput sequencing (MBD-seq) to assay DNA methylation in mouse samples. Consistent with previous reports, aging is associated with increase in methylation at CpG islands that likely overlap regulatory regions of genes that have been implicated in cancers (e.g., C1ql3, Srd5a2 and Ptk7). The differentially methylated regions in mice have high sequence conservation in humans and the pattern of methylation is also largely conserved between the two species. Based on human ENCODE data, these sites are targeted by polycomb proteins, including EZH2. Chromatin immunoprecipitation confirmed that these regions interact with EZH2 in mice as well, and there may be reduction in EZH2 occupancy with age at C1ql3. This adds to the growing evidence that EZH2 is part of the protein machinery that shapes the aging epigenome. The conservation in both sequence and methylation patterns of the age-dependent CpGs indicate that the epigenetic clock is a fundamental feature of aging in mammals.

Obesity is a highly prevalent, chronic disorder that has been increasing in incidence in young patients. Both epigenetic and genetic aberrations may play a role in the pathogenesis of obesity. Therefore, in-depth epigenomic and genomic analyses will advance our understanding of the detailed molecular mechanisms underlying obesity and aid in the selection of potential biomarkers for obesity in youth. Here, we performed microarray-based DNA methylation and gene expression profiling of peripheral white blood cells obtained from six young, obese individuals and six healthy controls. We observed that the hierarchical clustering of DNA methylation, but not gene expression, clearly segregates the obese individuals from the controls, suggesting that the metabolic disturbance that occurs as a result of obesity at a young age may affect the DNA methylation of peripheral blood cells without accompanying transcriptional changes. To examine the genome-wide differences in the DNA methylation profiles of young obese and control individuals, we identified differentially methylated CpG sites and investigated their genomic and epigenomic contexts. The aberrant DNA methylation patterns in obese individuals can be summarized as relative gains and losses of DNA methylation in gene promoters and gene bodies, respectively. We also observed that the CpG islands of obese individuals are more susceptible to DNA methylation compared to controls. Our pilot study suggests that the genome-wide aberrant DNA methylation patterns of obese individuals may advance not only our understanding of the epigenomic pathogenesis but also early screening of obesity in youth.

Non LTR retrotransposons (EhLINEs and EhSINEs) occupy 11% of the Entamoeba histolytica genome. Since promoter DNA methylation at cytosines has been correlated with transcriptional silencing of transposable elements in model organisms we checked whether this was the case in EhLINE1. We located promoter activity in a 841bp fragment at 5'-end of this element by luciferase reporter assay. From RNAseq and RT-PCR analyses we selected a transcriptionally active and silent copy to study cytosine DNA methylation of the promoter region by bisulfite sequencing. None of the cytosines were methylated in either copy. Further, we looked at methylation status of a few selected cytosines in all 5'-intact EhLINE1 copies by single nucleotide incorporation opposite cytosine in bisulfite-treated DNA, where dGTP would be incorporated if the cytosine was methylated. Again we did not find evidence of cytosine methylation, indicating that expression status of this element was not correlated with promoter DNA methylation. To test for any role of cytosine methylation in transcriptional regulation of the E. histolytica Hsp70 gene in which the promoter is fully methylated under normal growth conditions, we checked methylation status and found that the promoter remained fully methylated during heat-shock as well, although transcription was greatly enhanced by heat-shock, showing that cytosine methylation is not a repressive mark for EhHsp70. Our data present direct evidence that promoter methylation, a common mode of transposon silencing, is unlikely to be involved in transcriptional regulation of EhLINE1, and reinforce the conclusion that promoter DNA methylation may not be a major contributor to transcriptional regulation in E. histolytica.

© 2017 The Author(s). Published with license by Taylor & Francis Group, LLC © 2017, © Florence K. Crary-Dooley, Mitchell E. Tam, Keith W. Dunaway, Irva Hertz-Picciotto, Rebecca J. Schmidt, and Janine M. LaSalle. DNA methylation is an epigenetic mark at the interface of genetic and environmental factors relevant to human disease. Quantitative assessments of global DNA methylation levels have therefore become important tools in epidemiology research, particularly for understanding effects of environmental exposures in complex diseases. Among the available methods of quantitative DNA methylation measurements, bisulfite sequencing is considered the gold standard, but whole-genome bisulfite sequencing (WGBS) has previously been considered too costly for epidemiology studies with high sample numbers. Pyrosequencing of repetitive sequences within bisulfite-treated DNA has been routinely used as a surrogate for global DNA methylation, but a comparison of pyrosequencing to WGBS for accuracy and reproducibility of methylation levels has not been performed. This study compared the global methylation levels measured from uniquely mappable (non-repetitive) WGBS sequences to pyrosequencing assays of several repeat sequences and repeat assay-matched WGBS data and determined uniquely mappable WGBS data to be the most reproducible and accurate measurement of global DNA methylation levels. We determined sources of variation in repetitive pyrosequencing assays to be PCR amplification bias, PCR primer selection bias in methylation levels of targeted sequences, and inherent variability in methylation levels of repeat sequences. Low-coverage, uniquely mappable WGBS showed the strongest correlation between replicates of all assays. By using multiplexing by indexed bar codes, the cost of WGBS can be lowered significantly to improve the accuracy of global DNA methylation assessments for human studies.

The chemokine receptor CXCR4 plays a key role in the bone marrow microenvironment maintenance and in the hematopoietic stem and progenitor cells migration. In addition, CXCR4 is expressed in a broad spectrum of solid tumors where its methylation state has been recently proposed as a biomarker for cancer prognosis. To evaluate methylation status of CXCR4 promoter we developed a sensitive, accurate, specific and cost-effective two-step PCR method that does not require any specific equipment other than a conventional real-time PCR instrument. The principle of the technique relies on a novel normalization strategy which allows the detection and quantification of small methylation differences among pre-amplified DNA samples deriving from low amount of starting material. In addition, the analysis of melting curve profiles of PCR products provides additional information about the methylation status of CpG sites in between the primers. Finally, the principle of this technique can potentially be adapted for the investigation of the methylation status of any other DNA region.

In this study, the methylation of mitochondrial genome in the immature testis of Chinese mitten crab Eriocheir sinensis of the Yangtze River system was determined for the first time using MeDIP-seq. Our methylated DNA fragments covered more than 99% of the mitochondrial genome in E. sinensis loaded from GenBank. There were 8 mutated bases and 42 SNPs in the crab mitochondrial genome. The methylation presented in all genes as well as in an A + T region, but less in intergenic regions in the mitochondrial genome. However, the level of methylation of most genes coding proteins and the A + T region were high. But, the majority of genes encoding tRNAs were hypomethylated, and both the rRNA genes also showed methylation of low or median frequency. Especially, the level of methylation of the intergenic regions is the lowest. Those features indicated that the methylation of DNA may play an important role in gene expressing regulation in the mitochondrial genome of immature testis in E. sinensis.

There is substantial evidence of impaired DNA repair activities in Alzheimer's disease (AD) neurons and peripheral tissues, inducing some investigators to speculate that this could partially result from promoter hypermethylation of DNA repair genes, resulting in gene silencing in those tissues. In the present study a screening cohort composed by late-onset AD (LOAD) patients and healthy matched controls was evaluated with a commercially available DNA methylation array for the assessment of the methylation levels of a panel of 22 genes involved in major DNA repair pathways in blood DNA. We then applied a cost-effective PCR based methylation-sensitive high-resolution melting (MS-HRM) technique, in order to evaluate the promoter methylation levels of the following DNA repair genes: OGG1, PARP1, MRE11A, BRCA1, MLH1, and MGMT. The analysis was performed in blood DNA from 56 LOAD patients and 55 matched controls, including the samples previously assessed with the DNA methylation array as validating samples. Both approaches revealed that all the investigated genes were largely hypomethylated in LOAD and control blood DNA, and no difference between groups was observed. Collectively, present data do not support an increased promoter methylation of some of the major DNA repair genes in blood DNA of AD patients.

Background/Aims: This study investigated the gene expression and DNA methylation of selected DNA repair genes (MBD4, TDG, MLH1, MLH3) and DNMT1 in human bladder cancer in the context of pathophysiological and prognostic significance. Methods: To determine the relationship between the gene expression pattern, global methylation and promoter methylation status, we performed real-time PCR to quantify the mRNA of selected genes in 50 samples of bladder cancer and adjacent non-cancerous tissue. The methylation status was analyzed by methylation-specific polymerase chain reaction (MSP) or digestion of genomic DNA with a methylation-sensitive restriction enzyme and PCR with gene-specific primers (MSRE-PCR). The global DNA methylation level was measured using the antibody-based 5-mC detection method. Results: The relative levels of mRNA for MBD4, MLH3, and MLH1 were decreased in 28% (14/50), 34% (17/50) and 36% (18/50) of tumor samples, respectively. The MBD4 mRNA expression was decreased in 46% of non-muscle invasive tumors (Ta/T1) compared with 11% found in muscle invasive tumors (T2-T4) (P

We report a new multiplexed strategy for the electrochemical detection of regional DNA methylation across multiple regions. Using the sequence dependent affinity of bisulfite treated DNA towards gold surfaces, the method integrates the high sensitivity of a micro-fabricated multiplex device comprising a microarray of gold electrodes, with the powerful multiplexing capability of multiplex-PCR. The synergy of this combination enables the monitoring of the methylation changes across several genomic regions simultaneously from as low as 500 pg μl−1 of DNA with no sequencing requirement.

The methylation status of sperm DNA differs between individual bulls. However, the relationship between methylation status and bull sperm parameters is not well elucidated. The present study investigated genome-wide methylation profiles at 450,000 CpG sites in bull spermatozoa by using a human DNA methylation microarray. Semen samples from three adult Japanese Black bulls with different in vitro fertilization (IVF) results and from a young Holstein bull through sexual maturation (at ages 10, 10.5, 15, and 25 months) were used for the analysis. The heatmap displaying the results of microarray analysis shows inter- and intra-individual differences in methylation profiles. After setting a cut-off of 0.2 for differences between ages (10, 10.5 vs. 15, 25 months) or between IVF results (developed to the blastocyst-stage, > 20% vs. < 10%), different methylation levels were detected at approximately 100 CpGs. We confirmed the different DNA methylation levels of CpG sites by using combined bisulfite restriction analysis (COBRA); five of the CpG sites reflected methylation levels similar to those detected by the microarray. One of the CpG sites was thought to reflect an age-related increase in methylation levels, which was confirmed by COBRA and bisulfite sequencing. However, the relationship between methylation status and IVF results could not be shown here. In conclusion, methylation profiles of individual and age-related alterations in bull spermatozoa can be revealed using a human microarray, and methylation changes in some CpG sites can be easily visualized using COBRA. Combined analysis of DNA methylation levels and sperm parameters could be considered an effective approach for assessing bull fertility in the future.

We report a single-cell bisulfite sequencing (scBS-seq) method that can be used to accurately measure DNA methylation at up to 48.4% of CpG sites. Embryonic stem cells grown in serum or in 2i medium displayed epigenetic heterogeneity, with '2i-like' cells present in serum culture. Integration of 12 individual mouse oocyte datasets largely recapitulated the whole DNA methylome, which makes scBS-seq a versatile tool to explore DNA methylation in rare cells and heterogeneous populations.

COHCAP (City of Hope CpG Island Analysis Pipeline) is an algorithm to analyze single-nucleotide resolution DNA methylation data produced by either an Illumina methylation array or targeted bisulfite sequencing. The goal of the COHCAP algorithm is to identify CpG islands that show a consistent pattern of methylation among CpG sites. COHCAP is currently the only DNA methylation package that provides integration with gene expression data to identify a subset of CpG islands that are most likely to regulate downstream gene expression, and it can generate lists of differentially methylated CpG islands with ∼50% concordance with gene expression from both cell line data and heterogeneous patient data. For example, this article describes known breast cancer biomarkers (such as estrogen receptor) with a negative correlation between DNA methylation and gene expression. COHCAP also provides visualization for quality control metrics, regions of differential methylation and correlation between methylation and gene expression. This software is freely available at https://sourceforge.net/projects/cohcap/.

DNA Methylation is an epigenetic phenomenon in which methyl groups are added to the cytosines, thereby altering the physio-chemical properties of the DNA region and influencing gene expression. Aberrant DNA methylation in a set of genes or across the genome results in many epigenetic diseases including cancer. In this paper, we use entropy to analyze the extent and distribution of DNA methylation in Tumor Suppressor Genes (TSG's) and Oncogenes related to a specific type of cancer (viz.) KIRC (Kidney-renal-clear-cell-carcinoma). We apply various mathematical transformations to enhance the different regions in DNA methylation distribution and compare the resultant entropies for healthy and tumor samples. We also obtain the sensitivity and specificity of classification for the different mathematical transformations. Our findings show that it is not just the measure of methylation, but the distribution of the methylation levels in the genes that are significant in cancer.

Little is known about the contributions of DNA methylation/demethylation to plant aging and senescence. We used Arabidopsis thaliana to study how increasing age of an annual plant species influences DNA methylation. Based on methylation-sensitive DNA fragmentation assay, it could be concluded that aging A. thaliana was accompanied by DNA demethylation. Bisulfite sequencing reveals that cytosine methylation within the Actin2 3’ untranslated region and internal transcribed spacer with 5.8S rRNA (ITS1-5.8SrRNA-ITS2) DNA regions decreased with A. thaliana growth and aging. We show that transcription of methyltransferase genes, chromomethyltransferase AtCMT3 and methyltransferse AtMETI, significantly decreased during development and aging of the A. thaliana plants, whereas expression of demethylase genes - repressor of silencing AtROS1, demeter AtDME, and demeter-like AtDML2 and AtDML3 - increased at least at some stages of plant development. The data obtained in the present study suggest that plant DNA regions may undergo demethylation during plant aging via reduction of DNA methylation processes and activation of active DNA demethylation.

With a purpose of accurate and simultaneous determination of DNA methylation from multiple loci in multiple samples, here, we are demonstrating a method to aid rapid DNA methylation detection of genomic sequences. Using genomic DNA of peripheral blood from 14 healthy individuals, DNA methylation in 465 CpG sites from 12 loci of genes (ADAM22, ATF2, BCR, CD83, CREBBP, IL12B, IL17RA, MAP2K2, RBM38, TGFBR2, TGFBR3, and WNT5A) was analysed by targeted next generation bisulfite sequencing. Analysed region for three genes, BCR, IL17RA and RBM38 showed an absolute mean DNA methylation of 25.6%, 89.2% and 38.9% respectively. Other nine gene loci were unmethylated and exhibited

// Roberto Tirado-Magallanes 1, 2 , Khadija Rebbani 1 , Ricky Lim 1 , Sriharsa Pradhan 3 , Touati Benoukraf 1 1 Cancer Science Institute of Singapore, National University of Singapore, 117599 Singapore, Singapore 2 Computational Systems Biology Team, Institut de Biologie de l’Ecole Normale Superieure (IBENS), INSERM, Ecole Normale Superieure, PSL Research University, 75005 Paris, France 3 New England Biolab Inc., Ipswich, MA 01938, USA Correspondence to: Touati Benoukraf, email: benoukraf@nus.edu.sg Keywords: DNA methylation, bisulfite sequencing, gene regulation, chromatin modeling, imprinting Received: August 11, 2016 Accepted: November 07, 2016 Published: November 24, 2016 ABSTRACT The combination of DNA bisulfite treatment with high-throughput sequencing technologies has enabled investigation of genome-wide DNA methylation at near base pair level resolution, far beyond that of the kilobase-long canonical CpG islands that initially revealed the biological relevance of this covalent DNA modification. The latest high-resolution studies have revealed a role for very punctual DNA methylation in chromatin plasticity, gene regulation and splicing. Here, we aim to outline the major biological consequences of DNA methylation recently discovered. We also discuss the necessity of tuning DNA methylation resolution into an adequate scale to ease the integration of the methylome information with other chromatin features and transcription events such as gene expression, nucleosome positioning, transcription factors binding dynamic, gene splicing and genomic imprinting. Finally, our review sheds light on DNA methylation heterogeneity in cell population and the different approaches used for its assessment, including the contribution of single cell DNA analysis technology.

SummaryWe analyzed the in silico purified DNA methylation signatures of 82 mantle cell lymphomas (MCL) in comparison with cell subpopulations spanning the entire B cell lineage. We identified two MCL subgroups, respectively carrying epigenetic imprints of germinal-center-inexperienced and germinal-center-experienced B cells, and we found that DNA methylation profiles during lymphomagenesis are largely influenced by the methylation dynamics in normal B cells. An integrative epigenomic approach revealed 10,504 differentially methylated regions in regulatory elements marked by H3K27ac in MCL primary cases, including a distant enhancer showing de novo looping to the MCL oncogene SOX11. Finally, we observed that the magnitude of DNA methylation changes per case is highly variable and serves as an independent prognostic factor for MCL outcome.

DNA methylation is globally reprogrammed after fertilization, and as a result, the parental genomes have similar DNA-methylation profiles after implantation except at the germline differentially methylated regions (gDMRs). We and others have previously shown that human blastocysts might contain thousands of transient maternally methylated gDMRs (transient mDMRs), whose maternal methylation is lost in embryonic tissues after implantation. In this study, we performed genome-wide allelic DNA methylation analyses of purified trophoblast cells from human placentas and, surprisingly, found that more than one-quarter of the transient-in-embryo mDMRs maintained their maternally biased DNA methylation. RNA-sequencing-based allelic expression analyses revealed that some of the placenta-specific mDMRs were associated with expression of imprinted genes (e.g., TIGAR, SLC4A7, PROSER2-AS1, and KLHDC10), and three imprinted gene clusters were identified. This approach also identified some X-linked gDMRs. Comparisons of the data with those from other mammals revealed that genomic imprinting in the placenta is highly variable. These findings highlight the incomplete erasure of germline DNA methylation in the human placenta; understanding this erasure is important for understanding normal placental development and the pathogenesis of developmental disorders with imprinting effects.

DNA methylation is essential for normal mammalian development and plays critical roles in various biological processes, including genomic imprinting, X-chromosome inactivation and repression of transposable elements. Although DNA methylation patterns are relatively stable in somatic cells, global reprogramming of DNA methylation occurs during mammalian preimplantation development. Advances in DNA methylation profiling techniques have been revealing the DNA methylation dynamics in mammalian embryos. Recently, we and other groups reported genome-scale DNA methylation analyses of human oocytes and preimplantation embryos, highlighting both the similarities and differences in the DNA methylation dynamics between humans and mice. In this review, we introduce the current knowledge of DNA methylation dynamics during early mammalian development. We also discuss the possibility of the application of genome-scale DNA methylation analysis techniques to human gametes and embryos for diagnostic purposes.

Aberrant gene expression in cancer is due in part to irregular patterns of DNA methylation and miRNA target gene down regulation. Using data from The Cancer Genome Atlas (TGCA), we investigated co-localized mRNA, miRNA and DNA methylation data across 15 cancer types, focusing on cases where evidence for direct regulation was strong. Restricting attention to regions where miRNA markers co-localize within a corresponding mRNA transcript, we checked expression data from 2839 samples for 354 mRNAs, 389 miRNAs and 13,809 DNA methylation probes for correlations greater than an absolute 0.6. We identified 32 genes, 34 miRNAs and 143 DNA methylation site probes comprising 180 "triplet combinations" that together provide evidence of co-localized genomic regulation in cancer. The five triplet combinations showing the highest prevalence across tissue types were found in four genes, HOXC5, PDE2A, SH3TC2 and TP63. Of the total 32 genes, eight among two tumor types (Kidney Renal Clear Cell Carcinoma (KIRC, 4) and Low Grade Glioma (LGG, 4)) were significantly associated with survival time (p 0.002). Together, the data presented in this paper provide evidence toward our primary hypothesis, that both genes and miRNAs strongly correlated with methylation level are more likely to be associated with cancer outcomes.

The initiation of embryonic gene expression in ascidian embryos appears to be tightly regulated by the number of DNA replication cycles. DNA methylation is thought to contribute to the clock mechanism that counts the rounds of DNA replication. We used mass spectrometry and whole genome bisulfite sequencing to characterize DNA methylation changes that occur in early developmental stages of the ascidian, Ciona intestinalis. We found that global DNA methylation in early Ciona development was static, and a base-wise comparison between the genomes of consecutive developmental stages found no DNA demethylation that was related to zygotic gene activation. Additionally, 5hmC was hardly detected by mass spectrometry in the developmental samples, suggesting a lack of demethylation mediated by ten eleven translocation (TET) methylcytosine dioxygenase in C. intestinalis. We conclude that DNA methylation is not involved in regulating DNA replication-dependent transcriptional activation.

We report genome-wide single-base resolution maps of methylated cytosines and transcriptome change in Cd-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between Cd-exposed and Cd-free rice genomes. There are 2320 non-redundant differentially methylated regions detected in the genome. RNA sequencing revealed 2092 DNA methylation-modified genes differentially expressed under Cd exposure. More genes were found hypermethylated than those hypomethylated in CG, CHH and CHG (where H is A, C or T) contexts in upstream, gene body and downstream regions. Many of the genes were involved in stress response, metal transport and transcription factors. Most of the DNA methylation-modified genes were transcriptionally altered under Cd stress. A subset of loss of function mutants defective in DNA methylation and histone modification activities was used to identify transcript abundance of selected genes. Compared with wide type, mutation of MET1 and DRM2 resulted in general lower transcript levels of the genes under Cd stress. Transcripts of OsIRO2, OsPR1b and Os09g02214 in drm2 were significantly reduced. A commonly used DNA methylation inhibitor 5-azacytidine was employed to investigate whether DNA demethylation affected physiological consequences. 5-azacytidine provision decreased general DNA methylation levels of selected genes, but promoted growth of rice seedlings and Cd accumulation in rice plant.

DNA methylation plays an important role in the regulation of gene expression. In normal cells, transposable elements that constitute approximately 45% of the human genome are highly methylated to silence their expression. In cancer cells, transposable elements are hypomethylated; therefore, global DNA methylation level is considered as a biomarker for cancer diagnostics. In this study, a homogeneous assay for measuring global DNA methylation level based on bioluminescence resonance energy transfer (BRET) was developed using methyl-CpG binding domain (MBD)-fused luciferase. In this assay, the MBD-luciferase recognizes methylated CpG, thus, BRET between the luciferase and fluorescent DNA intercalating dye is detected. We demonstrated that the BRET signal depended on the DNA methylation level of the target DNA. Moreover, the BRET signal was correlated with the LINE1 DNA methylation level on human genomic DNA, as determined by the bisulfite method. These results indicate that the global DNA methylation level of human genomic DNA could be detected simply by measuring the BRET signal.

Aims: We compared four common methods for measuring DNA methylation levels and recommended the most efficient method in terms of cost and coverage. Materials & methods: The DNA methylation status of liver and stomach tissues was profiled using four different methods, whole-genome bisulphite sequencing (WG-BS), targeted bisulphite sequencing (Targeted-BS), methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA immunoprecipitation bisulphite sequencing (MeDIP-BS). We calculated DNA methylation levels using each method and compared the results. Results: MeDIP-BS yielded the most similar DNA methylation profile to WG-BS, with 20 times less data, suggesting remarkable cost savings and coverage efficiency compared with the other methods. Conclusion: MeDIP-BS is a practical cost-effective method for analyzing whole-genome DNA methylation that is highly accurate at base-pair resolution.

Pan-cancer analysis can identify cell- and tissue-specific genomic loci and regions with underlying biological functions. Here we present an online curated DNA Methylation Annotation Knowledgebase, DMAK, which includes the pan-cancer analysis results for differentially-methylated loci and regions by the Reduced Representation Bisulfite Sequencing profiling technology. DMAK contains 3 modules of curated information and analysis results on 688,445 CpG sites across 19 cancer and embryonic stem cell lines from ENCODE, and further analysis of survival associations with clinical sources retrieved from TCGA. The knowledgebase covers all identified differentially-methylated CpG sites and regions of interest, further annotated genomic information, together with tumor suppressor genes information and calculated methylation level. DMAK provides meaningful clues for deriving functional association network and related clinical association results based on protein-coding genes, including tumor suppressor genes, identified from differentially methylated regions of interest. Thus DMAK constitutes a comprehensive reference source for the current epigenetic research and clinical study.

DNA methylation has been proposed as one of the promising biomarkers for cancer diagnosis. In this study, we developed a DNA methylation detection system utilizing G-quadruplex and i-motif-forming sequences that requires neither sodium bisulfite treatment nor methylated DNA ligands. We hypothesized that G-quadruplex and i-motif structures would be stabilized by DNA methylation and arrest DNA polymerase activity during quantitative polymerase chain reaction (qPCR). The PCR products from VEGF, RET G-quadruplex, and i-motif-forming sequences were used as templates and analyzed by qPCR. Our results indicated that the initial elongation efficiency of PCR decreased with increasing DNA methylation levels in the G-quadruplex and i-motif-forming sequences. Moreover, we demonstrated that the initial elongation efficiency of PCR decreased with increased DNA methylation of the VEGF region on genomic DNA. These results indicated that DNA methylation of the G-quadruplex and i-motif-forming sequences on genomic DNA can be detected by qPCR.

A next generation bisulfite sequencing (NGBS) was used to study rDNA promoter methylation in human brain using postmortem samples of the parietal cortex. Qualitative analysis of patterns of CpG methylation was performed at the individual rDNA unit level. CpG site-specific differences in methylation frequency were observed with the core promoter harboring three out of four most methylated CpGs. Moreover, there was an overall trend towards co-methylation for all possible pairs of 26 CpG sites. The hypermethylated CpGs from the core promoter were also most likely to be co-methylated. Finally, although rare, non-CpG (CpH) methylation was detected at several sites with one of them confirmed using the PspGI-qPCR assay. Similar trends were observed in samples from control individuals as well as patients suffering of Alzheimer's disease (AD), mild cognitive impairment (MCI) or ataxia telangiectasia (AT). Taken together, while some methyl-CpG sites including those in the core promoter may have relatively greater inhibitory effect on rRNA transcription, co-methylation at multiple sites may be required for full and/or long lasting silencing of human rDNA.