Will large differences in the number of reads between samples decrease reliability of DESeq2 analysis?
1
0
Entering edit mode
Mike ▴ 10
@mike-12142
Last seen 2.6 years ago
Canada

I have some RNA-seq samples from mouse, 2 conditions with 4 replicates each, read quality is good and for each sample 85-90% of reads align. The number of aligned reads in millions are:

Condition 1 replicate 1: 100
Condition 1 replicate 2: 79
Condition 1 replicate 3: 52
Condition 1 replicate 4: 37
Condition 2 replicate 1: 56
Condition 2 replicate 2: 59
Condition 2 replicate 3: 31
Condition 2 replicate 4: 26

I know DESeq2 normalizes based on the number of reads but I'm wondering if the variability here will cause any problems?

deseq2 • 910 views
ADD COMMENT
0
Entering edit mode
@ryan-c-thompson-5618
Last seen 8 months ago
Scripps Research, La Jolla, CA

This looks like a reasonable, expected variation in total read counts. All the totals are within one order of magnitude. I don't think the differences in sequencing depth will be a problem for this data. I've used similar methods successfully on data sets with much wider ranges of total counts.

ADD COMMENT
0
Entering edit mode

Thanks that's good to know.

>This looks like a reasonable, expected variation in total read counts.

We sequenced on a NextSeq and expected the reads evenly distributed between the 8 samples and I was surprised by the variability, I expected some but not this much. I was planning to follow up with the sequencing facility because I was a bit concerned that the variability could be due to errors in quantification or poor quality control by them (we submit RNA and they do the rest). But you think the numbers seen here are normal?

ADD REPLY

Login before adding your answer.

Traffic: 853 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6