Genome-wide identification of 5-methylcytosine sites in bacterial genomes by high-throughput sequencing of MspJI restriction fragments.

Single-molecule Real-Time (SMRT) sequencing can easily identify sites of N6-methyladenine and N4-methylcytosine within DNA sequences, but similar identification of 5-methylcytosine sites is not as straightforward. In prokaryotic DNA, methylation typically occurs within specific sequence contexts, or motifs, that are a property of the methyltransferases that "write" these epigenetic marks. We present here a straightforward, cost-effective alternative to both SMRT and bisulfite sequencing for the determination of prokaryotic 5-methylcytosine methylation motifs. The method, called MFRE-Seq, relies on excision and isolation of fully methylated fragments of predictable size using MspJI-Family Restriction Enzymes (MFREs), which depend on the presence of 5-methylcytosine for cleavage. We demonstrate that MFRE-Seq is compatible with both Illumina and Ion Torrent sequencing platforms and requires only a digestion step and simple column purification of size-selected digest fragments prior to standard library preparation procedures. We applied MFRE-Seq to numerous bacterial and archaeal genomic DNA preparations and successfully confirmed known motifs and identified novel ones. This method should be a useful complement to existing methodologies for studying prokaryotic methylomes and characterizing the contributing methyltransferases.

Polymerization retardation isothermal amplification (PRIA): a strategy enables sensitively quantify genome-wide 5-methylcytosine oxides rapidly on handy instruments with nanoscale sample input.

The current methods for quantifying genome-wide 5-methylcytosine (5mC) oxides are still scarce, mostly restricted with two limitations: assay sensitivity is seriously compromised with cost, assay time and sample input; epigenetic information is irreproducible during polymerase chain reaction (PCR) amplification without bisulfite pretreatment. Here, we propose a novel Polymerization Retardation Isothermal Amplification (PRIA) strategy to directly amplify the minute differences between epigenetic bases and others by arranging DNA polymerase to repetitively pass large electron-withdrawing groups tagged 5mC-oxides. We demonstrate that low abundant 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxycytosine (5caC) in genomic DNA can be accurately quantified within 10 h with 100 ng sample input on a laboratory real-time quantitative PCR instrument, and even multiple samples can be analyzed simultaneously in microplates. The global levels of 5hmC and 5fC in mouse and human brain tissues, rat hippocampal neuronal tissue, mouse kidney tissue and mouse embryonic stem cells were quantified and the observations not only confirm the widespread presence of 5hmC and 5fC but also indicate their significant variation in different tissues and cells. The strategy is easily performed in almost all research and medical laboratories, and would provide the potential capability to other candidate modifications in nucleotides.

A Practical Guide to the Measurement and Analysis of DNA Methylation.

DNA methylation represents a fundamental epigenetic mark that is associated with transcriptional repression during development, maintenance of homeostasis, and disease. In addition to methylation-sensitive PCR and targeted deep-amplicon bisulfite sequencing to measure DNA methylation at defined genomic loci, numerous unsupervised techniques exist to quantify DNA methylation on a genome-wide scale, including affinity enrichment strategies and methods involving bisulfite conversion. Both affinity-enriched and bisulfite-converted DNA can serve as input material for array hybridization or sequencing using next-generation technologies. In this practical guide to the measurement and analysis of DNA methylation, the goal is to convey basic concepts in DNA methylation biology and explore genome-scale bisulfite sequencing as the current gold standard for assessment of DNA methylation. Bisulfite conversion chemistry and library preparation are discussed in addition to a bioinformatics approach to quality assessment, trimming, alignment, and methylation calling of individual cytosine residues. Bisulfite-converted DNA presents challenges for standard next-generation sequencing library preparation protocols and data-processing pipelines, but these challenges can be met with elegant solutions that leverage the power of high-performance computing systems. Quantification of DNA methylation, data visualization, statistical approaches to compare DNA methylation between sample groups, and examples of integrating DNA methylation data with other -omics data sets are also discussed. The reader is encouraged to use this article as a foundation to pursue advanced topics in DNA methylation measurement and data analysis, particularly the application of bioinformatics and computational biology principles to generate a deeper understanding of mechanisms linking DNA methylation to cellular function.

