SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq simulation Zimbobo Bioinformatics 17 08-06-2014 06:30 AM
ChIP-Seq: ChIP-Array: combinatory analysis of ChIP-seq/chip and microarray gene expre Newsbot! Literature Watch 0 05-19-2011 03:50 AM
ChIP-Seq: ChIP-chip versus ChIP-seq: Lessons for experimental design and data analysi Newsbot! Literature Watch 0 03-02-2011 03:50 AM
ChIP-Seq: ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip dat Newsbot! Literature Watch 0 05-13-2010 03:00 AM
PubMed: An integrated software system for analyzing ChIP-chip and ChIP-seq data. Newsbot! Literature Watch 0 11-04-2008 06:03 AM

Reply
 
Thread Tools
Old 08-06-2010, 07:44 PM   #1
golharam
Member
 
Location: Philadelphia, PA

Join Date: Dec 2009
Posts: 55
Default ChIP-Seq Simulation

Does anyone simulate a ChIP-Seq experiment prior to actually running the experiment? This is coming up because our collaborators want to simulate the experiment to optimize the parameters (aka read length, depth of coverage, etc) to make sure they are doing things correctly prior to sequencing. It makes sense considering we have several options on read length, among other things. So my question goes out to everyone, do you perform a simulation on ChIP-Seq experiments prior to running the experiment? We would like to do so for the following reasons: We are performing this experiment on an organism with a genome of about 35 MB, and the signal in our initial experiment is low, but we have a LOT of coverage, probably more than necessary.

1) What is the optimal read length to sequence to. Is 35, 50, 75, or 100bp enough? If 50 or 75 is enough to uniquely place 90% of the reads and 100bp doesn't give us any added benefit, then we would only like to sequence up to 75bp and save money.

2) How much starting material do we need. Is this really something we can answer prior to sequencing? We have enough based on our core facilities recommendations. Is that good enough?

3) How much coverage are we going to need? We expect this would dictate how much starting material we would need, and how many lanes to run.

Is there anything else or is all of this just overkill? I've read 3 or 4 papers regarding modeling ChIP-Seq in silico, but understanding the underlying code is proving to be difficult, and I'm not convinced its right for us.

I just want to make we are on the right track and would appreciate to hear what others are doing.
golharam is offline   Reply With Quote
Old 08-06-2010, 11:40 PM   #2
dawe
Senior Member
 
Location: 4530'25.22"N / 915'53.00"E

Join Date: Apr 2009
Posts: 258
Default

Hi there!

Quote:
Originally Posted by golharam View Post
Does anyone simulate a ChIP-Seq experiment prior to actually running the experiment?
AFAIK no, but if it was possible I would like to know :-)

Quote:
Originally Posted by golharam View Post
1) What is the optimal read length to sequence to. Is 35, 50, 75, or 100bp enough? If 50 or 75 is enough to uniquely place 90% of the reads and 100bp doesn't give us any added benefit, then we would only like to sequence up to 75bp and save money.
We currently use 36 bp in single read for human chip-seq and we are happy with it. The IP material is a small fraction of the entire genome and it is highly covered with short reads. You have a 35 Mb genome, which can be covered ~20x with a single average Illumina lane (36bp, SR), your IP will probably be much less. You should be able to have a good IP profile.

Quote:
Originally Posted by golharam View Post
2) How much starting material do we need. Is this really something we can answer prior to sequencing? We have enough based on our core facilities recommendations. Is that good enough?
I believe you should refer to your machine standards.

Quote:
Originally Posted by golharam View Post
3) How much coverage are we going to need? We expect this would dictate how much starting material we would need, and how many lanes to run.
As said before, I don't think you'll have any coverage issue here.

Quote:
Originally Posted by golharam View Post
Is there anything else or is all of this just overkill? I've read 3 or 4 papers regarding modeling ChIP-Seq in silico, but understanding the underlying code is proving to be difficult, and I'm not convinced its right for us.

I just want to make we are on the right track and would appreciate to hear what others are doing.
Can you post those papers?
BTW, the only problem I can imagine is related to a "low" signal (as you said). If your antibody is not specific, well, it may be hard to distinguish the IP from the Input background (needless to ask:
will you run an Input?). If you don't know about Ab performance, it may be useful to run a third sample with an aspecific Ab (such as IgG or something).

d
dawe is offline   Reply With Quote
Old 08-07-2010, 06:28 AM   #3
Nix
Member
 
Location: SLC, Utah

Join Date: Jun 2008
Posts: 60
Default

Hmm,

36bp is sufficient for chIP-seq, shoot for >10-20 million duplicate free uniquely mapped reads (1/2 - 1 lane on a GAIIx), and at minimum provide 5ng of correctly sized chIP product to your sequencing core. Lastly, run 3-4 biological replicas for each sample. Barcoding will cut the cost. Very shortly, reviewers are going to reject papers that lack these controls. See http://useq.sourceforge.net/usage.html for additional details.
Nix is offline   Reply With Quote
Old 08-07-2010, 12:56 PM   #4
golharam
Member
 
Location: Philadelphia, PA

Join Date: Dec 2009
Posts: 55
Default

PeakSeq enables systematic scoring of ChIP-seq experiments relative to controls
Nature Biotechnology, Vol. 27, No. 1. (04 January 2009), pp. 66-75.

Modeling ChIP Sequencing In Silico with Applications
PLoS Computational Biology, Vol. 4, No. 8. (22 August 2008), e1000158.

Design and analysis of ChIP-seq experiments for DNA-binding proteins
Nature Biotechnology, Vol. 26, No. 12. (16 November 2008), pp. 1351-1359.
golharam is offline   Reply With Quote
Reply

Tags
chip-seq

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 03:17 AM.


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