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

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

Long-range phasing of dynamic, tissue-specific and allele-specific regulatory elements

Epigenomic maps identify gene regulatory elements by their chromatin state. However, prevailing short-read sequencing methods cannot effectively distinguish alleles, evaluate the interdependence of elements in a locus or capture single-molecule dynamics. Here, we apply targeted nanopore sequencing to profile chromatin accessibility and DNA methylation on contiguous ~100-kb DNA molecules that span loci relevant to development, immunity and imprinting. We detect promoters, enhancers, insulators and transcription factor footprints on single molecules based on exogenous GpC methylation. We infer relationships among dynamic elements within immune loci, and order successive remodeling events during T cell stimulation. Finally, we phase primary sequence and regulatory elements across the H19/IGF2 locus, uncovering primate-specific features. These include a segmental duplication that stabilizes the imprinting control region and a noncanonical enhancer that drives biallelic IGF2 expression in specific contexts. Our study advances emerging strategies for phasing gene regulatory landscapes and reveals a mechanism that overrides IGF2 imprinting in human cells.

Generating specificity in genome regulation through transcription factor sensitivity to chromatin

Cell type-specific gene expression relies on transcription factors (TFs) binding DNA sequence motifs embedded in chromatin. Understanding how motifs are accessed in chromatin is crucial to comprehend differential transcriptional responses and the phenotypic impact of sequence variation. Chromatin obstacles to TF binding range from DNA methylation to restriction of DNA access by nucleosomes depending on their position, composition and modification. In vivo and in vitro approaches now enable the study of TF binding in chromatin at unprecedented resolution. Emerging insights suggest that TFs vary in their ability to navigate chromatin states. However, it remains challenging to link binding and transcriptional outcomes to molecular characteristics of TFs or the local chromatin substrate. Here, we discuss our current understanding of how TFs access DNA in chromatin and novel techniques and directions towards a better understanding of this critical step in genome regulation.

BANP opens chromatin and activates CpG-island-regulated genes

The majority of gene transcripts generated by RNA polymerase II in mammalian genomes initiate at CpG island (CGI) promoters^1,2, yet our understanding of their regulation remains limited. This is in part due to the incomplete information that we have on transcription factors, their DNA-binding motifs and which genomic binding sites are functional in any given cell type^3-5. In addition, there are orphan motifs without known binders, such as the CGCG element, which is associated with highly expressed genes across human tissues and enriched near the transcription start site of a subset of CGI promoters^6-8. Here we combine single-molecule footprinting with interaction proteomics to identify BTG3-associated nuclear protein (BANP) as the transcription factor that binds this element in the mouse and human genome. We show that BANP is a strong CGI activator that controls essential metabolic genes in pluripotent stem and terminally differentiated neuronal cells. BANP binding is repelled by DNA methylation of its motif in vitro and in vivo, which epigenetically restricts most binding to CGIs and accounts for differential binding at aberrantly methylated CGI promoters in cancer cells. Upon binding to an unmethylated motif, BANP opens chromatin and phases nucleosomes. These findings establish BANP as a critical activator of a set of essential genes and suggest a model in which the activity of CGI promoters relies on methylation-sensitive transcription factors that are capable of chromatin opening.

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.

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.

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.

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.

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.

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.