Cloning, characterization and expression of the gene coding for a cytosine-5-DNA methyltransferase recognizing GpC

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

One-pot trimodal mapping of unmethylated, hydroxymethylated, and open chromatin sites unveils distinctive 5hmC roles at dynamic chromatin loci

We present a method, named Mx-TOP, for profiling of three epigenetic regulatory layers-chromatin accessibility, general DNA modification, and DNA hydroxymethylation-from a single library. The approach is based on chemo-enzymatic covalent tagging of unmodified CG sites and hydroxymethylated cytosine (5hmC) along with GC sites in chromatin, which are then mapped using tag-selective base-resolution TOP-seq sequencing. Our in-depth validation of the approach revealed its sensitivity and informativity in evaluating chromatin accessibility and DNA modification interactions that drive transcriptional regulation. We employed the technology in a study of chromatin and DNA demethylation dynamics during in vitro neuronal differentiation. The study highlighted the involvement of gene body 5hmC in modulating an extensive decoupling between promoter accessibility and transcription. The importance of 5hmC in chromatin remodeling was further demonstrated by the observed resistance of the developmentally acquired open loci to the global 5hmC erasure in neuronal progenitors.

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Investigating pioneer factor activity and its coordination with chromatin remodelers using integrated synthetic oligo assay

Chromatin accessibility is regulated by pioneer factors (PFs) and chromatin remodelers (CRs). Here, we present a protocol, based on integrated synthetic oligonucleotide libraries in yeast, to systematically interrogate the nucleosome-displacing activities of PFs and their coordination with CRs. We describe steps for designing oligonucleotide sequences, constructing yeast libraries, measuring nucleosome configurations, and data analyses. This approach potentially can be adapted for use in higher eukaryotes to investigate the activities of many types of chromatin-associated factors. For complete details on the use and execution of this protocol, please refer to Yan et al.,1 and Chen et al.2.

Controllable assembly of dendritic DNA nanostructures for ultrasensitive detection of METTL3-METTL14 m6A methyltransferase activity in cancer cells and human breast tissues

N⁶-Methyladenosine (m⁶A) is a reversible chemical modification in eukaryotic messenger RNAs and long noncoding RNAs. The aberrant expression of RNA methyltransferase METTL3-METTL14 complex may change the m⁶A methylation level and cause multiple diseases including cancers. The conventional METTL3-METTL14 assays commonly suffer from time-consuming procedures and poor sensitivity. Herein, we develop a controllable amplification machinery based on MazF-activated terminal deoxynucleotidyl transferase (TdT)-assisted dendritic DNA structure assembly for ultrasensitive detection of METTL3-METTL14 complex activity in cancer cells and breast tissues. The presence of METTL3-METTL14 complex catalyzes the formation of m⁶A in detection probe, effectively preventing the cleavage of methylated detection probes by MazF. The methylated detection probes with 3'-OH termini can function as the primers for template-free polymerization catalyzed by TdT on magnetic beads (MBs), producing long chains of poly-thymidine (poly-T) sequences. Then poly-T sequences hybridize with signal probes that contain poly-adenine (poly-A) sequence, inducing TdT-mediated polymerization and the subsequent hybridization with more poly-A signal probes for generating dendritic DNA nanostructures assembled on MBs. After magnetic separation and elevated temperature treatment, the signal probes are disassembled from MBs to generate a high fluorescence signal. This method possesses excellent specificity and high sensitivity with a limit of detection (LOD) of 2.61 × 10^-15 M, and it can accurately quantify cellular METTL3-METTL14 complex at single-cell level. Furthermore, it can screen inhibitors, evaluate kinetic parameters, and discriminate breast cancer tissues from normal tissues.

Identification, expression, and purification of DNA cytosine 5-methyltransferases with short recognition sequences

Background

DNA methyltransferases (MTases) are enzymes that induce methylation, one of the representative epigenetic modifications of DNA, and are also useful tools for analyzing epigenomes. However, regarding DNA cytosine 5-methylation, MTases identified so far have drawbacks in that their recognition sequences overlap with those for intrinsic DNA methylation in mammalian cells and/or that the recognition sequence is too long for fine epigenetic mapping. To identify MTases with short recognition sequences that never overlap with the CG dinucleotide, we systematically investigated the 25 candidate enzymes identified using a database search, which showed high similarity to known cytosine 5-MTases recognizing short sequences.

Results

We identified MTases with six new recognition sequences, including TCTG, CC, CNG, TCG, GCY, and GGCA. Because the recognition sequence never overlapped with the CG dinucleotide, MTases recognizing the CC dinucleotide were promising.

Conclusions

In the current study, we established a procedure for producing active CC-methylating MTases and applied it to nucleosome occupancy and methylome sequencing to prove the usefulness of the enzyme for fine epigenetic mapping. MTases that never overlap with CG dinucleotides would allow us to profile multiple epigenomes simultaneously.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-022-00765-3.

Partitioned usage of chromatin remodelers by nucleosome-displacing factors

Nucleosome-displacing-factors (NDFs) in yeast, similar to pioneer factors in higher eukaryotes, can open closed chromatin and generate nucleosome-depleted regions (NDRs). NDRs in yeast are also affected by ATP-dependent chromatin remodelers (CRs). However, how NDFs and CRs coordinate in nucleosome invasion and NDR formation is still unclear. Here, we design a high-throughput method to systematically study the interplay between NDFs and CRs. By combining an integrated synthetic oligonucleotide library with DNA methyltransferase-based, single-molecule nucleosome mapping, we measure the impact of CRs on NDRs generated by individual NDFs. We find that CRs are dispensable for nucleosome invasion by NDFs, and they function downstream of NDF binding to modulate the NDR length. A few CRs show high specificity toward certain NDFs; however, in most cases, CRs are recruited in a factor-nonspecific and NDR length-dependent manner. Overall, our study provides a framework to investigate how NDFs and CRs cooperate to regulate chromatin opening.

Chromatin accessibility in yeast is regulated by nucleosome-displacing-factors (NDFs) and chromatin remodelers (CRs). Chen et al. show that NDFs first invade into nucleosomes and then recruit CRs to modulate the NDR length. NDF-specific and NDR length-dependent recruitment of CRs allow partitioned usage of CRs by NDFs.

Polycomb-lamina antagonism partitions heterochromatin at the nuclear periphery

The genome can be divided into two spatially segregated compartments, A and B, which partition active and inactive chromatin states. While constitutive heterochromatin is predominantly located within the B compartment near the nuclear lamina, facultative heterochromatin marked by H3K27me3 spans both compartments. How epigenetic modifications, compartmentalization, and lamina association collectively maintain heterochromatin architecture remains unclear. Here we develop Lamina-Inducible Methylation and Hi-C (LIMe-Hi-C) to jointly measure chromosome conformation, DNA methylation, and lamina positioning. Through LIMe-Hi-C, we identify topologically distinct sub-compartments with high levels of H3K27me3 and differing degrees of lamina association. Inhibition of Polycomb repressive complex 2 (PRC2) reveals that H3K27me3 is essential for sub-compartment segregation. Unexpectedly, PRC2 inhibition promotes lamina association and constitutive heterochromatin spreading into H3K27me3-marked B sub-compartment regions. Consistent with this repositioning, genes originally marked with H3K27me3 in the B compartment, but not the A compartment, remain largely repressed, suggesting that constitutive heterochromatin spreading can compensate for H3K27me3 loss at a transcriptional level. These findings demonstrate that Polycomb sub-compartments and their antagonism with lamina association are fundamental features of genome structure. More broadly, by jointly measuring nuclear position and Hi-C contacts, our study demonstrates how compartmentalization and lamina association represent distinct but interdependent modes of heterochromatin regulation.