Global DNA demethylation as an epigenetic marker of human brain metastases.

Brain metastases are the most common intracranial tumors in adults. They usually originate from: lung, breast, renal cell and gastrointestinal cancers, as well as melanoma. Prognosis for brain metastases is still poor and classical treatment combining surgery and radiation therapy should be strongly supported with molecular approaches. However, their successful application depends on a deep understanding of not only genetic, but also epigenetic background of the disease. That will result in an earlier and more precise diagnosis, successful treatment, as well as individualized estimation of clinical outcomes and prognosis. It has already been shown that the epigenetic machinery plays a crucial role in cancer biology, development, and progression. Therefore, we decided to look for metastasis through changes in the most studied epigenetic mark, 5-methylcytosine (m5C) in DNA. We performed global analysis of the m5C contents in DNA isolated from the brain metastatic tumor tissue and peripheral blood samples of the same patients, using thin layer chromatography separation of radioactively labeled nucleotides. We found that the m5C level in DNA from brain metastases: changes in the broad range, overlaps with that of blood, and negatively correlates with the increasing tumor grade. Because the amount of m5C in tumor tissue and blood is almost identical, the genomic DNA methylation can be a useful marker for brain metastases detection and differentiation. Our research creates a scope for future studies on epigenetic mechanisms in neuro-oncology and can lead to development of new diagnostic methods in clinical practice.

Ultrasensitive electrochemiluminescence immunosensor for 5-hydroxymethylcytosine detection based on Fe3O4@SiO2 nanoparticles and PAMAM dendrimers.

An ultrasensitive sandwiched electrochemiluminescence (ECL) immunosensor was developed for 5-hydroxymethylcytosine (5hmC) detection in genomic DNA by using Fe3O4@SiO2 core-shell magnetic nanomaterial as a immobilization matrix for anti-5hmC antibody, PAMAM conjugated avidin and Ru(bpy)2(phen-5-NH2)(PF6)2 as signal amplification unit. Importantly, Fe3O4@SiO2 nanoparticles were verified to not only possess enormous surface for loading antibody by amido link, but also exhibit excellent bioactivity. With the dual signal amplification strategy, the ECL immunosensor showed wide detection range from 0.1 to 30nM with low detection limit of 0.047nM (S/N = 3). Based on the specific immunoreaction, the developed method also illustrated excellent detection selectivity. The fabricated immunosensor was also applied to detect the 5hmC in genomic DNA of cancer tissue, which indicated that the immunosensor possess potential applications in clinical detection.

Probing RNA Modification Status at Single-Nucleotide Resolution in Total RNA.

RNA modifications, with over one hundred known so far, are commonly proposed to fine-tune the structure and function of RNA. While modifications in rRNA and tRNA are used to modulate RNA folding and decoding properties, little is known about the function of internal modifications in mRNA/lncRNA, which includes N(6)-methyl adenosine (m(6)A), 5-methyl cytosine (m(5)C), 2'-O-methylated nucleotides (Nm), pseudouridine (Ψ), and possible others. Functional studies of mRNA/lncRNA modifications have been hindered by the lack of methods for their identification at single-nucleotide resolution. Challenges for the determination of mRNA/lncRNA modifications at single-nucleotide resolution are mainly due to the low abundance of mRNA/lncRNA. Traditional deep sequencing methods cannot identify mRNA/lncRNA modifications, such as m(6)A, m(5)C, Nm, and Ψ, because reverse transcriptase is insensitive to their presence in cDNA synthesis. Antibody-based approach enables the identification of m(6)A regions in mRNA/lncRNA, but currently at ~100 nucleotide resolution. Here, we describe a method that accurately identifies m(6)A position and modification fraction in human mRNA and lncRNAs at single-nucleotide resolution, termed "Site-specific Cleavage And Radioactive-labeling followed by Ligation-assisted Extraction and Thin-layer chromatography (SCARLET)." This method combines two previously established techniques, site-specific cleavage and splint ligation, to probe the RNA modification status at any mRNA/lncRNA site in the total RNA pool. SCARLET can potentially analyze any nucleotide that maintains Watson-Crick base pairing in the transcriptome and determine whether it contains m(6)A, m(5)C, Nm, Ψ, or other modifications yet to be discovered. Precise determination of the position and modification fraction of RNA modifications reveals crucial parameters for functional investigation of RNA modifications.

