Sequencing depth before and aftr normalisation with DESEQ2 doesn't look as expected: Problem?
1
0
Entering edit mode
xavier • 0
@54f3711c
Last seen 3.2 years ago
France

Hi, I'm still new to DESEQ2 and I wonder if I used everything correctly. I have several RNA-seq (stranded, PE, 30 million reads, 27 million aligned, mouse). I'm trying to check what happens after normalisation and you can see below that there are some differences in sequencing depth that are partly corrected after normalisation. However the last 6 samples don't align with the others. Am I looking at the right parameter to see if it was correctly normalized? If yes, does that mean normalisation is not good enough? why and can I improve it? thanks for your help,

I used the code below after DESEQ2 to generate the graphs (the not normalized just does not have normalized true


df2=data.frame(Tissu=colData$condition, Individu=rownames(colData),
               prof=colSums(counts(dds, normalized = TRUE)))
# Dataframe with depth (sum of reads) 
# after normalisation.
png("Profondeur_apres_normalisation_tripl.png", 1000, 1000, pointsize=20) 
ggplot(df2, aes(x=factor(Individu), y=prof, fill =Tissu))+
  geom_bar(stat ="identity") + xlab("Condition") + ylab("Profondeur")+ 
  theme(axis.text.x =element_text(angle = 90))
dev.off()

enter image description here

normalisation DESEQ2 • 1.2k views
ADD COMMENT
1
Entering edit mode

colSums is (imho) only partially informative towards how "well" normalization works. Exactly-aligned colSums are only expected if you normalize naively with reads per million as this will make colSums more or less equal. As DESeq2 (and edgeR) additionally correct for composition your plot could simply indicate that composition is different in these last samples. If you want to really "see" the normalization then make some MA-plots. The bulk of data points should be centered somewhat at y=0. The easiest would be to simply run the DESeq2 pipeline as normal and then use plotMA on the results object.

ADD REPLY
2
Entering edit mode
swbarnes2 ★ 1.4k
@swbarnes2-14086
Last seen 3 hours ago
San Diego

I'd try looking at the PCA, see if that makes sense, instead of columsums. But yes, MA plots will have a tilt to them if your size normalization is off, since it will introduce a fold change bias.

ADD COMMENT

Login before adding your answer.

Traffic: 703 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