creating custom genome annotation for the ChIPQC package
1
1
Entering edit mode
@philip-lijnzaad-2499
Last seen 2.4 years ago
European Union

This is a summary of an e-mail exchange I had with Tom Caroll, about creating custom annotations for the ChIPQC package. Reproduced so that other can make use of it.

Philip Lijnzaad

--------

> Dear Tom and/or Rory,
> I am interested to use your ChIPQC package for our yeast data, but there is no standard yeast annotation that I can use.

On 08/07/2015 11:14 AM, tom carroll wrote:

You can provide your own annotation in the form of a named list of GRanges with the first element of the list containing a character vector of version.
I include an example below.  I did do a course in Boston you may find useful and contains information on custom annotation.http://bioconductor.org/help/course-materials/2014/BioC2014/Bioc2014_ChIPQC_Practical.pdf

##########################################
## You can also provide your own denovo List of GRanges but first
element must be a character with name "version".

library(GenomicRanges)
library(ChIPQC)
customAnnotation <- list(
  version="TomsCustom",
  Test=GRanges("chr19",IRanges(1,10000000)),
  AnotherTest=GRanges("chr19",IRanges(10000000,20000000))
  )
testExp <- ChIPQCsample("ETOH_Rep2.bamRange5.bam",annotation=customAnnotation,chromosomes=c("chr19"))
plotRegi(testExp)
##########

 

ChIPSeq QC ChIPQC • 3.5k views
ADD COMMENT
2
Entering edit mode
@thomas-carroll-7019
Last seen 2.1 years ago
United States/New York/The Rockefeller …

hi Philip,

Thanks for posting this here.

You can provide your own annotation in the form of a named list of GRanges with the first element of the list containing a character vector of version.

I include an example below.  I did a course in Boston you may find useful and contains information on custom annotation.

http://bioconductor.org/help/course-materials/2014/BioC2014/Bioc2014_ChIPQC_Practical.pdf

I can look to add a yeast annotation too. As long as we have a suitable TXDB package it should be simple although i am sure one can create our own TXDB for any genome not currently supported.

best,

Tom

 

 

 


##########################################
## You can also provide your own denovo List of GRanges but first
element must be a character with name "version".

library(GenomicRanges)

library(ChIPQC)
customAnnotation <- list(
  version="TomsCustom",
  Test=GRanges("chr19",IRanges(1,10000000)),
  AnotherTest=GRanges("chr19",IRanges(10000000,20000000))
  )
testExp <- ChIPQCsample("ETOH_Rep2.bamRange5.bam",annotation=customAnnotation,chromosomes=c("chr19"))
plotRegi(testExp)
##########

ADD COMMENT

Login before adding your answer.

Traffic: 446 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6