chip-seq analysis and peak annotation
2
0
Entering edit mode
4.2 years ago
tarek.mohamed ▴ 280

Hi All, I am analyzing chip-seq data. I did the alignment using bowtie using hg38 as a reference genome, and peaks calling using macs2. Now I am using R package "chipseeker" to annotate these peaks which needs a txdb object.

My question is can I use txdb of hg19 to annotate peaks, however I used hg38 genome as reference during alignment.

Thanks, Tarek

ChIP-Seq • 2.7k views
ADD COMMENT
1
Entering edit mode
4.2 years ago
EagleEye 7.0k

Try homer annotate peak to annotate your peaks with hg38 annotation (GTF). Remember you should always maintain the same genome/gene/transcript annotation throughout your analysis.

ADD COMMENT
0
Entering edit mode
4.2 years ago
Sinji ★ 3.1k

As another option, you can use ChIPseeker to annotate to the hg38 genome in the exact same way you are probably using hg19 (I assume you are following the Bioconductor guide). All you need to do is download and install this reference package and use it as your txdb object.

If you did not use a UCSC reference genome, then use a GTF file to generate your own txdb annotation using the GenomicFeatures R package.

ADD COMMENT

Login before adding your answer.

Traffic: 1987 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