Methylation plotter: a web tool for dynamic visualization of DNA methylation data
© Mallona et al.; licensee BioMed Central Ltd. 2014
Received: 27 March 2014
Accepted: 4 June 2014
Published: 7 June 2014
Methylation plotter is a Web tool that allows the visualization of methylation data in a user-friendly manner and with publication-ready quality. The user is asked to introduce a file containing the methylation status of a genomic region. This file can contain up to 100 samples and 100 CpGs. Optionally, the user can assign a group for each sample (i.e. whether a sample is a tumoral or normal tissue). After the data upload, the tool produces different graphical representations of the results following the most commonly used styles to display this type of data. They include an interactive plot that summarizes the status of every CpG site and for every sample in lollipop or grid styles. Methylation values ranging from 0 (unmethylated) to 1 (fully methylated) are represented using a gray color gradient. A practical feature of the tool allows the user to choose from different types of arrangement of the samples in the display: for instance, sorting by overall methylation level, by group, by unsupervised clustering or just following the order in which data were entered.
In addition to the detailed plot, Methylation plotter produces a methylation profile plot that summarizes the status of the scrutinized region, a boxplot that sums up the differences between groups (if any) and a dendrogram that classifies the data by unsupervised clustering. Coupled with this analysis, descriptive statistics and testing for differences at both CpG and group levels are provided.
The implementation is based in R/shiny, providing a highly dynamic user interface that generates quality graphics without the need of writing R code. Methylation plotter is freely available at http://gattaca.imppc.org:3838/methylation_plotter/.
KeywordsMethylation plot Methylation visualization R/shiny Lollipop plot
Cytosine methylation in CpG dinucleotides is an important mechanism involved in the regulation of multiple biological processes including pathological conditions [1–3]. While there is a wide range of methodologies to evaluate DNA methylation, bisulfite-treated DNA sequencing is the gold standard to determine DNA methylation at the single CpG level [1, 4, 5]. The functional implications of DNA methylation states are often determined by the CpG profile but at the regional level rather than by a single CpG site. Therefore, the interpretation and application of this sort of data require further analysis that is highly benefited by the implementation of visualization tools.
While some software tools to analyze and visually represent DNA methylation data have been published (reviewed in ), its use by wet lab users is often limited. On the other hand, popular spreadsheet tools like Excel are unable to generate lollipop plots by default. Even more, the Excel-based solutions perform poorly for repetitive tasks: in an automated analysis context, programmatic approaches are less error prone and more reproducible .
Specialized tools have been developed to work with converted bisulfite sequence files and to explore methylation trends, but are highly dependent on the operating system: MethTools, ) is Unix-based, and CpG Analyzer  or CpG PatternFinder ) run under Windows. MethDB  offers a web tool and thus is platform-independent, but is designed as a methylation data provider rather than a graphical tool. BiQ Analyzer [11, 12] and QUMA  provide web tools that plot lollipop-like graphics; however, they are rather devoted to clonal analysis, assessing the methylation status as a categorical variable (either methylated or unmethylated). Hence, a platform-independent tool to visualize continuous methylation data, as those produced by direct bisulfite sequencing or microarray platforms, is needed.
The interactive web application is written using shiny, an R framework that couples the R-based statistics computation and graphics generation to the rendering of a Web-based user interface . This technology allows to take advantage of the R power in an easy-to-use frontend. As the application is hosted in a remote server, the user does not need to consume local resources and just requires a Web browser to use the tool. User data is removed from the server as soon as the browser session terminates.
Results and discussion
Methylation plotter is an interactive application that allows rapid and easy generation of customized plots and statistical summaries of methylation data. The user is asked to upload a tab-separated file describing the status of up to 100 CpGs in up to 100 different samples as well the group each sample belongs to. The application generates an interactive plot that summarizes the status of every CpG site and for every sample in lollipop or grid styles. Methylation values ranging from 0 (unmethylated) to 1 (fully methylated) are represented using a gray color gradient.
Altogether, Methylation plotter provides descriptive statistics and basic non-parametric variance analysis (Kruskal-Wallis tests). For each sample, a data table summarizing the mean, standard deviation, minimum and maximum, and number of not available positions (NAs) is produced. The same descriptive statistics are produced for each CpG and, if the input data is ascribed to two or more groups, each CpG is tested for equality using the Kruskal-Wallis test.
All the figures are available to download as either raster (PNG) or vector format files (PDF), whereas statistical reports are served as spreadsheets (tab-separated values).
In summary, Methylation plotter is a user-friendly tool that produces ready-to-use plots and summary data required by most wet lab users analyzing DNA methylation. The application is freely accessible at http://gattaca.imppc.org:3838/methylation_plotter/.
Availability and requirements
Project name: Methylation plotter
Project home page: http://sourceforge.net/projects/methylationplotter
Operating system(s): Platform independent
Programming language: R/shiny
Other requirements: None
License: GPL v2
Any restrictions to use by non-academics: None
We thank Judith Flo for her excellent technical assistance and the members of the lab for comments and beta-testing of the tool. AD was supported in part by contract PTC2011-1091 from the Ministry of Economy and Knowledge. This work was supported by grants from the Ministry of Economy and Knowledge (SAF2011/23638), and the Generalitat de Catalunya (2009 SGR 1356).
- Esteller M:Cancer epigenomics: dna methylomes and histone-modification maps. Nat Rev Genet. 2007, 8 (4): 286-298. 10.1038/nrg2005.View ArticlePubMedGoogle Scholar
- Jones PA, Baylin SB:The epigenomics of cancer. Cell. 2007, 128 (4): 683-692. 10.1016/j.cell.2007.01.029.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki MM, Bird A:Dna methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008, 9 (6): 465-476. 10.1038/nrg2341.View ArticlePubMedGoogle Scholar
- Peinado MA, Jordà M:Methods for dna methylation analysis and applications in colon cancer. Mutat Res Fund Mol Mech Mutagen. 2010, 693 (1): 84-93.Google Scholar
- Laird PW:Principles and challenges of genome-wide dna methylation analysis. Nat Rev Genet. 2010, 11 (3): 191-203. 10.1038/nrg2732.View ArticlePubMedGoogle Scholar
- Mesirov JP:Computer science. Accessible reproducible research. Science. 2010, 327 (5964): 415-416. 10.1126/science.1179653.View ArticlePubMedGoogle Scholar
- Grunau C, Schattevoy R, Mache N, Rosenthal A:Methtools—a toolbox to visualize and analyze dna methylation data. Nucleic Acids Res. 2000, 28 (5): 1053-1058. 10.1093/nar/28.5.1053.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu Y, Manoharan HT, Pitot HC:Cpg analyzer, a windows-based utility program for investigation of dna methylation. Biotechniques. 2005, 39 (5): 656-10.2144/000112053.View ArticlePubMedGoogle Scholar
- Xu Y-H, Manoharan HT, Pitot HC:Cpg patternfinder: a windows®;-based utility program for easy and rapid identification of the cpg methylation status of dna. Biotechniques. 2007, 43: 334-342. 10.2144/000112537.View ArticlePubMedGoogle Scholar
- Grunau C, Renault E, Rosenthal A, Roizes G:Methdb—a public database for dna methylation data. Nucleic Acids Res. 2001, 29 (1): 270-274. 10.1093/nar/29.1.270.PubMed CentralView ArticlePubMedGoogle Scholar
- Bock C, Reither S, Mikeska T, Paulsen M, Walter J, Lengauer T:Biq analyzer: visualization and quality control for dna methylation data from bisulfite sequencing. Bioinformatics. 2005, 21 (21): 4067-4068. 10.1093/bioinformatics/bti652.View ArticlePubMedGoogle Scholar
- Lutsik P, Feuerbach L, Arand J, Lengauer T, Walter J, Bock C:Biq analyzer ht: locus-specific analysis of dna methylation by high-throughput bisulfite sequencing. Nucleic Acids Res. 2011, 39 (suppl 2): 551-556.View ArticleGoogle Scholar
- Kumaki Y, Oda M, Okano M:Quma: quantification tool for methylation analysis. Nucleic Acids Res. 2008, 36 (suppl 2): 170-175.View ArticleGoogle Scholar
- RStudio Inc.:Shiny: Web Application Framework for R. 2013, R package version 0.8.0. [http://CRAN.R-project.org/package=shiny].Google Scholar
- Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, Wu B, Doucet D, Thomas NJ, Wang Y, Vollmer E, Goldmann T, Seifart C, Jiang W, Barker DL, Chee MS, Floros J, Fan JB:High-throughput dna methylation profiling using universal bead arrays. Genome Res. 2006, 16 (3): 383-393. 10.1101/gr.4410706.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.