Metaseq: Average ChIP-seq signal over promoters error
3
0
Entering edit mode
9.1 years ago

Hello everyone,

I am new to python,

I am trying to plot chipseq signal over promoter using metaseq . The script available at metaseq site, It uses gtf file in the from of database. I don't have gtf file in the from of database, its giving error for that . Anyone have metaseq script for plotting chipseq signal over promoter which uses bam file data around 1.5GB (use the user defined data, not example data) and .gtf file

ChIP-Seq • 3.1k views
ADD COMMENT
1
Entering edit mode

Hi, I recommend you to use deepTools. You just need three commands and a bed file that you can download from UCSC:

$ bamCoverage -b file.bam -o file.bigwig
$ computeMatrix reference-point --scoreFile file.bigwig --regionsFile genes.bed -- beforeRegionStartLength 2000 --afterRegionStartLength 500 --referencePoint TSS -out matrix.tab.gz
$ profiler -m matrix.tab.gz -o profile.png

However, is not recommended to directly plot chip-seq signals as biases are frequently found. E.g. human and mouse promoters tend to have an enrichment of reads just because they are GC rich. Instead the log2ratio of ChIP-seq vs. input is preferable.

ADD REPLY
0
Entering edit mode

deeptools seems great. Note: in the code above profiler has become plotProfile

ADD REPLY
0
Entering edit mode

Dear all,

How to create a meta-gene profile using metaseq? Unfortunately, I didn't find any help for this.

ADD REPLY
2
Entering edit mode
9.1 years ago
Ryan Dale 5.0k

edit: For more standard cases like what you're asking about, try deepTools, ngsplot, or HOMER for a much simpler approach. Use metaseq when you need a lot of customization or are integrating with other tools.

metaseq's example 1 uses a database to show how to get TSSs. You can see this question and others for how to get TSSs or promoters using other methods, and just substitute that file for tsses.gtf in the example.

For a simpler example, see Plotting Reads Around Tss, which refers to the following script that you can modify for your data. Make sure to change the "processes" argument if you have multiple cores on your computer:

ADD COMMENT
0
Entering edit mode
9.1 years ago
Ming Tommy Tang ★ 3.8k

See my post here.

ADD COMMENT
0
Entering edit mode

Thanks for reply..

I tried with HTSeq but its giving below error

Traceback (most recent call last):
  File "streaming_throughRead.py", line 17, in <module>
    tsspos[window]+=p
  File "_HTSeq.pyx", line 526, in HTSeq._HTSeq.GenomicArray.__getitem__ (src/_HTSeq.c:9489)
  File "_HTSeq.pyx", line 372, in HTSeq._HTSeq.ChromVector.__getitem__ (src/_HTSeq.c:6649)
IndexError: start too small
ADD REPLY
0
Entering edit mode
9.1 years ago
chemcehn ▴ 210

Just try another program, e.g. HOMER, it is much simpler to use

ADD COMMENT

Login before adding your answer.

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