VAP: a versatile aggregate profiler for efficient genome-wide data representation and discovery

The analysis of genomic data such as ChIP-Seq usually involves representing the signal intensity level over genes or other genetic features. This is often illustrated as a curve (representing the aggregate profile of a group of genes) or as a heatmap (representing individual genes). However, no specific resource dedicated to easily generating such profiles is currently available. We therefore built the versatile aggregate profiler (VAP), designed to be used by experimental and computational biologists to generate profiles of genomic datasets over groups of regions of interest, using either an absolute or a relative method. Graphical representation of the results is automatically generated, and subgrouping can be performed easily, based on the orientation of the flanking annotations. The outputs include statistical measures to facilitate comparisons between groups or datasets. We show that, through its intuitive design and flexibility, VAP can help avoid misinterpretations of genomics data. VAP is highly efficient and designed to run on laptop computers by using a memory footprint control, but can also be easily compiled and run on servers. VAP is accessible at http://lab-jacques.recherche.usherbrooke.ca/vap/.

[1]  Bing Li,et al.  Histone H3 Lysine 36 Dimethylation (H3K36me2) Is Sufficient to Recruit the Rpd3s Histone Deacetylase Complex and to Repress Spurious Transcription* , 2009, Journal of Biological Chemistry.

[2]  Jun Song,et al.  CEAS: cis-regulatory element annotation system , 2006, Nucleic Acids Res..

[3]  J. Ahringer,et al.  Differential chromatin marking of introns and expressed exons by H3K36me3 , 2008, Nature Genetics.

[4]  Edith D. Wong,et al.  Saccharomyces Genome Database: the genomics resource of budding yeast , 2011, Nucleic Acids Res..

[5]  David Haussler,et al.  The Human Epigenome Browser at Washington University , 2011, Nature Methods.

[6]  Xin Li,et al.  Chemical-genomic dissection of the CTD code , 2010, Nature Structural &Molecular Biology.

[7]  C. Glass,et al.  Functional roles of enhancer RNAs for oestrogen-dependent transcriptional activation , 2013, Nature.

[8]  Pierre-Étienne Jacques,et al.  Yeast RNase III triggers polyadenylation-independent transcription termination. , 2009, Molecular cell.

[9]  ENCODEConsortium,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[10]  Lucas Lochovsky,et al.  ACT: aggregation and correlation toolbox for analyses of genome tracks , 2011, Bioinform..

[11]  Pierre-Étienne Jacques,et al.  The Euchromatic and Heterochromatic Landscapes Are Shaped by Antagonizing Effects of Transcription on H2A.Z Deposition , 2009, PLoS genetics.

[12]  F. Robert,et al.  Genome-wide replication-independent histone H3 exchange occurs predominantly at promoters and implicates H3 K56 acetylation and Asf1. , 2007, Molecular cell.

[13]  G. Hong,et al.  Nucleic Acids Research , 2015, Nucleic Acids Research.

[14]  Brian D. Strahl,et al.  Dimethylation of Histone H3 at Lysine 36 DemarcatesRegulatory and Nonregulatory ChromatinGenome-Wide , 2005, Molecular and Cellular Biology.

[15]  Michael R. Green,et al.  Dissecting the Regulatory Circuitry of a Eukaryotic Genome , 1998, Cell.

[16]  Matthew W Vaughn,et al.  The Fission Yeast Jmj2 Reverses Histone H3 Lysine 4 Trimethylation* , 2007, Journal of Biological Chemistry.

[17]  Mathieu Blanchette,et al.  Variant Histone H2A.Z Is Globally Localized to the Promoters of Inactive Yeast Genes and Regulates Nucleosome Positioning , 2005, PLoS biology.

[18]  Bing Li,et al.  Infrequently transcribed long genes depend on the Set2/Rpd3S pathway for accurate transcription. , 2007, Genes & development.

