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

Thread Tools
Old 12-13-2013, 01:39 AM   #1
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.

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

Join Date: Feb 2012
Posts: 62


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
Senior Member
Location: Germany

Join Date: May 2010
Posts: 150

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

frymor is offline   Reply With Quote

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