Mapping chromatin structure in vivo using DNA methyltransferases

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

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.

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.

Cooperative binding between distant transcription factors is a hallmark of active enhancers

Enhancers harbor binding motifs that recruit transcription factors (TFs) for gene activation. While cooperative binding of TFs at enhancers is known to be critical for transcriptional activation of a handful of developmental enhancers, the extent of TF cooperativity genome-wide is unknown. Here, we couple high-resolution nuclease footprinting with single-molecule methylation profiling to characterize TF cooperativity at active enhancers in the Drosophila genome. Enrichment of short micrococcal nuclease (MNase)-protected DNA segments indicates that the majority of enhancers harbor two or more TF-binding sites, and we uncover protected fragments that correspond to co-bound sites in thousands of enhancers. From the analysis of co-binding, we find that cooperativity dominates TF binding in vivo at the majority of active enhancers. Cooperativity is highest between sites spaced 50 bp apart, indicating that cooperativity occurs without apparent protein-protein interactions. Our findings suggest nucleosomes promoting cooperativity because co-binding may effectively clear nucleosomes and promote enhancer function.

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.

Enhancement of LacI binding in vivo

Transcription factors (TFs) bind to specific sequences in DNA to regulate transcription. Despite extensive measurements of TFs’ dissociation constant (K_d) in vitro, their apparent K_din vivo are usually unknown. LacI, a bacterial TF, is often used to artificially recruit proteins onto eukaryotic genomes. As LacI binds tightly to its recognition site (LacO) in vitro with a K_d about 10 picomolar (pM), it is often assumed that LacI also has high affinity to LacO in vivo. In this work, we measured LacI binding in living yeast cells using a fluorescent repressor operator system and found an apparent K_d of ∼0.6 μM, four orders of magnitude higher than that in vitro. By genetically altering (i) GFP-LacI structure, (ii) GFP-LacI stability, (iii) chromosome accessibility and (iv) LacO sequence, we reduced the apparent K_d to <10 nM. It turns out that the GFP tagging location and the fusion protein stability have a large effect on LacI binding, but surprisingly, chromosome accessibility only plays a mild role. These findings contribute to our quantitative understanding of the features that affect the apparent K_d of TF in cells. They also provide guidance for future design of more specific chromosomal recruitment through high-affinity TFs.

Absolute nucleosome occupancy map for the Saccharomyces cerevisiae genome

Mapping of nucleosomes, the basic DNA packaging unit in eukaryotes, is fundamental for understanding genome regulation because nucleosomes modulate DNA access by their positioning along the genome. A cell-population nucleosome map requires two observables: nucleosome positions along the DNA ("Where?") and nucleosome occupancies across the population ("In how many cells?"). All available genome-wide nucleosome mapping techniques are yield methods because they score either nucleosomal (e.g., MNase-seq, chemical cleavage-seq) or nonnucleosomal (e.g., ATAC-seq) DNA but lose track of the total DNA population for each genomic region. Therefore, they only provide nucleosome positions and maybe compare relative occupancies between positions, but cannot measure absolute nucleosome occupancy, which is the fraction of all DNA molecules occupied at a given position and time by a nucleosome. Here, we established two orthogonal and thereby cross-validating approaches to measure absolute nucleosome occupancy across the Saccharomyces cerevisiae genome via restriction enzymes and DNA methyltransferases. The resulting high-resolution (9-bp) map shows uniform absolute occupancies. Most nucleosome positions are occupied in most cells: 97% of all nucleosomes called by chemical cleavage-seq have a mean absolute occupancy of 90 ± 6% (±SD). Depending on nucleosome position calling procedures, there are 57,000 to 60,000 nucleosomes per yeast cell. The few low absolute occupancy nucleosomes do not correlate with highly transcribed gene bodies, but correlate with increased presence of the nucleosome-evicting chromatin structure remodeling (RSC) complex, and are enriched upstream of highly transcribed or regulated genes. Our work provides a quantitative method and reference frame in absolute terms for future chromatin studies.

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.

High-resolution digital profiling of the epigenome

The widespread adoption of short-read DNA sequencing as a digital epigenomic readout platform has motivated the development of genome-wide tools that achieve base-pair resolution. New methods for footprinting and affinity purification of nucleosomes, RNA polymerases, chromatin remodellers and transcription factors have increased the resolution of epigenomic profiling by two orders of magnitude, leading to new insights into how the chromatin landscape affects gene regulation. These digital epigenomic tools have also been applied to directly profile both turnover kinetics and transcription in situ. In this Review, we describe how these new genome-wide tools allow interrogation of diverse aspects of the epigenome.

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.