Construction of an APE1-Mediated Cascade Signal Amplification Platform for Homogeneously Sensitive and Rapid Measurement of DNA Methyltransferase in Escherichia coli Cells

DNA methylation is an essential genomic epigenetic behavior in both eukaryotes and prokaryotes. Deregulation of DNA methyltransferase (Dam MTase) can change the DNA methylation level and cause various diseases. Herein, we develop an apurinic/apyrimidinic endonuclease 1 (APE1)-mediated cascade signal amplification platform for homogeneously sensitive and rapid measurement of Dam MTase in Escherichia coli cells. This assay involves a partial double-stranded DNA (dsDNA) substrate and two hairpin signal probes (HP1 and HP2) that are modified with Cy5 and BHQ2 at two ends, respectively. When Dam MTase is present, it methylates the dsDNA substrate, and subsequently, endonuclease DpnI cleaves the methylated substrate, yielding trigger probe 1. Hybridization of trigger probe 1 with HP1 forms a partial dsDNA containing an apurinic/apyrimidinic (AP) site, which is cleaved by APE1 to induce the cyclic cleavage of HP1 and the production of abundant trigger probe 2. Subsequent hybridization of trigger probe 2 with HP2 forms a partial dsDNA with an AP site, inducing the cyclic cleavage of HP2 by APE1. Consequently, cyclic cleavage of HP1 and HP2 induces the generation of abundant Cy5 molecules, which are easily measured by single-molecule imaging. This assay can be performed homogeneously and rapidly within 2 h, which is the shortest among the reported amplification-based assays. Moreover, it exhibits good selectivity and high sensitivity, and it can discriminate Dam MTase from other enzymes and screen inhibitors. Importantly, it can accurately measure the Dam MTase activity in serum and E. coli cells, with promising applications in clinical diagnosis and drug discovery.

Phage-encoded ten-eleven translocation dioxygenase (TET) is active in C5-cytosine hypermodification in DNA

Chemical tailoring of canonical bases expands the functionality of DNA in the same manner that posttranscriptional and -translational modifications enhance functional diversity in RNA and proteins. We describe the activities of ten-eleven translocation dioxygenase (TET)-like iron(II)- and 2-oxo-glutarate–dependent 5mC dioxygenases that are encoded by several bacteriophages to enable hypermodification of C5-methyl cytosine bases in their DNA. Phage TETs act on methylation marks deposited within GpC sequences by functionally-associated cytosine 5-methyltransferases. The hydroxymethyl groups installed are further elaborated by tailoring enzymes, thereby decorating the phage DNA with diverse, complex modifications. These modifications are predicted to have protective roles against host defenses during viral infection.

TET/JBP (ten-eleven translocation/base J binding protein) enzymes are iron(II)- and 2-oxo-glutarate–dependent dioxygenases that are found in all kingdoms of life and oxidize 5-methylpyrimidines on the polynucleotide level. Despite their prevalence, few examples have been biochemically characterized. Among those studied are the metazoan TET enzymes that oxidize 5-methylcytosine in DNA to hydroxy, formyl, and carboxy forms and the euglenozoa JBP dioxygenases that oxidize thymine in the first step of base J biosynthesis. Both enzymes have roles in epigenetic regulation. It has been hypothesized that all TET/JBPs have their ancestral origins in bacteriophages, but only eukaryotic orthologs have been described. Here we demonstrate the 5mC-dioxygenase activity of several phage TETs encoded within viral metagenomes. The clustering of these TETs in a phylogenetic tree correlates with the sequence specificity of their genomically cooccurring cytosine C5-methyltransferases, which install the methyl groups upon which TETs operate. The phage TETs favor Gp5mC dinucleotides over the 5mCpG sites targeted by the eukaryotic TETs and are found within gene clusters specifying complex cytosine modifications that may be important for DNA packaging and evasion of host restriction.

Methylscaper: an R/shiny app for joint visualization of DNA methylation and nucleosome occupancy in single-molecule and single-cell data

SUMMARY

Differential DNA methylation and chromatin accessibility are associated with disease development, particularly cancer. Methods that allow profiling of these epigenetic mechanisms in the same reaction and at the single-molecule or single-cell level continue to emerge. However, a challenge lies in jointly visualizing and analyzing the heterogeneous nature of the data and extracting regulatory insight. Here, we present methylscaper, a visualization framework for simultaneous analysis of DNA methylation and chromatin accessibility landscapes. Methylscaper implements a weighted principal component analysis that orders DNA molecules, each providing a record of the chromatin state of one epiallele, and reveals patterns of nucleosome positioning, transcription factor occupancy, and DNA methylation. We demonstrate methylscaper's utility on a long-read, single-molecule methyltransferase accessibility protocol for individual templates (MAPit-BGS) dataset and a single-cell nucleosome, methylation, and transcription sequencing (scNMT-seq) dataset. In comparison to other procedures, methylscaper is able to readily identify chromatin features that are biologically relevant to transcriptional status while scaling to larger datasets.

AVAILABILITY AND IMPLEMENTATION

Methylscaper, is implemented in R (version > 4.1) and available on Bioconductor: https://bioconductor.org/packages/methylscaper/, GitHub: https://github.com/rhondabacher/methylscaper/, and Web: https://methylscaper.com.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

Chloroviruses

Chloroviruses are large dsDNA, plaque-forming viruses that infect certain chlorella-like green algae; the algae are normally mutualistic endosymbionts of protists and metazoans and are often referred to as zoochlorellae. The viruses are ubiquitous in inland aqueous environments throughout the world and occasionally single types reach titers of thousands of plaque-forming units per ml of native water. The viruses are icosahedral in shape with a spike structure located at one of the vertices. They contain an internal membrane that is required for infectivity. The viral genomes are 290 to 370 kb in size, which encode up to 16 tRNAs and 330 to ~415 proteins, including many not previously seen in viruses. Examples include genes encoding DNA restriction and modification enzymes, hyaluronan and chitin biosynthetic enzymes, polyamine biosynthetic enzymes, ion channel and transport proteins, and enzymes involved in the glycan synthesis of the virus major capsid glycoproteins. The proteins encoded by many of these viruses are often the smallest or among the smallest proteins of their class. Consequently, some of the viral proteins are the subject of intensive biochemical and structural investigation.

Single-ribonucleotide repair-mediated ligation-dependent cycling signal amplification for sensitive and specific detection of DNA methyltransferase

Specific and sensitive detection of DNA MTase activity can be achieved by a single-ribonucleotide repair-mediated ligation-dependent cycling signal amplification approach.

