SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to deal with samples partially without replicates in DESeq? zatoichi888 Bioinformatics 2 09-12-2013 11:12 PM
Spike in dataset to compare DESeq vs cuffdiff narges Bioinformatics 1 09-07-2013 03:47 PM
Question of using DESeq for non-replicate samples camelbbs RNA Sequencing 2 11-14-2012 04:29 PM
Spike-in for identification of samples hbn Illumina/Solexa 4 08-24-2011 07:40 AM
Automated size fractionation of samples or libraries james hadfield Sample Prep / Library Generation 1 08-17-2010 05:39 AM

Reply
 
Thread Tools
Old 12-13-2013, 01:39 AM   #1
frymor
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 150
Default DESeq in libraries with spike-in samples

Hi all,

in our experiment we have control triplicates as well as two different knockouts. for each of the conditions the total RNA as well as two different fractions of polysomal (P) and non-polysomal (NP) RNAs were sequenced.

for the detection of translation efficiency the biologist added Spike-In controls into the P and NP fractions, but not to the controls.

As far as I understand it, DESeq uses the estimation of the effective library size to calculate the differentially regulated genes differentially regulated.

In the libraries with the Spike-Ins I have an enhanced library size of reads with no relation to the reference genome used in the analysis.

How do I run the DESeq analysis on these files?

I mapped the fastq files to the genomic sequence of the Spike-In and know how many reads were mapped to it.

Do I need to somehow normalize my libraries with the Spike-Ins to the number of reads mapped to this genomic sequences?

When we sequences 20 million reads, we get a lot more reads from the organism of interest, when we don't have the Spike-In in the probe. On the other hand, When I map the fastq files to the reference genome, the reads from the Spike-In controls won't map and therefore won't be in the bam files.

DESeq normalize anyhow according to the library size. Does it means, that I don't need to worry about these differences?

Thanks in advance.

Assa
frymor is offline   Reply With Quote
Old 12-16-2013, 01:16 PM   #2
jparsons
Member
 
Location: SF Bay Area

Join Date: Feb 2012
Posts: 62
Default

Assa,

You should be concerned with the effect of the differential spike-in between the treatment and control samples. This effect may be quite minor or more serious, depending on the percentage of spike to total in your sample.

Thankfully, you have the data you need. Feed DESeq the count data with the spike-in excluded (mapped only to the genomic sequence). It will handle the library size normalization correctly at this point.

To give yourself some peace of mind, you can calculate the sizeFactors both ways. You should expect the P+spike factor to be lower than the P-spike factor.
jparsons is offline   Reply With Quote
Old 12-17-2013, 04:00 AM   #3
frymor
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 150
Default

The problem I have is, that there are no Spike-In in the control samples, but only in the two treatment conditions.

I my mapping I used the genome without the Spike-In sequences. So I guess, I don't really have to woory about "contamination" of my library sizes with residues from the Spike-Ins.

Thanks for the advice

Assa
frymor is offline   Reply With Quote
Reply

Tags
deseq, spike-in

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 06:24 AM.


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