Experimental Approaches for Target Profiling of RNA Cytosine Methyltransferases.

RNA cytosine methyltransferases (m(5)C-RMTs) constitute an important class of RNA-modifying enzymes, methylating specific cytosines within particular RNA targets in both coding and noncoding RNAs. Almost all organisms express at least one m(5)C-RMT, and vertebrates often express different types or variants of m(5)C-RMTs in different cell types. Deletion or mutation of particular m(5)C-RMTs is connected to severe pathological manifestations ranging from developmental defects to infertility and mental retardation. Some m(5)C-RMTs show spatiotemporal patterns of expression and activity requiring careful experimental design for their analysis in order to capture their context-dependent targets. An essential step for understanding the functions of both the enzymes and the modified cytosines is defining the one-to-one connection between particular m(5)C-RMTs and their target cytosines. Recent technological and methodological advances have provided researchers with new tools to comprehensively explore RNA cytosine methylation and methyltransferases. Here, we describe three complementary approaches applicable for both discovery and validation of candidate target sites of specific m(5)C-RMTs.

RNA 5-Methylcytosine Analysis by Bisulfite Sequencing.

Cells have developed molecular machineries, which can chemically modify DNA and RNA nucleosides. One particular and chemically simple modification, (cytosine-5) methylation (m(5)C), has been detected both in RNA and DNA suggesting universal use of m(5)C for the function of these nucleotide polymers. m(5)C can be reproducibly mapped to abundant noncoding RNAs (transfer RNA, tRNA and ribosomal RNA, rRNA), and recently, also nonabundant RNAs (including mRNAs) have been reported to carry this modification. Quantification of m(5)C content in total RNA preparations indicates that a limited number of RNAs carry this modification and suggests specific functions for (cytosine-5) RNA methylation. What exactly is the biological function of m(5)C in RNA? Before attempting to address this question, m(5)C needs to be mapped specifically and reproducibly, preferably on a transcriptome-wide scale. To facilitate the detection of m(5)C in its sequence context, RNA bisulfite sequencing (RNA-BisSeq) has been developed. This method relies on the efficient chemical deamination of nonmethylated cytosine, which can be read out as single nucleotide polymorphism (nonmethylated cytosine as thymine vs. methylated cytosine as cytosine), when differentially comparing cDNA libraries to reference sequences after DNA sequencing. Here, the basic protocol of RNA-BisSeq, its current applications and limitations are described.

Detection and mapping of 5-methylcytosine and 5-hydroxymethylcytosine with nanopore MspA.

Precise and efficient mapping of epigenetic markers on DNA may become an important clinical tool for prediction and identification of ailments. Methylated CpG sites are involved in gene expression and are biomarkers for diseases such as cancer. Here, we use the engineered biological protein pore Mycobacterium smegmatis porin A (MspA) to detect and map 5-methylcytosine and 5-hydroxymethylcytosine within single strands of DNA. In this unique single-molecule tool, a phi29 DNA polymerase draws ssDNA through the pore in single-nucleotide steps, and the ion current through the pore is recorded. Comparing current levels generated with DNA containing methylated CpG sites to current levels obtained with unmethylated copies of the DNA reveals the precise location of methylated CpG sites. Hydroxymethylation is distinct from methylation and can also be mapped. With a single read, the detection efficiency in a quasirandom DNA strand is 97.5 ± 0.7% for methylation and 97 ± 0.9% for hydroxymethylation.