DNA methylation is a predominant epigenetic modification that plays crucial roles in various cellular processes. DNA methyltransferase (MTase) is responsible for DNA methylation, and its dysregulation may induce aberrant methylation patterns that are closely related to cancers. Conventional methods for DNA MTase assay are usually cumbersome and laborious with poor sensitivity. Alternatively, some signal amplification strategies are employed to improve the sensitivity, but they suffer from poor specificity and consequently limited sensitivity due to the nonspecific amplification. Herein, we develop for the first time a new fluorescence method to specifically and sensitively detect DNA MTase activity on the basis of single-ribonucleotide repair-mediated ligation-dependent cycling signal amplification. In the presence of DNA MTase, the hairpin substrate is methylated and cleaved by endonuclease Dpn I, releasing a 24-nt cleavage product. The 24-nt cleavage product may function as a primer and adjacently hybridize with the ligation probes (LP1 and LP2) to form the template (LP1–LP2) for strand displacement amplification (SDA), initiating the single-ribonucleotide repair-mediated cyclic ligation-dependent SDA to produce a large number of reporter probes. The reporter probe can subsequently hybridize with the signal probe that is modified with FAM and BHQ1 to form a stable double-stranded DNA (dsDNA) duplex with a ribonucleotide mismatch. Ribonuclease HII (RNase HII) can excise the single ribonucleotide, resulting in the cyclic cleavage of signal probes and the generation of an enhanced fluorescence signal. Taking advantage of the high specificity of RNase HII-catalyzed single-ribonucleotide excision and the high amplification efficiency of cyclic ligation-dependent SDA, this assay exhibits the highest sensitivity reported so far with a detection limit of 4.8 × 10^–6 U mL^–1 and a large dynamic range of 5 orders of magnitude. Moreover, this method can be used for the discrimination of Dam MTase from other DNA MTases, the accurate quantification of Dam MTase activity in E. coli cells, and the screening of Dam MTase inhibitors, providing a new paradigm for biomedical research and clinical diagnosis.

Nucleosome Occupancy and Methylome Sequencing (NOMe-seq)

Various methodologies are available to interrogate specific components of epigenetic mechanisms such as DNA methylation or nucleosome occupancy at both the locus-specific and the genome-wide level. It has become increasingly clear, however, that comprehension of the functional interactions between epigenetic mechanisms is critical for understanding how cellular transcription programs are regulated or deregulated during normal and disease development. The Nucleosome Occupancy and Methylome sequencing (NOMe-seq) assay allows us to directly measure the relationship between DNA methylation and nucleosome occupancy by taking advantage of the methyltransferase M.CviPI, which methylates unprotected GpC dinucleotides to create a footprint of chromatin accessibility. This assay generates dual nucleosome occupancy and DNA methylation information at a single-DNA molecule resolution using as little as 200,000 cells and in as short as 15 min reaction time. DNA methylation levels and nucleosome occupancy status of genomic regions of interest can be subsequently interrogated by cloning PCR-amplified bisulfite DNA and sequencing individual clones. Alternatively, NOMe-seq can be combined with next-generation sequencing in order to generate an integrated global map of DNA methylation and nucleosome occupancy, which allows for comprehensive examination as to how these epigenetic components correlate with each other.

Defining Regulatory Elements in the Human Genome Using Nucleosome Occupancy and Methylome Sequencing (NOMe-Seq)

NOMe-seq (nucleosome occupancy and methylome sequencing) identifies nucleosome-depleted regions that correspond to promoters, enhancers, and insulators. The NOMe-seq method is based on the treatment of chromatin with the M.CviPI methyltransferase, which methylates GpC dinucleotides that are not protected by nucleosomes or other proteins that are tightly bound to the chromatin (GpC^m does not occur in the human genome and therefore there is no endogenous background of GpC^m). Following bisulfite treatment of the M.CviPI-methylated chromatin (which converts unmethylated Cs to Ts and thus allows the distinction of GpC from GpC^m) and subsequent genomic sequencing, nucleosome-depleted regions can be ascertained on a genome-wide scale. The bisulfite treatment also allows the distinction of CpG from C^mpG (most endogenous methylation occurs at CpG dinucleotides) and thus the endogenous methylation status of the genome can also be obtained in the same sequencing reaction. Importantly, open chromatin is expected to have high levels of GpC^m but low levels of C^mpG; thus, each of the two separate methylation analyses serve as independent (but opposite) measures which provide matching chromatin designations for each regulatory element.NOMe-seq has advantages over ChIP-seq for identification of regulatory elements because it is not reliant upon knowing the exact modifications on the surrounding nucleosomes. Also, NOMe-seq has advantages over DHS (DNase hypersensitive site)-seq, FAIRE (Formaldehyde-Assisted Isolation of Regulatory Elements)-seq, and ATAC (Assay for Transposase-Accessible Chromatin)-seq because it also gives positioning information for several nucleosomes on either side of each open regulatory element. Here, we provide a detailed protocol for NOMe-seq that begins with the isolation of chromatin, followed by methylation of GpCs with M.CviPI and treatment with bisulfite, and ending with the creation of next generation sequencing libraries. We also include sequencing QC analysis metrics and bioinformatics steps that can be used to identify nucleosome-depleted regions throughout the genome.

CAME: identification of chromatin accessibility from nucleosome occupancy and methylome sequencing

Motivation

Chromatin accessibility plays a key role in epigenetic regulation of gene activation and silencing. Open chromatin regions allow regulatory elements such as transcription factors and polymerases to bind for gene expression while closed chromatin regions prevent the activity of transcriptional machinery. Recently, Methyltransferase Accessibility Protocol for individual templates-Bisulfite Genome Sequencing (MAPit-BGS) and nucleosome occupancy and methylome sequencing (NOMe-seq) have been developed for simultaneously profiling chromatin accessibility and DNA methylation on single molecules. Therefore, there is a great demand in developing computational methods to identify chromatin accessibility from MAPit-BGS and NOMe-seq.

Results

In this article, we present CAME (Chromatin Accessibility and Methylation), a seed-extension based approach that identifies chromatin accessibility from NOMe-seq. The efficiency and effectiveness of CAME were demonstrated through comparisons with other existing techniques on both simulated and real data, and the results show that our method not only can precisely identify chromatin accessibility but also outperforms other methods.

Availability and Implementation

CAME is implemented in java and the program is freely available online at http://sourceforge.net/projects/came/.

Contacts

jechoi@gru.edu or khryu@dblab.chungbuk.ac.kr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Experimental mitochondria-targeted DNA methylation identifies GpC methylation, not CpG methylation, as potential regulator of mitochondrial gene expression

Like the nucleus, mitochondria contain their own DNA and recent reports provide accumulating evidence that also the mitochondrial DNA (mtDNA) is subjective to DNA methylation. This evidence includes the demonstration of mitochondria-localised DNA methyltransferases and demethylases, and the detection of mtDNA methylation as well as hydroxymethylation. Importantly, differential mtDNA methylation has been linked to aging and diseases, including cancer and diabetes. However, functionality of mtDNA methylation has not been demonstrated. Therefore, we targeted DNA methylating enzymes (modifying cytosine in the CpG or GpC context) to the mtDNA. Unexpectedly, mtDNA gene expression remained unchanged upon induction of CpG mtDNA methylation, whereas induction of C-methylation in the GpC context decreased mtDNA gene expression. Intriguingly, in the latter case, the three mtDNA promoters were differentially affected in each cell line, while cellular function seemed undisturbed. In conclusion, this is the first study which directly addresses the potential functionality of mtDNA methylation. Giving the important role of mitochondria in health and disease, unravelling the impact of mtDNA methylation adds to our understanding of the role of mitochondria in physiological and pathophysiological processes.

Methylation-assisted bisulfite sequencing to simultaneously map 5fC and 5caC on a genome-wide scale for DNA demethylation analysis

