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.

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.

TRIM29 suppresses TWIST1 and invasive breast cancer behavior

TRIM29 (ATDC) exhibits a contextual function in cancer, but seems to exert a tumor-suppressor role in breast cancer. Here, we show that TRIM29 is often silenced in primary breast tumors and cultured tumor cells as a result of aberrant gene hypermethylation. RNAi-mediated silencing of TRIM29 in breast tumor cells increased their motility, invasiveness, and proliferation in a manner associated with increased expression of mesenchymal markers (N-cadherin and vimentin), decreased expression of epithelial markers (E-cadherin and EpCAM), and increased expression and activity of the oncogenic transcription factor TWIST1, an important driver of the epithelial-mesenchymal transition (EMT). Functional investigations revealed an inverse relationship in the expression of TRIM29 and TWIST1, suggesting the existence of a negative regulatory feedback loop. In support of this relationship, we found that TWIST1 inhibited TRIM29 promoter activity through direct binding to a region containing a cluster of consensus E-box elements, arguing that TWIST1 transcriptionally represses TRIM29 expression. Analysis of a public breast cancer gene-expression database indicated that reduced TRIM29 expression was associated with reduced relapse-free survival, increased tumor size, grade, and metastatic characteristics. Taken together, our results suggest that TRIM29 acts as a tumor suppressor in breast cancer through its ability to inhibit TWIST1 and suppress EMT.

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.

Fibroblast growth factor receptor 2 homodimerization rapidly reduces transcription of the pluripotency gene Nanog without dissociation of activating transcription factors

Nanog or Gata6-positive cells co-exist and are convertible within the inner cell mass of murine blastocysts and embryonic stem (ES) cells. Previous studies demonstrate fibroblast growth factor receptor 2 (FGFR2) triggers Nanog gene down-regulation and differentiation to primitive endoderm (PE); however, the underlying mechanisms responsible for reversible and fluctuating cell fate are poorly understood. Using an inducible FGFR2 dimerization system in ES cells, we demonstrate that FGFR2 activation rapidly down-regulated Nanog gene transcription through activation of the Mek pathway and subsequently differentiated ES cells into PE cells. FGFR2 rather selectively repressed the Nanog gene with minimal effect on other pluripotency genes, including Oct4 and Sox2. We determined the Nanog promoter region containing minimum Oct4/Sox2 binding sites was sufficient for this transcriptional down-regulation by FGFR2, when the reporter transgenes were integrated with insulators. Of interest, FGFR2-mediated Nanog transcriptional reduction occurred without dissociation of RNA polymerase II, p300, Oct4, Sox2, and Tet1 from the Nanog proximal promoter region and with no increase in repressive histone methylation marks or DNA methylation, implying the gene repression is in the early and transient phase. Furthermore, addition of a specific FGFR inhibitor readily reversed this Nanog repression status. These findings illustrate well how FGFR2 induces rapid but reversible Nanog repression within ES cells.

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.

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.

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.

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.

Bisulfite sequencing of DNA

Exact positions of 5-methylcytosine (m(5)C) on a single strand of DNA can be determined by bisulfite genomic sequencing (BGS). Treatment with bisulfite ion preferentially deaminates unmethylated cytosines, which are then converted to uracil upon desulfonation. Amplifying regions of interest from deaminated DNA and sequencing products cloned from amplicons permits determination of methylation at single-nucleotide resolution along single DNA molecules, which is not possible with other methylation analysis techniques. This unit describes a BGS technique suitable for most DNA sources, including formaldehyde-fixed tissue. Considerations for experimental design and common sources of error are discussed.