Understanding the connection between epigenetic DNA methylation and nucleosome positioning from computer simulations

Cytosine methylation is one of the most important epigenetic marks that regulate the process of gene expression. Here, we have examined the effect of epigenetic DNA methylation on nucleosomal stability using molecular dynamics simulations and elastic deformation models. We found that methylation of CpG steps destabilizes nucleosomes, especially when these are placed in sites where the DNA minor groove faces the histone core. The larger stiffness of methylated CpG steps is a crucial factor behind the decrease in nucleosome stability. Methylation changes the positioning and phasing of the nucleosomal DNA, altering the accessibility of DNA to regulatory proteins, and accordingly gene functionality. Our theoretical calculations highlight a simple physical-based explanation on the foundations of epigenetic signaling.

DNA methylation and nucleosome occupancy regulate the cancer germline antigen gene MAGEA11

MAGEA11 is a cancer germline (CG) antigen and androgen receptor co-activator. Its expression in cancers other than prostate, and its mechanism of activation, has not been reported. In silico analyses reveal that MAGEA11 is frequently expressed in human cancers, is increased during tumor progression, and correlates with poor prognosis and survival. In prostate and epithelial ovarian cancers (EOC), MAGEA11 expression was associated with promoter and global DNA hypomethylation, and with activation of other CG genes. Pharmacological or genetic inhibition of DNA methyltransferases (DNMTs) and/or histone deacetylases (HDACs) activated MAGEA11 in a cell line specific manner. MAGEA11 promoter activity was directly repressed by DNA methylation, and partially depended on Sp1, as pharmacological or genetic targeting of Sp1 reduced MAGEA11 promoter activity and endogenous gene expression. Importantly, DNA methylation regulated nucleosome occupancy specifically at the -1 positioned nucleosome of MAGEA11. Methylation of a single Ets site near the transcriptional start site (TSS) correlated with -1 nucleosome occupancy and, by itself, strongly repressed MAGEA11 promoter activity. Thus, DNA methylation regulates nucleosome occupancy at MAGEA11, and this appears to function cooperatively with sequence-specific transcription factors to regulate gene expression. MAGEA11 regulation is highly instructive for understanding mechanisms regulating CG antigen genes in human cancer.

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.

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.

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.

Distributed probing of chromatin structure in vivo reveals pervasive chromatin accessibility for expressed and non-expressed genes during tissue differentiation in C. elegans

Background

Tissue differentiation is accompanied by genome-wide changes in the underlying chromatin structure and dynamics, or epigenome. By controlling when, where, and what regulatory factors have access to the underlying genomic DNA, the epigenome influences the cell's transcriptome and ultimately its function. Existing genomic methods for analyzing cell-type-specific changes in chromatin generally involve two elements: (i) a source for purified cells (or nuclei) of distinct types, and (ii) a specific treatment that partitions or degrades chromatin by activity or structural features. For many cell types of great interest, such assays are limited by our inability to isolate the relevant cell populations in an organism or complex tissue containing an intertwined mixture of other cells. This limitation has confined available knowledge of chromatin dynamics to a narrow range of biological systems (cell types that can be sorted/separated/dissected in large numbers and tissue culture models) or to amalgamations of diverse cell types (tissue chunks, whole organisms).

Results

Transgene-driven expression of DNA/chromatin modifying enzymes provides one opportunity to query chromatin structures in expression-defined cell subsets. In this work we combine in vivo expression of a bacterial DNA adenine methyltransferase (DAM) with high throughput sequencing to sample tissue-specific chromatin accessibility on a genome-wide scale. We have applied the method (DALEC: Direct Asymmetric Ligation End Capture) towards mapping a cell-type-specific view of genome accessibility as a function of differentiated state. Taking advantage of C. elegans strains expressing the DAM enzyme in diverse tissues (body wall muscle, gut, and hypodermis), our efforts yield a genome-wide dataset measuring chromatin accessibility at each of 538,000 DAM target sites in the C. elegans (diploid) genome.

Conclusions

Validating the DALEC mapping results, we observe a strong association between observed coverage by nucleosomes and low DAM accessibility. Strikingly, we observed no extended regions of inaccessible chromatin for any of the tissues examined. These results are consistent with "local choreography" models in which differential gene expression is driven by intricate local rearrangements of chromatin structure rather than gross impenetrability of large chromosomal regions.

In vivo methylation of mtDNA reveals the dynamics of protein-mtDNA interactions