Error rates for nanopore discrimination among cytosine, methylcytosine, and hydroxymethylcytosine along individual DNA strands.

Cytosine, 5-methylcytosine, and 5-hydroxymethylcytosine were identified during translocation of single DNA template strands through a modified Mycobacterium smegmatis porin A (M2MspA) nanopore under control of phi29 DNA polymerase. This identification was based on three consecutive ionic current states that correspond to passage of modified or unmodified CG dinucleotides and their immediate neighbors through the nanopore limiting aperture. To establish quality scores for these calls, we examined ~3,300 translocation events for 48 distinct DNA constructs. Each experiment analyzed a mixture of cytosine-, 5-methylcytosine-, and 5-hydroxymethylcytosine-bearing DNA strands that contained a marker that independently established the correct cytosine methylation status at the target CG of each molecule tested. To calculate error rates for these calls, we established decision boundaries using a variety of machine-learning methods. These error rates depended upon the identity of the bases immediately 5' and 3' of the targeted CG dinucleotide, and ranged from 1.7% to 12.2% for a single-pass read. We estimate that Q40 values (0.01% error rates) for methylation status calls could be achieved by reading single molecules 5-19 times depending upon sequence context.

Enzymatic analysis of Tet proteins: key enzymes in the metabolism of DNA methylation.

One of the most exciting recent advances in the epigenetic field is the discovery that 5-methylcytosine (5mC) in DNA can be iteratively oxidized by a family of proteins known as Tet proteins to generate 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). These 5mC derivatives can be further processed by thymine-DNA glycosylase (TDG) followed by base excision repair or by replication-dependent dilution leading to DNA demethylation. Given the similarity between 5mC and its oxidation derivatives, many of the conventional techniques used for 5mC analysis cannot distinguish between 5mC and 5hmC/5fC/5caC. Here, we describe 2D-TLC and mass spectrometry methods that we have successfully used in differentiating 5mC from its oxidative derivatives as well as in characterizing the enzymatic activity of Tet proteins both in vitro and in vivo.

Bisulfite methylation profiling of large genomes.

Bisulfite conversion of genomic DNA differentiates cytosines from 5-methylcytosines and, thus, identifies DNA methylation patterns at the single-base level. Here, we review recent developments incorporating high-throughput sequencing of bisulfite-converted DNA for target-specific analyses and genome-wide mapping of plant and mammalian methylomes. These developments include the analysis of human embryonic stem cell and fetal fibroblast methylomes at single-base resolution, which supports the presence of non-CG DNA methylation in wild-type embryonic stem cells and induced pluripotent stem cells. New developments in nanopore sequencing technologies may lead to directly detecting 5-methylcytosine independently of bisulfite conversion, but the current accuracy of this approach remains a limitation. Furthermore, recent investigations detecting 5-hydroxymethylcytosine within mammalian DNA may add yet another level of complexity to the epigenetic code of the methylome.

Capillary electrophoresis-laser induced fluorescence analysis of endogenous damage in mitochondrial and genomic DNA.

Reactive oxygen molecules are formed in vivo as by-products of normal aerobic metabolism. All organisms dependent on oxygen are inevitably exposed to these species so that DNA damage can occur in both genomic and mitochondrial DNA (mtDNA). In order to determine endogenous DNA damage we have developed an analytical method that involves the isolation and hydrolysis of genomic DNA or mtDNA, the labeling of modified and unmodified nucleotides and micellar electrokinetic chromatography with laser-induced fluorescence detection. With this method we have found etheno-adenine, thymine glycol, uracil, hypoxanthine, and 5-methylcytosine. These were identified by the addition of internal standards to the genomic or mtDNA. There are a large number of other signals in the electropherograms of mtDNA that we have never found in genomic DNA analysis because they are at lower concentration in the genome. In the DNA of untreated patients with chronic lymphocytic leukemia (CLL), uracil and high levels of etheno-adenine were found, which can be explained by antioxidant enzyme alterations and oxidative stress in the CLL lymphocytes.