The PinkThing for analysing ChIP profiling data in their genomic context
© Nielsen et al.; licensee BioMed Central Ltd. 2013
Received: 14 December 2012
Accepted: 12 March 2013
Published: 4 April 2013
Current epigenetic research makes frequent use of whole-genome ChIP profiling for determining the in vivo binding of proteins, e.g. transcription factors and histones, to DNA. Two important and recurrent questions for these large scale analyses are: 1) What is the genomic distribution of a set of binding sites? and 2) Does this genomic distribution differ significantly from another set of sites?
We exemplify the functionality of the PinkThing by analysing a ChIP profiling dataset of cohesin binding sites. We show the subset of cohesin sites with no CTCF binding have a characteristic genomic distribution different from the set of all cohesin sites.
The PinkThing is a web application for fast and easy analysis of the context of genomic loci, such as peaks from ChIP profiling experiments. The output of the PinkThing analysis includes: categorisation of position relative to genes (intronic, exonic, 5’ near, 3’ near 5’ far, 3’ far and distant), distance to the closest annotated 3’ and 5’ end of genes, direction of transcription of the nearest gene, and the option to include other genomic elements like ESTs and CpG islands. The PinkThing enables easy statistical comparison between experiments, i.e. experimental versus background sets, reporting over- and underrepresentation as well as p-values for all comparisons. Access and use of the PinkThing is free and open (without registration) to all users via the website: http://pinkthing.cmbi.ru.nl
Transcription factor (TF) binding sites can be identified in vivo using the emerging technologies for ChIP profiling such as ChIP-on-chip and ChIP-sequencing. These experiments locate hundreds to tens of thousands TF binding sites, which subsequently have to be validated and analysed for biological function. Certain initial analyses of TF binding sites have already become commonplace: mapping of the binding sites to the genome, detecting closest genes to the binding sites, categorising binding sites by their position relative to the genes, and the subsequent GO analysis of the genes closest to the binding sites. The same type of analyses apply to other regions found by ChIP profiling, e.g. with specific histone modification patterns or DNA hyper- or hypomethylated loci. The PinkThing gathers all these analyses into one single user-friendly tool, standardises the characterisation of genomic locations, uses up-to-date ENSEMBL gene annotation, and at the same time enables statistically sound comparisons at each step of the analysis. Although it is possible to perform similar analysis using bioinformatic packages like Taverna  or Galaxy , these packages require installation, initialisation and scripting (Taverna) or a certain level of statistics skills from the user (e.g. Galaxy). Another more specialised tool is the HyperBrowser , which is based on Galaxy and designed to incorporate functions to query and correlate annotation along the genome. The versatility of this tool is impressive, but it comes at the cost of usability. Compared to the PinkThing, the HyperBrowser has a complex interface that takes time to learn to use efficiently. A more accessible tool is GREAT  which provides a web interface for analysis, with a focus specifically on cis-regulatory regions for human, mouse and zebrafish.
With the PinkThing tool the genomic analyses with statistics are available and ready-to-use directly from the website by a simple upload of a file containing the coordinates of genomic regions of interest. In addition, the PinkThing provides the option of supplying an appropriate background distribution to be used for comparisons, e.g. the set of all binding sites of a ChIP-seq experiment when examining a subset of these sites. Furthermore, the PinkThing analysis of genomic distributions supports all species that are annotated in Ensembl. PinkThing has already been successfully used for the analysis of ChIP profiling data in a wide variety of data, including [5–17].
Through the PinkThing web page, the user uploads the genomic sites of interest using a standard format for genomic locations (BED format). PinkThing compares the sites to the Ensembl gene annotation and optionally to Ensembl CpG islands, ESTs and regulatory features. The initial results include basic statistics of the genomic distribution of the uploaded positions:
histograms of distances to the 5’ and the 3’ ends of genes,
a barplot and a pie chart of the genomic distribution relative to Ensembl genes/CpG islands/ESTs (Figure 1a).
a barplot of the genomic distribution.
Every output plot is clickable for download of an SVG (Scalable Vector Graphics) version of the plot for publication. Subsequent analyses can be accessed from the results page and include:
transcription start sites (TSS) within the surrounding regions,
the conservation score (GERP ) for these regions,
comparisons to other sets of genomic regions,
retrieval of the reference DNA sequence for the uploaded regions,
GO annotation and GO enrichment analysis of the nearest genes (only available for human and mouse).
As an example, we considered the CTCF (CCCTC binding factor) and cohesin binding sites as determined using ChIP-seq in a study by Schmidt et al.. In this study, Schmidt et al. also discovered a subset of cohesin binding sites that did not overlap with CTCF binding sites. Using the PinkThing we compared the genomic distribution of the set of cohesin sites with the set of cohesin sites not bound by CTCF (cohesin-no-CTCF sites). We chose the cutoffs for ‘near’ and ‘far’ categories in our analysis to be 1000bp and 3000bp respectively, to match the categories presented in their paper  and examined PinkThing-determined categorisation and genomic distribution of the sites (Figure 1a). The categorisation of the subset cohesin-no-CTCF is very informative when viewed in relation to the background set of all cohesin binding sites: The PinkThing provides the relative sizes of each of the categories for the two sets of sites (Figure 1b), the Brandt-Snedecor test statistic for comparison of the two distributions, as well as p-values for whether the differences in distributions are significant per category (Figure 1c) as well as a barplot showing the frequencies of the two distributions grouped per category (Figure 1d).
The Brandt-Snedecor test indicates that the genomic distribution of cohesin and cohesin-no-CTCF sites are significantly different (p<2.2e−16) and the individual category comparison shows the cohesin-no-CTCF sites are significantly overrepresented in promoter regions (5’ near) and exon regions, and underrepresented in 3’ regions as compared to the set of all cohesin sites (Figure 1c). The Gene Ontology (GO) result from the PinkThing (Additional file 1: noCTCF-vs-cohesin.xls) indicates that the cohesin-no-CTCF binding sites are, relative to all cohesin sites, overrepresented (adjusted p<10e−3) in the following GO categories: nucleic acid binding (GO:0003676), transcription regulator activity (GO:0030528), regulation of metabolic process (GO:0019222) (See Supplementary material: noCTCF-vs-cohesin-GO.xls). The difference in genomic distributions suggests a difference in function, which was confirmed by Schmidt et al. who showed that the cohesin-no-CTCF sites are enriched for DNA binding sites of tissue-specific transcription factors, and further explored this fact in a follow-up study .
Lookup and data retrieval from the Ensembl database was implemented using the Ensembl Perl API http://www.ensembl.org/info/docs/api/core/core\_tutorial.html.
Statistical tests and generation of plots and figures on the PinkThing website are implemented in R  and bioperl . Detailed description of plots and their interpretation is in the PinkThing manual available from the PinkThing website.
For each imported locus, PinkThing determines the distance to the closest gene. Positions that fall into multiple categories are assigned to the highest ranking category in the order: exon, intron, 5’ near, 3’ near, 5’ far, 3’ far and distant. As the surrounding sequence may contain overlapping gene annotations, the gene annotations are merged such that each position corresponds to a unique genome category. The limits for “near” and “far” may be chosen by the user, with default values being 5kb and 25kb, respectively.
Comparison of category distributions
The initial statistical test is the Brandt-Snedecor test for homogeneity of k binomial distributions . The test assumes i) that the observations are independent; ii) that the set chosen as background distribution is at least twice the size as the set chosen as sample distribution and iii) that the k categories are exhaustive and mutually exclusive. Our implementation of the Brandt-Snedecor test does not test whether there is an overlap of actual observations in the two sets of the comparison. The test compares the distribution of observations over categories, indifferent to the identity of the observation. In our example above we show how we apply the comparison between a set of sites (cohesin binding sites) and a subset of those sites (cohesin-noCTCF sites), but the test can also be applied to compare two disparate sets of sites in the genome.
To determine whether each of the genomic categories is significantly overrepresented, we apply the log likelihood G-test for independence  per category, comparing each category against the union of the other categories.
The PinkThing uses Ontologizer  to perform GO statistics on the set of genes that are closest to the uploaded set of sites. The output includes both the graphical Ontologizer visualization of overrepresented categories as well as a table with p-values of all terms. The GO statistics can either be calculated against the whole genome as background or against the genes corresponding to another uploaded set. Currently the GO ontology analysis is available for human and mouse.
Determination of the annotation context of genomic loci is an indispensable foundation for their functional analysis. To allow distributions of genomic categories to be compared within the same genome assembly, the implementation of categorisation in PinkThing consistently assigns exactly one category to any genomic position, thus providing consistent and reproducible analysis. The PinkThing categorisation into consistent genomic categories allows the comparison with a background or reference distribution. This comparison is essential when examining distributions where the categories vary in size and abundance dependent on the context of the genome (the species) and the experiment (e.g. selecting for specific genomic elements). PinkThing makes it easy to compare results from a specific experiment with a chosen background distribution, thus increasing the value and confidence in interpreting results.
PinkThing is a collection of the most common genomic analyses related to genomic context, combined in a simple point and click web interface. With no prerequisites other than obtaining a dataset, PinkThing provides easy access to sound statistical analysis of genomic location data.
Availability and requirements
Access and use of the PinkThing is free and open (without registration) to all users via the website: http://pinkthing.cmbi.ru.nl
Project name: The PinkThing for analysing ChIP profiling data in their genomic context
Project home page: http://www.bioinformatics.org/websvn/listing.php?repname=pinkthing
Operating systems: Usage of web tool is platform independent, access is available via all major web browsers at http://pinkthing.cmbi.ru.nl.
Code repository (SVN): http://www.bioinformatics.org/websvn/listing.php?repname=pinkthing
License: GNU GPL
The two data sets with cohesin binding sites with and without CTCF are available for download from the front page of the PinkThing and included with the additional files for this article (Additional file 2: all cohesin binding sites, cohesin.bed; Additional file 3: cohesin binding sites without CTCF, no_CTCF_cohesin.bed).
Go to http://pinkthing.cmbi.ru.nland select Ensembl version 53.
Then enter the desired cutoff for near and far positions as 1000 and 3000 respectively.
Upload the two files one at a time by selecting the file location and click ‘upload and calculate’.
To perform the comparison, from the start page of the PinkThing, select the two uploaded files in the box ‘Compare genomic distributions’, choosing the cohesin sites with no CTCF binding as the sample file (foreground) and choosing the set of all cohesin binding sites as sample space(background).
Click ‘Compare’ to view the results.
To find overrepresented ontologies, from the start page of the PinkThing, select the two uploaded files in the respected boxes for Sample space (all cohesin binding sites) and sample file (cohesin binding sites without CTCF) and click Ontologize.
The authors would like to thank Willem Welboren, Arjen Brinkman, and the rest of the PinkThing users who provided valuable feedback during development.
This work was supported by HEROIC, an Integrated Project funded by the European Union under the 6th Framework Program (LSHG-CT-2005-018883).
- Hull D, Wolstencroft K, Stevens R, Goble C, Pocock MR, Li P, Oinn T: Taverna: a tool for building and running workflows of services. Nucleic Acids Res. 2006, 34: (Web Server issue):W729–W732 [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1538887&tool=pmcentrez&rendertype=abstract]Google Scholar
- Goecks J, Nekrutenko A, Taylor J, Galaxy Team T: Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010, 11 (8): R86-10.1186/gb-2010-11-8-r86. [http://genomebiology.com/2010/11/8/R86]PubMedPubMed CentralView ArticleGoogle Scholar
- Sandve GK, Gundersen S, Rydbeck H, Glad IK, Holden L, Holden M, Clancy T, Ferkingstad E, Johansen M, Nygaard V, T ϕ stesen E, Frigessi A, Hovig E, Liest ϕ l K: The Genomic HyperBrowser: inferential genomics at the sequence level. Genome Biol. 2010, 11 (12): :R121-[http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3046481&tool=pmcentrez&rendertype=abstract]View ArticleGoogle Scholar
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G: GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010, 28 (5): 495-501. 10.1038/nbt.1630. [http://www.ncbi.nlm.nih.gov/pubmed/20436461]PubMedView ArticleGoogle Scholar
- Nielsen R, Pedersen TA, Hagenbeek D, Moulos P, Siersbaek R, Megens E, Denissov S, B ϕ rgesen M, Francoijs KJ, Mandrup S, Stunnenberg HG: Genome-wide profiling of PPARgamma: RXR and RNA polymerase II occupancy reveals temporal activation of distinct metabolic pathways and changes in RXR dimer composition during adipogenesis. Genes Dev. 2008, 22 (21): 2953-2967. 10.1101/gad.501108. [http://genesdev.cshlp.org/cgi/content/abstract/22/21/2953]PubMedPubMed CentralView ArticleGoogle Scholar
- Sanderson L: Exploring the activation and function of PPARα and PPARβ/δ using genomics. Phd thesis. Wageningen 2009, [http://library.wur.nl/WebQuery/wda/lang/1918369]Google Scholar
- Boyd M, Hansen M, Jensen TGK, Perearnau A, Olsen AK, Bram LL, Bak M, Tommerup N, Olsen Jr, Troelsen JT: Genome-wide analysis of CDX2 binding in intestinal epithelial cells (Caco-2). J Biol Chem. 2010, 285 (33): 25115-25125. 10.1074/jbc.M109.089516. [http://www.jbc.org/cgi/content/abstract/285/33/25115]PubMedPubMed CentralView ArticleGoogle Scholar
- van der Meer DLM, Degenhardt T, Väisänen S, de Groot PJ, Heinäniemi M, de Vries SC, Müller M, Carlberg C, Kersten S: Profiling of promoter occupancy by PPARalpha in human hepatoma cells via ChIP-chip analysis. Nucleic Acids Res. 2010, 38 (9): 2839-2850. 10.1093/nar/gkq012. [http://nar.oxfordjournals.org/cgi/content/abstract/38/9/2839]PubMedPubMed CentralView ArticleGoogle Scholar
- Vermeulen M, Eberl HC, Matarese F, Marks H, Denissov S, Butter F, Lee KK, Olsen JV, Hyman AA, Stunnenberg HG: Quantitative interaction proteomics and genome-wide profiling of epigenetic histone marks and their readers. Cell. 2010, 142 (6): 967-980. 10.1016/j.cell.2010.08.020. [http://dx.doi.org/10.1016/j.cell.2010.08.020]PubMedView ArticleGoogle Scholar
- Rada-Iglesias A, Bajpai R, Swigut T, Brugmann SA, Flynn RA, Wysocka J: A unique chromatin signature uncovers early developmental enhancers in humans. Nature. 2010, [http://www.ncbi.nlm.nih.gov/pubmed/21160473]Google Scholar
- Rao NAS, Mccalman MT, Moulos P, Francoijs Kj, Chatziioannou A, Kolisis FN, Alexis MN, Mitsiou DJ, Stunnenberg HG: Coactivation of GR and NFKB alters the repertoire of their binding sites and target genes. Genome Res. 2011, Sep, 21 (9): 1404-16. 10.1101/gr.118042.110. http://dx.doi.org/10.1101/gr.118042.110. Epub 2011 Jul 12.PubMedPubMed CentralView ArticleGoogle Scholar
- Boergesen M, Pedersen TAs, van Heeringen SJ, Hagenbeek D, Bindesb ϕ ll C, Caron S, Lalloyer F, Steffensen KR, Nebb HI, Gustafsson JAk, Stunnenberg HG, Staels B, Mandrup S: Genome-wide profiling of liver X receptor, retinoid X receptor, and peroxisome proliferator-activated receptor α in mouse liver reveals extensive sharing of binding sites. Mol Cell Biol. 2012, 32 (4): 852-867. 10.1128/MCB.06175-11. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3272984&tool=pmcentrez&rendertype=abstract]PubMedPubMed CentralView ArticleGoogle Scholar
- Nitzsche A, Paszkowski-Rogacz M, Matarese F, Janssen-Megens EM, Hubner NC, Schulz H, de Vries I, Ding L, Huebner N, Mann M, Stunnenberg HG, Buchholz F: RAD21 cooperates with pluripotency transcription factors in the maintenance of embryonic stem cell identity. PloS one. 2011, 6 (5): :e19470-[http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3093395&tool=pmcentrez&rendertype=abstract]View ArticleGoogle Scholar
- Costessi A, Mahrour N, Tijchon E, Stunnenberg R, Stoel Ma, Jansen PW, Sela D, Martin-Brown S, Washburn MP, Florens L, Conaway JW, Conaway RC, Stunnenberg HG: The tumour antigen PRAME is a subunit of a Cul2 ubiquitin ligase and associates with active NFY promoters. EMBO J. 2011, 30 (18): 3786-3798. 10.1038/emboj.2011.262. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3173790&tool=pmcentrez&rendertype=abstract]PubMedPubMed CentralView ArticleGoogle Scholar
- Nielsen FGG, Markus KGt, Friborg RMl, Favrholdt LM, Stunnenberg HG, Huynen M: CATCHprofiles: clustering and alignment tool for ChIP profiles. PloS one. 2012, 7: e28272-10.1371/journal.pone.0028272. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3251562&tool=pmcentrez&rendertype=abstract]PubMedPubMed CentralView ArticleGoogle Scholar
- Yu CY, Mayba O, Lee JV, Tran J, Harris C, Speed TP, Wang JC: Genome-wide analysis of glucocorticoid receptor binding regions in adipocytes reveal gene network involved in triglyceride homeostasis. PloS one. 2010, 5 (12): e15188-10.1371/journal.pone.0015188. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3004788&tool=pmcentrez&rendertype=abstract]PubMedPubMed CentralView ArticleGoogle Scholar
- Kramer JM, Kochinke K, Oortveld MAW, Marks H, Kramer D, de Jong EK, Asztalos Z, Westwood JT, Stunnenberg HG, Sokolowski MB, Keleman K, Zhou H, van Bokhoven H, Schenck A: Epigenetic regulation of learning and memory by Drosophila EHMT/G9a. PLoS Biol. 2011, 9: e1000569-10.1371/journal.pbio.1000569. [http://dx.plos.org/10.1371/journal.pbio.1000569]PubMedPubMed CentralView ArticleGoogle Scholar
- Cooper GM, Stone Ea, Asimenos G, Green ED, Batzoglou S, Sidow A: Distribution and intensity of constraint in mammalian genomic sequence. Genome Res. 2005, 15 (7): 901-913. 10.1101/gr.3577405. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1172034&tool=pmcentrez&rendertype=abstract]PubMedPubMed CentralView ArticleGoogle Scholar
- Schmidt D, Schwalie PC, Ross-Innes CS, Hurtado A, Brown GD, Carroll JS, Flicek P, Odom DT: A CTCF-independent role for cohesin in tissue-specific transcription. Genome Res. 2010, 20 (5): 578-588. 10.1101/gr.100479.109. [http://genome.cshlp.org/cgi/content/abstract/gr.100479.109v1]PubMedPubMed CentralView ArticleGoogle Scholar
- Faure AJ, Schmidt D, Watt S, Schwalie PC, Wilson MD, Xu H, Ramsay RG, Odom DT, Flicek P: Cohesin regulates tissue-specific expression by stabilizing highly occupied cis-regulatory modules. Genome Res. 2012, 22 (11): 2163-2175. 10.1101/gr.136507.111. [http://genome.cshlp.org/content/early/2012/10/02/gr.136507.111]PubMedPubMed CentralView ArticleGoogle Scholar
- The R Development CoreTeam: R: A Language and environment for statistical computing. Vienna, Austria. 2010, [http://www.r-project.org]Google Scholar
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JGR, Korf I, Lapp H, Lehväslaiho H, Matsalla C, Mungall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E: The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002, 12 (10): 1611-1618. 10.1101/gr.361602. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=187536&tool=pmcentrez&rendertype=abstract]PubMedPubMed CentralView ArticleGoogle Scholar
- Sachs L: Angewandte Statistik (German Edition). 2003, Springer Verlag, Berlin, [http://www.amazon.com/Angewandte-Statistik-German-Lothar-Sachs/dp/3540405550]Google Scholar
- Sokal RR, Rohlf FJ: Biometry: The Principles and Practice of Statistics in Biological Research. 1995, Freeman, [http://books.google.com/books?id=N6KCNw5NHNkC&pgis=1]Google Scholar
- Grossmann S, Bauer S, Robinson PN, Vingron M: Improved detection of overrepresentation of Gene-Ontology annotations with parent child analysis. Bioinformatics (Oxford, England). 2007, 23 (22): 3024-3031. 10.1093/bioinformatics/btm440. [http://www.ncbi.nlm.nih.gov/pubmed/17848398]View ArticleGoogle Scholar