Question: a problem with peak annotation after ChIP-seq
0
gravatar for gpdus67
5.1 years ago by
gpdus670
Korea, Republic Of
gpdus670 wrote:

Hi. I downloaded GEO data(GSE22276) and try doing downstream analysis. By the way I got a problem in the step of annotation with HOMER toolannotatePeaks.pl). Input data is the output of previous step(in HOMER,it's findPeaks.pl) in the format of bed.

but result came out with anntotion column blank, NA and same for  following colums.

     

input file has chromosome(chr1), start and end position in first three positions, which is necessary, I have no idea about why annotation failed.

any idea for this?

 

Best wishes.

chip-seq • 3.3k views
ADD COMMENTlink modified 3.4 years ago by dcwest0 • written 5.1 years ago by gpdus670
2
gravatar for geek_y
5.1 years ago by
geek_y10k
Barcelona
geek_y10k wrote:

toolannotatePeaks.pl takes a BED file and a genome file. Are you providing the proper genome file/name ? Check if the chromosome names are same ? Like 'Chr1' or '1' in both Genome and BED file ? Post the exact command used and first few lines of BED file.

ADD COMMENTlink modified 5.1 years ago • written 5.1 years ago by geek_y10k

I solved the problem. in bed file, chromosome names are like 'chr1' while '1' in genome and gtf files. i made them same and it works.  I really appreciate your help. thanks.

ADD REPLYlink written 5.1 years ago by gpdus670
0
gravatar for Rituriya
3.4 years ago by
Rituriya30
Canada
Rituriya30 wrote:

Hi All,

I'm facing similar issue from HOMER tool, reported by gpdus67. I have checked my peaks file and the genome file both have 'chr1' only. Still I always end up getting blank values in Annotated peak file. Is there anything to do with makeTagDirectory output? Command used:

./annotatePeaks.pl path-to-dir/peaks.txt hg18 > annotatedpeaks.txt

What am I doing wrong?

Peaks file also contains a lot of commented information in the beginning. Does it needs to be removed? Few lines from my peaks file:

#PeakID chr     start   end     strand  Normalized Tag Count    focus ratio     findPeaks Score Fold Change vs Local    p-value vs Local        Clonal Fold Change
chr1-1  chr1    91625448        91625674        +       258.5   0.584   378.000000      84.31   0.00e+00        0.49
chr1-2  chr1    154452940       154453166       +       255.1   0.597   372.000000      54.33   0.00e+00        0.49
chr1-4  chr1    121186277       121186503       +       214.0   0.561   314.000000      13.15   1.48e-225       0.54

I am tired of getting empty annotated peak file. Please help.

Thanks in advance.

ADD COMMENTlink written 3.4 years ago by Rituriya30

It must be removed, HOMER is peculiar about this.

ADD REPLYlink written 3.4 years ago by dally190
0
gravatar for Guangchuang Yu
3.4 years ago by
Guangchuang Yu2.3k
China/Guangzhou/Southern Medical University
Guangchuang Yu2.3k wrote:

you can try ChIPseeker.

ADD COMMENTlink written 3.4 years ago by Guangchuang Yu2.3k
0
gravatar for dcwest
3.4 years ago by
dcwest0
dcwest0 wrote:

Hi gpdus67 and Rituriya, I'm not sure if this will help, but I had the same issue too--- I ended up needing to convert my .txt peak file from a Macs-based .txt file into an acceptable .txt file that annotate.pl could handle. So, in Homer, I did:

changeNewLine.pl path/peakfile.txt

I did this for all of my peak files... then it was able to annotate my stuff. Hope this helps! ~Diana

ADD COMMENTlink written 3.4 years ago by dcwest0
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: 1401 users visited in the last hour