RNAseq uneven PCR bias "how to take into best take into account uneven PCR bias ?"
0
0
Entering edit mode
7.2 years ago

Hello, I am a student doing a project with some conceptual difficulties (I do realise such question optimally shouldn't occur). The project aims to compare 8 RNA-seq samples (4 vs 4) testing for differentially expressed genes. The problem is that 2 of the samples (1 for each categories) weren’t amplified as much as the others (an additional 5 cycles).

Obviously, such treatment is sub-optimal for comparative studies, I was wondering what was the best remedy for this problem short of resequencing. I was considering DESEQ2 for it’s outlier capacities or dropping the samples (doing a 3 vs 3) Additional information: the transcriptome is de novo using the samples.

RNA-Seq rna-seq • 1.8k views
ADD COMMENT
1
Entering edit mode

I think that if you remove PCR duplicates and use DESeq2 you shouldn't have a bias

ADD REPLY
3
Entering edit mode

The user then runs the risk of removing actual signal, If you can unequivocally see on a PCA that those two samples are very different from their respective groups, I'd maybe include it as a batch effect in the DESeq2 model

ADD REPLY
2
Entering edit mode

Of course you would have a bias! Removing duplicates would incorrectly decrease the inferred expression level of highly expressed genes, from which duplicates are normal and expected. That is why rmdup is not recommended for RNA-seq (or any NGS assay in which duplicate generation is expected irrespective of PCR duplicates).

ADD REPLY
0
Entering edit mode

but is there in your opinion a correct solution ? Is the batch effect integration an acceptable way in your opinion ?

ADD REPLY
1
Entering edit mode

It's a good solution

ADD REPLY

Login before adding your answer.

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