Active DNA demethylation is mediated by ten-eleven translocation (TET) proteins that progressively oxidize 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). We have developed a methylation-assisted bisulfite sequencing (MAB-seq) method that enables direct genome-scale mapping and quantification of 5fC and 5caC marks together at single-base resolution. In bisulfite sequencing (BS), unmethylated cytosine residues (Cs), 5fCs and 5caCs, are converted to uracil and cannot be discriminated from each other. The pretreatment of the DNA with the CpG methylation enzyme M.SssI, which converts only the Cs to 5mCs, protects Cs but not 5fCs and 5caCs, which enables direct detection of 5fCs and 5caCs as uracils. Here we also describe an adapted version of the protocol to perform reduced-representation MAB-seq (RRMAB-seq) that provides increased coverage on CpG-rich regions, thus reducing the execution costs and increasing the feasibility of the technique. The main advantage of MAB-seq is to reduce the number of chemical/enzymatic DNA treatments required before bisulfite treatment and to avoid the need for prohibitive sequencing coverage, thus making it more reliable and affordable than subtractive approaches. The method presented here is the ideal tool for studying DNA demethylation dynamics in any biological system. Overall timing is ∼3 d for library preparation.

High Fractional Occupancy of a Tandem Maf Recognition Element and Its Role in Long-Range β-Globin Gene Regulation

Enhancers and promoters assemble protein complexes that ultimately regulate the recruitment and activity of RNA polymerases. Previous work has shown that at least some enhancers form stable protein complexes, leading to the formation of enhanceosomes. We analyzed protein-DNA interactions in the murine β-globin gene locus using the methyltransferase accessibility protocol for individual templates (MAPit). The data show that a tandem Maf recognition element (MARE) in locus control region (LCR) hypersensitive site 2 (HS2) reveals a remarkably high degree of occupancy during differentiation of mouse erythroleukemia cells. Most of the other transcription factor binding sites in LCR HS2 or in the adult β-globin gene promoter regions exhibit low fractional occupancy, suggesting highly dynamic protein-DNA interactions. Targeting of an artificial zinc finger DNA-binding domain (ZF-DBD) to the HS2 tandem MARE caused a reduction in the association of MARE-binding proteins and transcription complexes at LCR HS2 and the adult βmajor-globin gene promoter but did not affect expression of the βminor-globin gene. The data demonstrate that a stable MARE-associated footprint in LCR HS2 is important for the recruitment of transcription complexes to the adult βmajor-globin gene promoter during erythroid cell differentiation.

Expression of exogenous DNA methyltransferases: application in molecular and cell biology

DNA methyltransferases might be used as powerful tools for studies in molecular and cell biology due to their ability to recognize and modify nitrogen bases in specific sequences of the genome. Methylation of the eukaryotic genome using exogenous DNA methyltransferases appears to be a promising approach for studies on chromatin structure. Currently, the development of new methods for targeted methylation of specific genetic loci using DNA methyltransferases fused with DNA-binding proteins is especially interesting. In the present review, expression of exogenous DNA methyltransferase for purposes of in vivo analysis of the functional chromatin structure along with investigation of the functional role of DNA methylation in cell processes are discussed, as well as future prospects for application of DNA methyltransferases in epigenetic therapy and in plant selection.

Human lymphoid translocation fragile zones are hypomethylated and have accessible chromatin

Chromosomal translocations are a hallmark of hematopoietic malignancies. CG motifs within translocation fragile zones (typically 20 to 600 bp in size) are prone to chromosomal translocation in lymphomas. Here we demonstrate that the CG motifs in human translocation fragile zones are hypomethylated relative to the adjacent DNA. Using a methyltransferase footprinting assay on isolated nuclei (in vitro), we find that the chromatin at these fragile zones is accessible. We also examined in vivo accessibility using cellular expression of a prokaryotic methylase. Based on this assay, which measures accessibility over a much longer time interval than is possible with in vitro methods, these fragile zones were found to be more accessible than the adjacent DNA. Because DNA within the fragile zones can be methylated by both cellular and exogenous methyltransferases, the fragile zones are predominantly in a duplex DNA conformation. These observations permit more-refined models for why these zones are 100- to 1,000-fold more prone to undergo chromosomal translocation than the adjacent regions.

Integrated DNA methylation and chromatin structural analysis at single-molecule resolution

Chromatin limits the accessibility of DNA to trans-acting factors in transcription, replication, and repair. Although transcriptional variation between cells in a population may contribute to survival and disease, most assays of chromatin structure recover only population averages. We have developed DNA methyltransferases (MTases) as probing agents of DNA accessibility in chromatin, either expressed in vivo in budding yeast or as recombinant enzymatic probes of nuclei isolated from mammalian cells. In this chapter, we focus on the use of recombinant MTase (M) M.CviPI to probe chromatin accessibility in nuclei isolated from mammalian cell lines and animal tissue. This technique, named methylation accessibility protocol for individual templates (MAPit), reports protein-DNA interactions at single-molecule resolution. The single-molecule readout allows identification of chromatin subpopulations and rare epigenetic variants within a cell population. Furthermore, the use of M.CviPI in mammalian systems gives a comprehensive view of both chromatin structure and endogenous DNA methylation in a single assay.

Reconfiguration of nucleosome-depleted regions at distal regulatory elements accompanies DNA methylation of enhancers and insulators in cancer

It is well established that cancer-associated epigenetic repression occurs concomitant with CpG island hypermethylation and loss of nucleosomes at promoters, but the role of nucleosome occupancy and epigenetic reprogramming at distal regulatory elements in cancer is still poorly understood. Here, we evaluate the scope of global epigenetic alterations at enhancers and insulator elements in prostate and breast cancer cells using simultaneous genome-wide mapping of DNA methylation and nucleosome occupancy (NOMe-seq). We find that the genomic location of nucleosome-depleted regions (NDRs) is mostly cell type specific and preferentially found at enhancers in normal cells. In cancer cells, however, we observe a global reconfiguration of NDRs at distal regulatory elements coupled with a substantial reorganization of the cancer methylome. Aberrant acquisition of nucleosomes at enhancer-associated NDRs is associated with hypermethylation and epigenetic silencing marks, and conversely, loss of nucleosomes with demethylation and epigenetic activation. Remarkably, we show that nucleosomes remain strongly organized and phased at many facultative distal regulatory elements, even in the absence of a NDR as an anchor. Finally, we find that key transcription factor (TF) binding sites also show extensive peripheral nucleosome phasing, suggesting the potential for TFs to organize NDRs genome-wide and contribute to deregulation of cancer epigenomes. Together, our findings suggest that “decommissioning” of NDRs and TFs at distal regulatory elements in cancer cells is accompanied by DNA hypermethylation susceptibility of enhancers and insulator elements, which in turn may contribute to an altered genome-wide architecture and epigenetic deregulation in malignancy.

Local depletion of DNA methylation identifies a repressive p53 regulatory region in the NEK2 promoter

Genome-scale mapping suggests that the function of DNA methylation varies with genomic context beyond transcriptional repression. However, the use of DNA-demethylating agents (e.g. 5-aza-2'-deoxycytidine (5aza-dC)) to study epigenetic regulation often focuses on gene activation and ignores repression elicited by 5aza-dC. Here, we show that repression of NEK2, which encodes the never in mitosis A (NIMA)-related kinase, by 5aza-dC is context-specific as NEK2 transcript levels were reduced in HCT116 colon cancer cells but not in isogenic p53(-/-) cells. Bisulfite sequencing showed that DNA methylation was restricted to the distal region of the NEK2 promoter. Demethylation by 5aza-dC was associated with increased accessibility to micrococcal nuclease, i.e. nucleosome depletion. Conversely, methyltransferase accessibility protocol for individual templates (MAPit) methylation footprinting showed that nucleosome occupancy and DNA methylation at the distal promoter were significantly increased in p53(-/-) cells, suggesting dynamic regulation of chromatin structure at this region by p53 in HCT116 cells. Stabilization of endogenous p53 by doxorubicin or ectopic expression of p53, but not a p53 DNA-binding mutant, decreased NEK2 expression. Chromatin immunoprecipitation demonstrated direct and specific association of p53 with the distal NEK2 promoter, which was enhanced by doxorubicin. Luciferase reporters confirmed that this region is required for p53-mediated repression of NEK2 promoter activity. Lastly, modulation of p53 abundance altered nucleosome occupancy and DNA methylation at its binding region. These results identify NEK2 as a novel p53-repressed gene, illustrate that its repression by 5aza-dC is specific and associated with nucleosome reorganization, and provide evidence that identification of partially methylated regions can reveal novel p53 target genes.

Multiplex mapping of chromatin accessibility and DNA methylation within targeted single molecules identifies epigenetic heterogeneity in neural stem cells and glioblastoma

Human tumors are comprised of heterogeneous cell populations that display diverse molecular and phenotypic features. To examine the extent to which epigenetic differences contribute to intratumoral cellular heterogeneity, we have developed a high-throughput method, termed MAPit-patch. The method uses multiplexed amplification of targeted sequences from submicrogram quantities of genomic DNA followed by next generation bisulfite sequencing. This provides highly scalable and simultaneous mapping of chromatin accessibility and DNA methylation on single molecules at high resolution. Long sequencing reads from targeted regions maintain the structural integrity of epigenetic information and provide substantial depth of coverage, detecting for the first time minority subpopulations of epigenetic configurations formerly obscured by existing genome-wide and population-ensemble methodologies. Analyzing a cohort of 71 promoters of genes with exons commonly mutated in cancer, MAPit-patch uncovered several differentially accessible and methylated promoters that are associated with altered gene expression between neural stem cell (NSC) and glioblastoma (GBM) cell populations. In addition, considering each promoter individually, substantial epigenetic heterogeneity was observed across the sequenced molecules, indicating the presence of epigenetically distinct cellular subpopulations. At the divergent MLH1/EPM2AIP1 promoter, a locus with three well-defined, nucleosome-depleted regions (NDRs), a fraction of promoter copies with inaccessible chromatin was detected and enriched upon selection of temozolomide-tolerant GBM cells. These results illustrate the biological relevance of epigenetically distinct subpopulations that in part underlie the phenotypic heterogeneity of tumor cell populations. Furthermore, these findings show that alterations in chromatin accessibility without accompanying changes in DNA methylation may constitute a novel class of epigenetic biomarker.

Epigenetic diversity of Kaposi's sarcoma-associated herpesvirus

Spontaneous lytic reactivation of Kaposi’s sarcoma–associated herpesvirus (KSHV) occurs at a low rate in latently infected cells in disease and culture. This suggests imperfect epigenetic maintenance of viral transcription programs, perhaps due to variability in chromatin structure at specific loci across the population of KSHV episomal genomes. To characterize this locus-specific chromatin structural diversity, we used MAPit single-molecule footprinting, which simultaneously maps endogenous CG methylation and accessibility to M.CviPI at GC sites. Diverse chromatin structures were detected at the LANA, RTA and vIL6 promoters. At each locus, chromatin ranged from fully closed to fully open across the population. This diversity has not previously been reported in a virus. Phorbol ester and RTA transgene induction were used to identify chromatin conformations associated with reactivation of lytic transcription, which only a fraction of episomes had. Moreover, certain chromatin conformations correlated with CG methylation patterns at the RTA and vIL6 promoters. This indicated that some of the diverse chromatin conformations at these loci were epigenetically distinct. Finally, by comparing chromatin structures from a cell line infected with constitutively latent virus, we identified products of lytic replication. Our findings show that epigenetic drift can restrict viral propagation by chromatin compaction at latent and lytic promoters.

Genome-wide mapping of nucleosome positioning and DNA methylation within individual DNA molecules

DNA methylation and nucleosome positioning work together to generate chromatin structures that regulate gene expression. Nucleosomes are typically mapped using nuclease digestion requiring significant amounts of material and varying enzyme concentrations. We have developed a method (NOMe-seq) that uses a GpC methyltransferase (M.CviPI) and next generation sequencing to generate a high resolution footprint of nucleosome positioning genome-wide using less than 1 million cells while retaining endogenous DNA methylation information from the same DNA strand. Using a novel bioinformatics pipeline, we show a striking anti-correlation between nucleosome occupancy and DNA methylation at CTCF regions that is not present at promoters. We further show that the extent of nucleosome depletion at promoters is directly correlated to expression level and can accommodate multiple nucleosomes and provide genome-wide evidence that expressed non-CpG island promoters are nucleosome-depleted. Importantly, NOMe-seq obtains DNA methylation and nucleosome positioning information from the same DNA molecule, giving the first genome-wide DNA methylation and nucleosome positioning correlation at the single molecule, and thus, single cell level, that can be used to monitor disease progression and response to therapy.

Simultaneous single-molecule detection of endogenous C-5 DNA methylation and chromatin accessibility using MAPit

Bisulfite genomic sequencing provides a single-molecule view of cytosine methylation states. After deamination, each cloned molecule contains a record of methylation within its sequence. The full power of this technique is harnessed by treating nuclei with an exogenous DNMT prior to DNA extraction. This exogenous methylation marks regions of accessibility and footprints nucleosomes, as well as other DNA-binding proteins. Thus, each cloned molecule records not only the endogenous methylation present (at CG sites, in mammals), but also the exogenous (GC, when using the Chlorella virus protein M.CviPI). We term this technique MAPit, methylation accessibility protocol for individual templates.

A novel lineage of myoviruses infecting cyanobacteria is widespread in the oceans

Viruses infecting bacteria (phages) are thought to greatly impact microbial population dynamics as well as the genome diversity and evolution of their hosts. Here we report on the discovery of a novel lineage of tailed dsDNA phages belonging to the family Myoviridae and describe its first representative, S-TIM5, that infects the ubiquitous marine cyanobacterium, Synechococcus. The genome of this phage encodes an entirely unique set of structural proteins not found in any currently known phage, indicating that it uses lineage-specific genes for virion morphogenesis and represents a previously unknown lineage of myoviruses. Furthermore, among its distinctive collection of replication and DNA metabolism genes, it carries a mitochondrial-like DNA polymerase gene, providing strong evidence for the bacteriophage origin of the mitochondrial DNA polymerase. S-TIM5 also encodes an array of bacterial-like metabolism genes commonly found in phages infecting cyanobacteria including photosynthesis, carbon metabolism and phosphorus acquisition genes. This suggests a common gene pool and gene swapping of cyanophage-specific genes among different phage lineages despite distinct sets of structural and replication genes. All cytosines following purine nucleotides are methylated in the S-TIM5 genome, constituting a unique methylation pattern that likely protects the genome from nuclease degradation. This phage is abundant in the Red Sea and S-TIM5 gene homologs are widespread in the oceans. This unusual phage type is thus likely to be an important player in the oceans, impacting the population dynamics and evolution of their primary producing cyanobacterial hosts.

DNA methyltransferase accessibility protocol for individual templates by deep sequencing

A single-molecule probe of chromatin structure can uncover dynamic chromatin states and rare epigenetic variants of biological importance that bulk measures of chromatin structure miss. In bisulfite genomic sequencing, each sequenced clone records the methylation status of multiple sites on an individual molecule of DNA. An exogenous DNA methyltransferase can thus be used to image nucleosomes and other protein-DNA complexes. In this chapter, we describe the adaptation of this technique, termed Methylation Accessibility Protocol for individual templates, to modern high-throughput sequencing, which both simplifies the workflow and extends its utility.

DNA methylation directly silences genes with non-CpG island promoters and establishes a nucleosome occupied promoter

Despite the fact that 45% of all human gene promoters do not contain CpG islands, the role of DNA methylation in control of non-CpG island promoters is controversial and its relevance in normal and pathological processes is poorly understood. Among the few studies which investigate the correlation between DNA methylation and expression of genes with non-CpG island promoters, the majority do not support the view that DNA methylation directly leads to transcription silencing of these genes. Our reporter assays and gene reactivation by 5-aza-2'-deoxycytidine, a DNA demethylating agent, show that DNA methylation occurring at CpG poor LAMB3 promoter and RUNX3 promoter 1(RUNX3 P1) can directly lead to transcriptional silencing in cells competent to express these genes in vitro. Using Nucleosome Occupancy Methylome- Sequencing, NOMe-Seq, a single-molecule, high-resolution nucleosome positioning assay, we demonstrate that active, but not inactive, non-CpG island promoters display a nucleosome-depleted region (NDR) immediately upstream of the transcription start site (TSS). Furthermore, using NOMe-Seq and clonal analysis, we show that in RUNX3 expressing 623 melanoma cells, RUNX3 P1 has two distinct chromatin configurations: one is unmethylated with an NDR upstream of the TSS; another is methylated and nucleosome occupied, indicating that RUNX3 P1 is monoallelically methylated. Together, these results demonstrate that the epigenetic signatures comprising DNA methylation, histone marks and nucleosome occupancy of non-CpG island promoters are almost identical to CpG island promoters, suggesting that aberrant methylation patterns of non-CpG island promoters may also contribute to tumorigenesis and should therefore be included in analyses of cancer epigenetics.

WIF1 is a frequent target for epigenetic silencing in squamous cell carcinoma of the cervix

Aberrant activation of the Wnt/β-catenin signaling axis is a prominent oncogenic mechanism in numerous cancers including cervical cancer. Wnt inhibitory factor-1 (WIF1) is a secreted protein that binds Wnt and antagonizes Wnt activity. While the WIF1 gene is characterized as a target for epigenetic silencing in some tumor types, WIF1 expression has not been examined in human cervical tissue and cervical cancer. Here, we show that WIF1 is unmethylated and its gene product is expressed in normal cervical epithelium and some cultured cervical tumor lines. In contrast, several cervical cancer lines contained dense CpG methylation within the WIF1 gene, and expression of both WIF1 transcript and protein was restored by culturing cells in the presence of the global DNA demethylating agent 5-aza-2'-deoxycytidine. Using single-molecule MAPit methylation footprinting, we observed differences in chromatin structure within the WIF1 promoter region between cell lines that express and those that do not express WIF1, consistent with transcriptional activity and repression, respectively. The WIF1 promoter was aberrantly methylated in ∼60% (10 of 17) high-grade highly undifferentiated squamous cell cervical tumors examined, whereas paired normal tissue showed significantly lower levels of CpG methylation. WIF1 protein was not detectable by immunohistochemistry in tumors with quantitatively high levels of WIF1 methylation. Of note, WIF1 protein was not detectable in two of the seven unmethylated cervical tumors examined, suggesting other mechanisms may contribute WIF1 repression. Our findings establish the WIF1 gene as a frequent target for epigenetic silencing in squamous cell carcinoma of the cervix.

Simultaneous single-molecule mapping of protein-DNA interactions and DNA methylation by MAPit

Sites of protein binding to DNA are inferred from footprints or spans of protection against a probing reagent. In most protocols, sites of accessibility to a probe are detected by mapping breaks in DNA strands. As discussed in this unit, such methods obscure molecular heterogeneity by averaging cuts at a given site over all DNA strands in a sample population. The DNA methyltransferase accessibility protocol for individual templates (MAPit), an alternative method described in this unit, localizes protein-DNA interactions by probing with cytosine-modifying DNA methyltransferases followed by bisulfite sequencing. Sequencing individual DNA products after amplification of bisulfite-converted sequences permits assignment of the methylation status of every enzyme target site along a single DNA strand. Use of the GC-methylating enzyme M.CviPI allows simultaneous mapping of chromatin accessibility and endogenous CpG methylation. MAPit is therefore the only footprinting method that can detect subpopulations of molecules with distinct patterns of protein binding or chromatin architecture and correlate them directly with the occurrence of endogenous methylation. Additional advantages of MAPit methylation footprinting as well as considerations for experimental design and potential sources of error are discussed.

Identification of restriction endonucleases sensitive to 5-cytosine methylation at non-CpG sites, including expanded (CAG)n/(CTG)n repeats

Most epigenetic studies assess methylation of 5'-CpG-3' sites but recent evidence indicates that non-CpG cytosine methylation occurs at high levels in humans and other species. This is most prevalent at 5'-CHG-3', where H = A, C or T, and it preferentially occurs at 5'-CpA-3' and 5'-CpT-3' sites. With the goal of facilitating the detection of non-CpG methylation, the restriction endonucleases ApeKI, BbvI, EcoP15I, Fnu 4HI, MwoI and TseI were assessed for their sensitivity to 5-methylcytosine at GpCpA, GpCpT, GpCpC or GpCpG sites, where methylation is catalyzed by the DNA 5-cytosine 5'-GpC-3' methyltransferase M.CviPI. We tested a variety of sequences including various plasmid-based sites, a cloned disease-associated (CAG)83•(CTG)83 repeat and in vitro synthesized tracts of only (CAG)500•(CTG)500 or (CAG)800•(CTG)800. The repeat tracts are enriched for the preferred CpA and CpT motifs. We found that none of the tested enzymes can cleave their recognition sequences when they are 5'-GpC-3' methylated. A genomic site known to convert its non-CpG methylation levels upon C2C12 differentiation was confirmed through the use of these enzymes. These enzymes can be useful in rapidly and easily determining the most common non-CpG methylation status in various sequence contexts, as well as at expansions of (CAG)n•(CTG)n repeat tracts associated with diseases like myotonic dystrophy and Huntington disease.

MethylViewer: computational analysis and editing for bisulfite sequencing and methyltransferase accessibility protocol for individual templates (MAPit) projects

Bisulfite sequencing is a widely-used technique for examining cytosine DNA methylation at nucleotide resolution along single DNA strands. Probing with cytosine DNA methyltransferases followed by bisulfite sequencing (MAPit) is an effective technique for mapping protein-DNA interactions. Here, MAPit methylation footprinting with M.CviPI, a GC methyltransferase we previously cloned and characterized, was used to probe hMLH1 chromatin in HCT116 and RKO colorectal cancer cells. Because M.CviPI-probed samples contain both CG and GC methylation, we developed a versatile, visually-intuitive program, called MethylViewer, for evaluating the bisulfite sequencing results. Uniquely, MethylViewer can simultaneously query cytosine methylation status in bisulfite-converted sequences at as many as four different user-defined motifs, e.g. CG, GC, etc., including motifs with degenerate bases. Data can also be exported for statistical analysis and as publication-quality images. Analysis of hMLH1 MAPit data with MethylViewer showed that endogenous CG methylation and accessible GC sites were both mapped on single molecules at high resolution. Disruption of positioned nucleosomes on single molecules of the PHO5 promoter was detected in budding yeast using M.CviPII, increasing the number of enzymes available for probing protein-DNA interactions. MethylViewer provides an integrated solution for primer design and rapid, accurate and detailed analysis of bisulfite sequencing or MAPit datasets from virtually any biological or biochemical system.

Hypomethylation of a LINE-1 promoter activates an alternate transcript of the MET oncogene in bladders with cancer

It was recently shown that a large portion of the human transcriptome can originate from within repetitive elements, leading to ectopic expression of protein-coding genes. However the mechanism of transcriptional activation of repetitive elements has not been definitively elucidated. For the first time, we directly demonstrate that hypomethylation of retrotransposons can cause altered gene expression in humans. We also reveal that active LINE-1s switch from a tetranucleosome to dinucleosome structure, acquiring H2A.Z- and nucleosome-free regions upstream of TSSs, previously shown only at active single-copy genes. Hypomethylation of a specific LINE-1 promoter was also found to induce an alternate transcript of the MET oncogene in bladder tumors and across the entire urothelium of tumor-bearing bladders. These data show that, in addition to contributing to chromosomal instability, hypomethylation of LINE-1s can alter the functional transcriptome and plays a role not only in human disease but also in disease predisposition.

A surprisingly large portion of our transcriptome originates within repetitive elements, most commonly LINE-1s. However, the mechanism of activation has not been definitively shown. We directly demonstrate for the first time the causal relationship between DNA hypomethylation and transcriptional activation of LINE-1 promoters. Hypomethylation of specific LINE-1 promoters can alter the transcriptome, including activating an alternate transcript of the MET oncogene, not only in primary bladder tumors but also in premalignant urothelium across entire bladders with tumors. Our study has important implications for tumor biology, cancer detection, and treatment, and it also answers the long-standing question of whether hypomethylation of retrotransposons induces ectopic gene expression and influences disease susceptibility in humans, a phenomenon first described in agouti mice.

Preparation of genome-wide DNA fragment libraries using bisulfite in polyacrylamide gel electrophoresis slices with formamide denaturation and quality control for massively parallel sequencing by oligonucleotide ligation and detection

Bisulfite sequencing is widely used for analysis of DNA methylation status (i.e., 5-methylcytosine [5mC] vs. cytosine [C]) in CpG-rich or other loci in genomic DNA (gDNA). Such methods typically involve reaction of gDNA with bisulfite followed by polymerase chain reaction (PCR) amplification of specific regions of interest that, overall, converts C-->T (thymine) and 5mC-->C and then capillary sequencing to measure C versus T composition at CpG sites. Massively parallel sequencing by oligonucleotide ligation and detection (SOLiD) has recently enabled relatively low-cost whole genome sequencing, and it would be highly desirable to apply such massively parallel sequencing to bisulfite-converted whole genomes to determine DNA methylation status of an entire genome, which has heretofore not been reported. As an initial step toward achieving this goal, we have extended our ongoing interest in improving bisulfite conversion sample preparation to include a human genome-wide fragment library for SOliD. The current article features novel use of formamide denaturant during bisulfite conversion of a suitably constructed library directly in a band slice from polyacryamide gel electrophoresis (PAGE). To validate this new protocol for 5mC-protected fragment library conversion, which we refer to as Bis-PAGE, capillary-based size analysis and Sanger sequencing were carried out for individual amplicons derived from single-molecule PCR (smPCR) of randomly selected library fragments. smPCR/Capillary Sanger sequencing of approximately 200 amplicons unambiguously demonstrated greater than 99% C-->T conversion. All of these approximately 200 Sanger sequences were analyzed with a previously published web-accessible bioinformatics tool (methBLAST) for mapping to human chromosomes, the results of which indicated random distribution of analyzed fragments across all chromosomes. Although these particular Bis-PAGE conversion and quality control methods were exemplified in the context of a fragment library for SOLiD, the concepts can be generalized to include other genome-wide library constructions intended for DNA methylation analysis by alternative high-throughput or massively parallelized methods that are currently available.

Single-molecule analysis of chromatin: changing the view of genomes one molecule at a time

Wrapping DNA into chromatin provides a wealth of regulatory mechanisms that ensure normal growth and development in eukaryotes. Our understanding of chromatin structure, including nucleosomes and non-histone protein-DNA interactions, has benefited immensely from nuclease and chemical digestion techniques. DNA-bound proteins, such as histones or site-specific factors, protect DNA against nuclease cleavage and generate large nucleosomal or small regulatory factor footprints. Chromatin subject to distinct modes of regulation often coincides with sites of nuclease hypersensitivity or nucleosome positioning. An inherent limitation of cleavage-based analyses has been the inability to reliably analyze regions of interest when levels of digestion depart from single-hit kinetics. Moreover, cleavage-based techniques provide views that are averaged over all the molecules in a sample population. Therefore, in cases of occupancy of multiple regulatory elements by factors, one cannot define whether the factors are bound to the same or different molecules in the population. The recent development of DNA methyltransferase-based, single-molecule MAP-IT technology overcomes limitations of ensemble approaches and has opened numerous new avenues in chromatin research. Here, we review the strengths, limitations, applications and future prospects of MAP-IT ranging from structural issues to mechanistic questions in eukaryotic chromatin regulation.

DNA methyltransferase probing of chromatin structure within populations and on single molecules

Non-invasive methods for mapping chromatin structure are necessary for creating an accurate view of genome function and dynamics in vivo. Ectopic induction of cytosine-5 DNA methyltransferases (C5 MTases) in Saccharomyces cerevisiae is a powerful technique for probing chromatin structure with minimal disruption to yeast physiology. Accessibility of MTases to their cognate sites is impaired based on the strength and span of the protein-DNA interaction to be probed. Methylated cytosines that resist chemical deamination are detected positively by the PCR-based technique of bisulfite genomic sequencing. PCR amplicons can be sequenced directly yielding an average m(5)C frequency or accessibility of each target site within the population, a technique termed methyltransferase accessibility protocol (MAP). More recently, the sequencing of cloned molecules in MAP for individual templates (MAPit) enables assignment of the methylation status of each target site along a continuous DNA strand from a single cell. The unique capability to score methylation at multiple sites in single molecules permits detection of inherent structural variability in chromatin. Here, MAPit analysis of the repressed and induced PHO5 promoter of budding yeast, using a C5 MTase with dinucleotide recognition specificity, reveals considerable cell-to-cell heterogeneity in chromatin structure. Substantial variation is observed in the extent to which the MTase gains entry to each of the nucleosomes positioned at PHO5, suggesting differences in their intrinsic thermodynamic stability in vivo. MAPit should be readily adaptable to the analysis of chromatin structure and non-histone protein-DNA interactions in a variety of model systems.

Single-molecule and population probing of chromatin structure using DNA methyltransferases

Probing chromatin structure with DNA methyltransferases offers advantages over more commonly used nuclease-based and chromatin immunoprecipitation methods for detection of nucleosomes and non-histone protein-DNA interactions. Here, we describe two related methods in which the readout of MTase accessibility is obtained by assaying 5-methylcytosine in DNA through the PCR-based technique of bisulfite genomic sequencing. The methyltransferase accessibility protocol (MAP) determines the relative frequency at which the enzyme accesses each of its target sites over an entire population of PCR amplified product. While MAP yields much quantitative information about relative accessibility of a region of chromatin, a complementary single-molecule view of methyltransferase accessibility, termed MAP for individual templates (MAP-IT), is provided by analysis of cloned PCR products. Absolute rather than relative methylation frequencies in a region are obtained by summing the methylation status at each site over a cohort of clones. Moreover, as the integrity of individual molecules is maintained in MAP-IT, unique information about the distribution of multiple footprints along continuous regions is gleaned. In principle, the population MAP and single-molecule MAP-IT strategies can be used to analyze chromatin structure in a variety of model systems. Here, we describe the application of MAP in living Saccharomyces cerevisiae cells and MAP-IT in the analysis of a mammalian tumor suppressor gene in nuclei. This application of MAP-IT provides the first means to simultaneously determine CpG methylation of mammalian genes and their overlying chromatin structure in the same single DNA molecule.

Cloning of Nt.CviQII nicking endonuclease and its cognate methyltransferase: M.CviQII methylates AG sequences

Chlorella virus NY-2A has a large, highly methylated dsDNA genome (45% of the cytosines are 5-methylcytosine and 37% of the adenines are N(6)-methyladenine). Here, we report the cloning, expression, and characterization of the NY-2A-encoded CviQII nicking-modification (N-M) system. The nicking endonuclease, Nt.CviQII, recognizes R downward arrowAG (R=A or G, downward arrow indicating cleavage site) sequences and cleaves the phosphodiester bond 5' to the adenosine. Because of the difficulty in cloning and expressing the wild-type Nt.CviQII, C-terminal truncation mutants were generated and full-length Nt.CviQII was reconstructed by intein-mediated peptide ligation. The truncation mutants and the reconstructed full-length Nt.CviQII have the same recognition and cleavage specificity as the native enzyme. Full-length and truncated Nt.CviQII produced by a cell-free transcription/translation system have similar reaction rates. The methyltransferase, M.CviQII, was also cloned and expressed. It modifies the adenine in AG doublets of DNA in vitro and in vivo in Escherichia coli. To our knowledge, M.CviQII is the first adenine methyltransferase that recognizes a dinucleotide. Therefore, M.CviQII may be a useful reagent for blocking endonuclease cleavage when restriction sites overlap with AG sequences.

Active PHO5 chromatin encompasses variable numbers of nucleosomes at individual promoters

Transcriptional activation is often associated with chromatin remodeling. However, little is known about the dynamics of remodeling of nucleosome arrays in vivo. Upon induction of Saccharomyces cerevisiae PHO5, a novel kinetic assay of DNA methyltransferase accessibility showed that nucleosomes adjacent to the histone-free upstream activating sequence (UASp1) are disrupted earlier and at higher frequency in the cell population than are those more distal. Individually cloned molecules, each representing the chromatin state of a full promoter from a single cell, revealed multiple promoter classes with either no remodeling or variable numbers of disrupted nucleosomes. Individual promoters in the remodeled fraction were highly enriched for contiguous blocks of disrupted nucleosomes, the majority of which overlapped the UAS region. These results support a probabilistic model in which chromatin remodeling at PHO5 spreads from sites of transactivator association with DNA and attenuates with distance.

Mapping chromatin structure in vivo using DNA methyltransferases

Cytosine-5 DNA methyltransferases (C5 DMTases) are effective reagents for analyzing chromatin and footprinting DNA-bound factors in vivo. Cytosine methylation in accessible regions is assayed positively by the PCR-based technique of bisulfite sequencing. In this article, we outline two complementary uses for the DNA methyltransferase CviPI (M.CviPI, GC specificity) in probing chromatin organization. First, we describe the use of the naturally occurring, free enzyme as a diffusible probe to map changes in nucleosome structure and to footprint factor interactions at cis-regulatory sequences. In a second application, termed targeted gene methylation (TAGM), the DMTase is targeted via in-frame fusion to a DNA-binding factor. The rapid accumulation of DNA methylation enables highly sensitive detection of factor binding. Both strategies can be applied with any C5 DMTase, such as M.SssI, which also possesses a short-recognition specificity (CG). A description of methods for constructing C5 DMTase-expressing strains of Saccharomyces cerevisiae and analyzing chromatin regions is provided. We also include comprehensive protocols for the isolation and bisulfite treatment of genomic DNA as well as the subsequent bisulfite sequencing steps. Data demonstrating the efficacy of both DMTase probing techniques, theoretical considerations, and experimental analyses are presented at GAL1 and PHO5.

The use of prokaryotic DNA methyltransferases as experimental and analytical tools in modern biology

Prokaryotic DNA methyltransferases (MTases) are used as experimental and research tools in molecular biology and molecular genetics due to their ability to recognize and transfer methyl groups to target bases in specific DNA sequences. As a practical tool, prokaryotic DNA MTases can be used in recombinant DNA technology for in vitro alteration and enhancing of cleavage specificity of restriction endonucleases. The ability of prokaryotic DNA MTases to methylate cytosine residues in specific sequences, which are also methylated in eukaryotic DNA, makes it possible to use them as analytical reagent for determination of the site-specific level of methylation in eukaryotic DNA. In vivo DNA methylation by prokaryotic DNA MTases is used in different techniques for probing chromatin structure and protein-DNA interactions. Additional prospects are opened by development of the methods of DNA methylation targeted to predetermined DNA sequences by fusion of DNA MTases to DNA binding proteins. This review will discuss the application of prokaryotic DNA MTases of Type II in the methods and approaches mentioned above.

Cloning of CviPII nicking and modification system from chlorella virus NYs-1 and application of Nt.CviPII in random DNA amplification

The cloning and expression of the CviPII DNA nicking and modification system encoded by chlorella virus NYs-1 is described. The system consists of a co-linear MTase encoding gene (cviPIIM) and a nicking endonuclease encoding gene (cviPIINt) separated by 12 nt. M.CviPII possesses eight conserved amino acid motifs (I to VIII) typical of C5 MTases, but, like another chlorella virus MTase M.CviJI, lacks conserved motifs IX and X. In addition to modification of the first cytosine in CCD (D = A, G or T) sequences, M.CviPII modifies both the first two cytosines in CCAA and CCCG sites as well. Nt.CviPII has significant amino acid sequence similarity to Type II restriction endonuclease CviJI that recognizes an overlapping sequence (RG--CY). Nt.CviPII was expressed in Escherichia coli with or without a His-tag in a host pre-modified by M.CviPII. Recombinant Nt.CviPII recognizes the DNA sequence CCD and cleaves the phosphodiester bond 5' of the first cytosine while the other strand of DNA at this site is not affected. Nt.CviPII displays site preferences with CCR (R = A or G) sites preferred over CCT sites. Nt.CviPII is active from 16 to 65 degrees C with a temperature optimum of 30-45 degrees C. Nt.CviPII can be used to generate single-stranded DNAs (ssDNAs) for isothermal strand-displacement amplification. Nt.CviPII was used in combination with Bst DNA polymerase I large fragment to rapidly amplify anonymous DNA from genomic DNA or from a single bacterial colony.