Home

metagene: A package to produce Metafeature plots

This repository contains R functions used to do multiple ChIP-Seq experiments comparisons.

This package produces Metagene-like plots to compare the behavior of DNA-interacting proteins at selected groups of features. A typical analysis can be done in viscinity of transcription start sites (TSS) of genes or at any regions of interest (such as enhancers). Multiple combinations of group of features and/or group of bam files can be compared in a single analysis. Bootstraping analysis is used to compare the groups and locate regions with statistically different enrichment profiles. In order to increase the sensitivity of the analysis, alignment data is used instead of peaks produced with peak callers (i.e.: MACS2 or PICS). The metagene package uses bootstrap to obtain a better estimation of the mean enrichment and the confidence interval for every group of samples.

Examples using real biological data

Compare the enrichment of the H3K27ac histone mark in active and bivalent transcription start sites (TSS)

Compare the enrichment of the H3K27ac and H4K27me3 histone marks in active transcription start sites (TSS)

How to do it…

Extract ChIP-Seq profiles values from metagene output

Authors

Charles Joly Beauparlant, Fabien Claude Lamaze, Rawane Samb, Astrid Louise Deschenes and [Arnaud Droit] (http://ca.linkedin.com/in/drarnaud "Arnaud Droit").

See Arnaud Droit Lab website.

License

This package and the underlying metagene code are distributed under the Artistic license 2.0. You are free to use and redistribute this software.

For more information on Artistic 2.0 License see http://opensource.org/licenses/Artistic-2.0

Bugs/Feature requests

If you have any bugs or feature requests, let me know. Thanks!