scholarly journals NanoMethViz: An R/Bioconductor package for visualizing long-read methylation data

2021 ◽  
Vol 17 (10) ◽  
pp. e1009524
Author(s):  
Shian Su ◽  
Quentin Gouil ◽  
Marnie E. Blewitt ◽  
Dianne Cook ◽  
Peter F. Hickey ◽  
...  

A key benefit of long-read nanopore sequencing technology is the ability to detect modified DNA bases, such as 5-methylcytosine. The lack of R/Bioconductor tools for the effective visualization of nanopore methylation profiles between samples from different experimental groups led us to develop the NanoMethViz R package. Our software can handle methylation output generated from a range of different methylation callers and manages large datasets using a compressed data format. To fully explore the methylation patterns in a dataset, NanoMethViz allows plotting of data at various resolutions. At the sample-level, we use dimensionality reduction to look at the relationships between methylation profiles in an unsupervised way. We visualize methylation profiles of classes of features such as genes or CpG islands by scaling them to relative positions and aggregating their profiles. At the finest resolution, we visualize methylation patterns across individual reads along the genome using the spaghetti plot and heatmaps, allowing users to explore particular genes or genomic regions of interest. In summary, our software makes the handling of methylation signal more convenient, expands upon the visualization options for nanopore data and works seamlessly with existing methylation analysis tools available in the Bioconductor project. Our software is available at https://bioconductor.org/packages/NanoMethViz.

2021 ◽  
Author(s):  
Shian Su ◽  
Quentin Gouil ◽  
Marnie E. Blewitt ◽  
Dianne Cook ◽  
Peter F. Hickey ◽  
...  

AbstractMotivationA key benefit of long-read nanopore sequencing technology is the ability to detect modified DNA bases, such as 5-methylcytosine. Tools for effective visualization of data generated by this platform to assess changes in methylation profiles between samples from different experimental groups remains a challenge.ResultsTo make visualization of methylation changes more straightforward, we developed the R/Bioconductor package NanoMethViz. Our software can handle methylation calls generated from a range of different methylation callers and manages large datasets using a compressed data format. To fully explore the methylation patterns in a dataset, NanoMethViz allows plotting of data at various resolutions. At the sample-level, we use multidimensional scaling to look at the relationships between methylation profiles in an unsupervised way. We visualize methylation profiles of classes of features such as genes or CpG islands by scaling them to relative positions and aggregating their profiles. At the finest resolution, we visualize methylation patterns across individual reads along the genome using the spaghetti plot, allowing users to explore particular genes or genomic regions of interest.In summary, our software makes the handling of methylation signal more convenient, expands upon the visualization options for nanopore data and works seamlessly with existing methylation analysis tools available in the Bioconductor project. Our software is available at https://bioconductor.org/packages/NanoMethViz.


Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. 653-653 ◽  
Author(s):  
Ying Qu ◽  
Andreas Lennartsson ◽  
Verena I. Gaidzik ◽  
Stefan Deneberg ◽  
Sofia Bengtzén ◽  
...  

