This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Metaseq: Average ChIP-seq signal over promoters error

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
• 4,190 views • link
updated by Ram Ram • written by Sudhir Jadhao Sudhir Jadhao

Dear all,

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

• 0 views • link
updated by Ram Ram • written by hlsz.laszlo • hlsz.laszlo

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.

• 0 views • link
updated by Ram Ram • written by Fidel Fidel

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

• 0 views • link
written by endrebak endrebak

3 answers

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:

• 0 views • link
updated by Ram Ram • written by Ryan Dale Ryan Dale

See my post here.

• 0 views • link
updated by Ram Ram • written by Ming Tommy Tang Ming Tommy Tang

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
• 0 views • link
updated by Ram Ram • written by Sudhir Jadhao Sudhir Jadhao

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

• 0 views • link
updated by Ram Ram • written by chemcehn chemcehn

Log in to answer this question.

More posts like this