[19]  J. Corden,et al.  RNA polymerase II C-terminal domain: Tethering transcription to transcript and template. , 2013, Chemical reviews.

[20]  Data production leads,et al.  An integrated encyclopedia of DNA elements in the human genome , 2012 .

[21]  Dirk Eick,et al.  The RNA polymerase II carboxy-terminal domain (CTD) code. , 2013, Chemical reviews.

[22]  Raymond K. Auerbach,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[23]  Andrew B. Nobel,et al.  The Set2/Rpd3S Pathway Suppresses Cryptic Transcription without Regard to Gene Length or Transcription Frequency , 2009, PloS one.

[24]  Pierre-Étienne Jacques,et al.  A Positive Feedback Loop Links Opposing Functions of P-TEFb/Cdk9 and Histone H2B Ubiquitylation to Regulate Transcript Elongation in Fission Yeast , 2012, PLoS genetics.

[25]  Luc Gaudreau,et al.  Histone H2A.Z is essential for estrogen receptor signaling. , 2009, Genes & development.

[26]  Pierre-Étienne Jacques,et al.  DSIF and RNA Polymerase II CTD Phosphorylation Coordinate the Recruitment of Rpd3S to Actively Transcribed Genes , 2010, PLoS genetics.

[27]  F. Robert,et al.  Control of Chromatin Structure by Spt6: Different Consequences in Coding and Regulatory Regions , 2010, Molecular and Cellular Biology.

[28]  G. Ast,et al.  Chromatin organization marks exon-intron structure , 2009, Nature Structural &Molecular Biology.

[29]  Pierre-Étienne Jacques,et al.  Systematic identification of fragile sites via genome-wide location analysis of γ-H2AX , 2010, Nature Structural &Molecular Biology.

[30]  G. Bourque,et al.  The Majority of Primate-Specific Regulatory Sequences Are Derived from Transposable Elements , 2013, PLoS genetics.

[31]  F. Robert,et al.  The writers, readers, and functions of the RNA polymerase II C-terminal domain code. , 2013, Chemical reviews.

[32]  Clifford A. Meyer,et al.  Cistrome: an integrative platform for transcriptional regulation studies , 2011, Genome Biology.

[33]  Megan F. Cole,et al.  Genome-wide Map of Nucleosome Acetylation and Methylation in Yeast , 2005, Cell.

[34]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[35]  Pierre-Étienne Jacques,et al.  A universal RNA polymerase II CTD cycle is orchestrated by complex interplays between kinase, phosphatase, and isomerase enzymes along genes. , 2012, Molecular cell.

[36]  Tao Ye,et al.  seqMINER: an integrated ChIP-seq data interpretation platform , 2010, Nucleic acids research.

[37]  Johannes Söding,et al.  Uniform transitions of the general RNA polymerase II transcription complex , 2010, Nature Structural &Molecular Biology.

[38]  Mark Johnston,et al.  Linking cell cycle to histone modifications: SBF and H2B monoubiquitination machinery and cell-cycle regulation of H3K79 dimethylation. , 2009, Molecular cell.

[39]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[40]  Petra C. Schwalie,et al.  A CTCF-independent role for cohesin in tissue-specific transcription. , 2010, Genome research.

[41]  Bing Li,et al.  Combined Action of PHD and Chromo Domains Directs the Rpd3S HDAC to Transcribed Chromatin , 2007, Science.

[42]  Dindial Ramotar,et al.  The Peptidyl Prolyl Isomerase Rrd1 Regulates the Elongation of RNA Polymerase II during Transcriptional Stresses , 2011, PloS one.

[43]  Krishnamurthy Natarajan,et al.  Set2 methylation of histone H3 lysine 36 suppresses histone exchange on transcribed genes , 2012, Nature.

[44]  N. D. Clarke,et al.  Dynamic Deposition of Histone Variant H3.3 Accompanies Developmental Remodeling of the Arabidopsis Transcriptome , 2012, PLoS genetics.