Abstract Abstract 653 DNA methylation is involved in multiple biologic processes including normal cell differentiation and tumorigenesis. In AML, methylation patterns have been shown to differ significantly from normal hematopoietic cells. Most studies of DNA methylation in AML have previously focused on CpG islands within the promoter of genes, representing only a very small proportion of the DNA methylome. In this study, we performed genome-wide methylation analysis of 62 AML patients with CN-AML and CD34 positive cells from healthy controls by Illumina HumanMethylation450K Array covering 450.000 CpG sites in CpG islands as well as genomic regions far from CpG islands. Differentially methylated CpG sites (DMS) between CN-AML and normal hematopoietic cells were calculated and the most significant enrichment of DMS was found in regions more than 4kb from CpG Islands, in the so called open sea where hypomethylation was the dominant form of aberrant methylation. In contrast, CpG islands were not enriched for DMS and DMS in CpG islands were dominated by hypermethylation. DMS successively further away from CpG islands in CpG island shores (up to 2kb from CpG Island) and shelves (from 2kb to 4kb from Island) showed increasing degree of hypomethylation in AML cells. Among regions defined by their relation to gene structures, CpG dinucleotide located in theoretic enhancers were found to be the most enriched for DMS (Chi χ2<0.0001) with the majority of DMS showing decreased methylation compared to CD34 normal controls. To address the relation to gene expression, GEP (gene expression profiling) by microarray was carried out on 32 of the CN-AML patients. Totally, 339723 CpG sites covering 18879 genes were addressed on both platforms. CpG methylation in CpG islands showed the most pronounced anti-correlation (spearman ρ =-0.4145) with gene expression level, followed by CpG island shores (mean spearman rho for both sides' shore ρ=-0.2350). As transcription factors (TFs) have shown to be crucial for AML development, we especially studied differential methylation of an unbiased selection of 1638 TFs. The most enriched differential methylation between CN-AML and normal CD34 positive cells were found in TFs known to be involved in hematopoiesis and with Wilms tumor protein-1 (WT1), activator protein 1 (AP-1) and runt-related transcription factor 1 (RUNX1) being the most differentially methylated TFs. The differential methylation in WT 1 and RUNX1 was located in intragenic regions which were confirmed by pyro-sequencing. AML cases were characterized with respect to mutations in FLT3, NPM1, IDH1, IDH2 and DNMT3A. Correlation analysis between genome wide methylation patterns and mutational status showed statistically significant hypomethylation of CpG Island (p<0.0001) and to a lesser extent CpG island shores (p<0.001) and the presence of DNMT3A mutations. This links DNMT3A mutations for the first time to a hypomethylated phenotype. Further analyses correlating methylation patterns to other clinical data such as clinical outcome are ongoing. In conclusion, our study revealed that non-CpG island regions and in particular enhancers are the most aberrantly methylated genomic regions in AML and that WT 1 and RUNX1 are the most differentially methylated TFs. Furthermore, our data suggests a hypomethylated phenotype in DNMT3A mutated AML. Disclosures: No relevant conflicts of interest to declare.


Genes ◽  
2018 ◽  
Vol 9 (10) ◽  
pp. 500
Author(s):  
Juan A. Subirana ◽  
Xavier Messeguer

Repetitive genome regions have been difficult to sequence, mainly because of the comparatively small size of the fragments used in assembly. Satellites or tandem repeats are very abundant in nematodes and offer an excellent playground to evaluate different assembly methods. Here, we compare the structure of satellites found in three different assemblies of the Caenorhabditis elegans genome: the original sequence obtained by Sanger sequencing, an assembly based on PacBio technology, and an assembly using Nanopore sequencing reads. In general, satellites were found in equivalent genomic regions, but the new long-read methods (PacBio and Nanopore) tended to result in longer assembled satellites. Important differences exist between the assemblies resulting from the two long-read technologies, such as the sizes of long satellites. Our results also suggest that the lengths of some annotated genes with internal repeats which were assembled using Sanger sequencing are likely to be incorrect.


2018 ◽  
Vol 4 (11) ◽  
pp. eaau6986 ◽  
Author(s):  
Lu Wang ◽  
Patrick A. Ozark ◽  
Edwin R. Smith ◽  
Zibo Zhao ◽  
Stacy A. Marshall ◽  
...  

The tet methylcytosine dioxygenase 2 (TET2) enzyme catalyzes the conversion of the modified DNA base 5-methylcytosine to 5-hydroxymethylcytosine. TET2 is frequently mutated or dysregulated in multiple human cancers, and loss of TET2 is associated with changes in DNA methylation patterns. Here, using newly developed TET2-specific antibodies and the estrogen response as a model system for studying the regulation of gene expression, we demonstrate that endogenous TET2 occupies active enhancers and facilitates the proper recruitment of estrogen receptor α (ERα). Knockout of TET2 by CRISPR-CAS9 leads to a global increase of DNA methylation at enhancers, resulting in attenuation of the estrogen response. We further identified a positive feedback loop between TET2 and ERα, which further requires MLL3 COMPASS at these enhancers. Together, this study reveals an epigenetic axis coordinating a transcriptional program through enhancer activation via DNA demethylation.


1993 ◽  
Vol 13 (9) ◽  
pp. 5538-5548
Author(s):  
Y C Choi ◽  
C B Chae

In contrast to many other genes containing a CpG island, the testis-specific H2B (TH2B) histone gene exhibits tissue-specific methylation patterns in correlation with gene activity. Characterization of the methylation patterns within a 20-kb segment containing the TH2A and TH2B genes in comparison with that in a somatic histone cluster revealed that: (i) the germ cell-specific unmethylated domain of the TH2A and TH2B genes is defined as a small region surrounding the CpG islands of the TH2A and TH2B genes and (ii) somatic histone genes are unmethylated in both liver and germ cells, like other genes containing CpG islands, whereas flanking sequences are methylated. Transfection of in vitro-methylated TH2B, somatic H2B, and mouse metallothionein I constructs into F9 embryonal carcinoma cells revealed that the CpG islands of the TH2A and TH2B genes were demethylated like those of the somatic H2A and H2B genes and the metallothionein I gene. The demethylation of those CpG islands became significantly inefficient at a high number of integrated copies and a high density of methylated CpG dinucleotides. In contrast, three sites in the somatic histone cluster, of which two sites are located in the long terminal repeat of an endogenous retrovirus-like sequence, were efficiently demethylated even at a high copy number and a high density of methylated CpG dinucleotides. These results suggest two possible mechanisms for demethylation in F9 cells and methylation of CpG islands of the TH2A and TH2B genes at the postblastula stage during embryogenesis.


2016 ◽  
Author(s):  
J. Alakuijala ◽  
Z. Szabadka
Keyword(s):  

1993 ◽  
Vol 13 (9) ◽  
pp. 5538-5548 ◽  
Author(s):  
Y C Choi ◽  
C B Chae

In contrast to many other genes containing a CpG island, the testis-specific H2B (TH2B) histone gene exhibits tissue-specific methylation patterns in correlation with gene activity. Characterization of the methylation patterns within a 20-kb segment containing the TH2A and TH2B genes in comparison with that in a somatic histone cluster revealed that: (i) the germ cell-specific unmethylated domain of the TH2A and TH2B genes is defined as a small region surrounding the CpG islands of the TH2A and TH2B genes and (ii) somatic histone genes are unmethylated in both liver and germ cells, like other genes containing CpG islands, whereas flanking sequences are methylated. Transfection of in vitro-methylated TH2B, somatic H2B, and mouse metallothionein I constructs into F9 embryonal carcinoma cells revealed that the CpG islands of the TH2A and TH2B genes were demethylated like those of the somatic H2A and H2B genes and the metallothionein I gene. The demethylation of those CpG islands became significantly inefficient at a high number of integrated copies and a high density of methylated CpG dinucleotides. In contrast, three sites in the somatic histone cluster, of which two sites are located in the long terminal repeat of an endogenous retrovirus-like sequence, were efficiently demethylated even at a high copy number and a high density of methylated CpG dinucleotides. These results suggest two possible mechanisms for demethylation in F9 cells and methylation of CpG islands of the TH2A and TH2B genes at the postblastula stage during embryogenesis.


2020 ◽  
Vol 2 (3) ◽  
Author(s):  
Batuhan Cakir ◽  
Martin Prete ◽  
Ni Huang ◽  
Stijn van Dongen ◽  
Pinar Pir ◽  
...  

Abstract In the last decade, single cell RNAseq (scRNAseq) datasets have grown in size from a single cell to millions of cells. Due to its high dimensionality, it is not always feasible to visualize scRNAseq data and share it in a scientific report or an article publication format. Recently, many interactive analysis and visualization tools have been developed to address this issue and facilitate knowledge transfer in the scientific community. In this study, we review several of the currently available scRNAseq visualization tools and benchmark the subset that allows to visualize the data on the web and share it with others. We consider the memory and time required to prepare datasets for sharing as the number of cells increases, and additionally review the user experience and features available in the web interface. To address the problem of format compatibility we have also developed a user-friendly R package, sceasy, which allows users to convert their own scRNAseq datasets into a specific data format for visualization.


2013 ◽  
Vol 42 (4) ◽  
pp. 2235-2244 ◽  
Author(s):  
Jan O. Haerter ◽  
Cecilia Lövkvist ◽  
Ian B. Dodd ◽  
Kim Sneppen

Abstract Inheritance of 5-methyl cytosine modification of CpG (CG/CG) DNA sequences is needed to maintain early developmental decisions in vertebrates. The standard inheritance model treats CpGs as independent, with methylated CpGs maintained by efficient methylation of hemimethylated CpGs produced after DNA replication, and unmethylated CpGs maintained by an absence of de novo methylation. By stochastic simulations of CpG islands over multiple cell cycles and systematic sampling of reaction parameters, we show that the standard model is inconsistent with many experimental observations. In contrast, dynamic collaboration between CpGs can provide strong error-tolerant somatic inheritance of both hypermethylated and hypomethylated states of a cluster of CpGs, reproducing observed stable bimodal methylation patterns. Known recruitment of methylating enzymes by methylated CpGs could provide the necessary collaboration, but we predict that recruitment of demethylating enzymes by unmethylated CpGs strengthens inheritance and allows CpG islands to remain hypomethylated within a sea of hypermethylation.


Sign in / Sign up

Export Citation Format

Share Document