To characterize the organization of mtDNA–protein complexes (known as nucleoids) in vivo, we have probed the mtDNA surface exposure using site-specific DNA methyltransferases targeted to the mitochondria. We have observed that DNA methyltransferases have different accessibility to different sites on the mtDNA based on the levels of protein occupancy. We focused our studies on selected regions of mtDNA that are believed to be major regulatory regions involved in transcription and replication. The transcription termination region (TERM) within the tRNA^Leu(UUR) gene was consistently and strongly protected from methylation, suggesting frequent and high affinity binding of mitochondrial transcription termination factor 1 (mTERF1) to the site. Protection from methylation was also observed in other regions of the mtDNA, including the light and heavy strand promoters (LSP, HSP) and the origin of replication of the light strand (OL). Manipulations aiming at increasing or decreasing the levels of the mitochondrial transcription factor A (TFAM) led to decreased in vivo methylation, whereas manipulations that stimulated mtDNA replication led to increased methylation. We also analyzed the effect of ATAD3 and oxidative stress in mtDNA exposure. Our data provide a map of human mtDNA accessibility and demonstrate that nucleoids are dynamically associated with proteins.

Genome-wide mapping of nucleosome positions in Schizosaccharomyces pombe

The majority of nuclear eukaryotic DNA is packaged into nucleosome cores where DNA is wrapped tightly around histone protein octamers. Such histone bound nucleosomal DNA is less accessible than the short linker DNA between nucleosome cores or the DNA in extended nucleosome free regions. Therefore, the positions of nucleosomes relative to a DNA sequence feature, like a transactivator binding site, a transcriptional start site or an origin of replication, can have profound effects on nuclear processes like transcription, replication, recombination and repair. Now that many DNA related processes are studied in a genome-wide manner, it is increasingly important to map the basic organization of their chromosomal DNA substrate, i.e., the positions of nucleosomes, on a genome-wide scale as well. To this end, the protection of nucleosomal DNA from digestion with micrococcal nuclease (MNase) is used as an assay for the presence of a nucleosome. The MNase protected DNA fragments, so called mononucleosomal DNA, can be mapped genome-wide by hybridization to microarrays. This method has been established for Saccharomyces cerevisiae, and we present here the adaptation of the method for Schizosaccharomyces pombe. As an independent method to validate genome-wide data for individual loci, we also include a protocol for the determination of locus specific nucleosome positioning by indirect end labeling.

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.

Application of DNA methyltransferases in targeted DNA methylation

DNA methylation is an essential epigenetic modification. In bacteria, it is involved in gene regulation, DNA repair, and control of cell cycle. In eukaryotes, it acts in concert with other epigenetic modifications to regulate gene expression and chromatin structure. In addition to these biological roles, DNA methyltransferases have several interesting applications in biotechnology, which are the main focus of this review, namely, (1) in vivo footprinting: as several bacterial DNA methyltransferases cannot methylate DNA bound to histone proteins, the pattern of DNA methylation after expression of DNA methyltransferases in the cell allows determining nucleosome positioning; (2) mapping the binding specificity of DNA binding proteins: after fusion of a DNA methyltransferase to a DNA-binding protein and expression of the fusion protein in a cell, the DNA methylation pattern reflects the DNA-binding specificity of the DNA-binding protein; and (3) targeted gene silencing: after fusion of a DNA methyltransferase to a suitable DNA-binding domain, DNA methylation can be directed to promoter regions of target genes. Thereby, gene expression can be switched off specifically, efficiently, and stably, which has a number of potential medical applications.

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.

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.

Footprinting of mammalian promoters: use of a CpG DNA methyltransferase revealing nucleosome positions at a single molecule level

Promoters are molecular ‘modules’, which are controlled as individual entities yet are often analyzed by nuclease digestion methodologies which, a priori, destroy this modularity. About 40% of mammalian genes contain CpG islands in their promoters and exonic regions, which are normally unmethylated. We developed a footprinting strategy to map the chromatin structure at unmethylated CpG islands by treatment of isolated nuclei with the CpG-specific DNA methyltransferase SssI (M.SssI), followed by genomic bisulfite sequencing of individual progeny DNA molecules. This gave single molecule resolution over the promoter region and allowed for the physical linkage between binding sites on individual promoter molecules to be maintained. Comparison of the p16 promoters in two human cell lines, J82 and LD419, expressing the p16 gene at 25-fold different levels showed that the two cell lines contain remarkably different, heterogeneously positioned nucleosomes over the promoter region, which were not distinguishable by standard methods using nucleases. Our high resolution approach gives a ‘digitized’ visualization of each promoter providing information regarding nucleosome occupancy and may be utilized to define transcription factor binding and chromatin remodeling.