Question: chip-seq analysis and peak annotation
0
gravatar for tarek.mohamed
2.8 years ago by
tarek.mohamed260
tarek.mohamed260 wrote:

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 • 1.9k views
ADD COMMENTlink modified 2.8 years ago by Sinji3.0k • written 2.8 years ago by tarek.mohamed260
1
gravatar for EagleEye
2.8 years ago by
EagleEye6.6k
Sweden
EagleEye6.6k wrote:

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 COMMENTlink modified 2.8 years ago • written 2.8 years ago by EagleEye6.6k
0
gravatar for Sinji
2.8 years ago by
Sinji3.0k
UT Southwestern Medical Center
Sinji3.0k wrote:

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 COMMENTlink modified 2.8 years ago • written 2.8 years ago by Sinji3.0k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2190 users visited in the last hour