SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
any software to draw gene track together with ChIP-seq peak crazyhottommy Bioinformatics 3 06-07-2013 01:10 PM
average profiles over gene bodies bogdan Bioinformatics 0 06-06-2012 01:59 PM
ChIP-Seq: A new algorithm for quantifying binding site pattern similarity with applic Newsbot! Literature Watch 0 03-28-2012 06:30 AM

Reply
 
Thread Tools
Old 11-05-2013, 07:23 PM   #1
polvo
Junior Member
 
Location: LA

Join Date: Nov 2013
Posts: 3
Default Quantifying ChIP-seq peak density in gene bodies

I am looking for a way to quantify all of the peak density/abundance of ChIP signal for a POL II ChIP-seq experiment in every refseq gene body that doesn't involve normalization for gene size and can be a variable span accommodate each gene from TSS to TTS. The only software I know that quantifies abundance requires you to either enter a fixed span for each genomic region or normalizes for gene length. Thank you.
polvo is offline   Reply With Quote
Old 11-06-2013, 06:23 AM   #2
crazyhottommy
Senior Member
 
Location: Gainesville

Join Date: Apr 2012
Posts: 140
Default

Have a look at the Homer software:
http://biowhat.ucsd.edu/homer/ngs/analyzeRNA.html

"Measuring Gene Expression in Exons vs. Gene Bodies.
Depending on the type of sequencing you are analyzing, you will want to quantify RNA from different parts of the gene. The "-count [...]" option controls which regions of the gene are used for analysis (use like "-count exons" or "-count genes"). The options below only pertain to 'rna' or a custom GTF file:
exons - Counts tags in exons only. Use this for most applications of RNA-Seq, such as polyA-RNA-seq or other techniques that aim to measure mRNA.
cds - Counts tags in coding regions only. This could be useful for quantifying ribosome coverage on coding sequences with techniques such as Ribo-Seq
introns - Counts tags on introns only.
5utr or 3utr - Count tags on 5' UTR and 3' UTR regions, respectively.
genes (default) - Counts tags on the full gene body (TSS to TTS). This is useful for GRO-Seq where we expect coverage across the entire transcript. Can also be used to quantify H3K36me3 or PolII ChIP-Seq."
crazyhottommy is offline   Reply With Quote
Old 11-06-2013, 06:28 AM   #3
crazyhottommy
Senior Member
 
Location: Gainesville

Join Date: Apr 2012
Posts: 140
Default

Quote:
Originally Posted by polvo View Post
I am looking for a way to quantify all of the peak density/abundance of ChIP signal for a POL II ChIP-seq experiment in every refseq gene body that doesn't involve normalization for gene size and can be a variable span accommodate each gene from TSS to TTS. The only software I know that quantifies abundance requires you to either enter a fixed span for each genomic region or normalizes for gene length. Thank you.
Or you may provide a bed file containing all the gene information from TSS to TSS
and use coverageBed from bedtools http://bedtools.readthedocs.org/en/l.../coverage.html
crazyhottommy is offline   Reply With Quote
Old 11-06-2013, 03:13 PM   #4
jwfoley
Senior Member
 
Location: Stanford

Join Date: Jun 2009
Posts: 179
Default

You could use UniPeak and annotate with known genes, as done in this paper: http://www.biomedcentral.com/1471-2164/14/720/abstract

But you're unlikely to find pol II very far into gene bodies anyway, unless you used the antibody for the serine 2-phosphorylated form. So if you start from entire gene-body annotations you may have problems.

Last edited by jwfoley; 11-06-2013 at 03:24 PM.
jwfoley is offline   Reply With Quote
Reply

Tags
chip-seq quantification

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:57 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO