This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RNA-seq dedupe PCR contamination before or after mapping

Will deduping (marking duplicates) with Picard before mapping affect my variant calling? Is it common practice to dedupe after mapping?

Reason why I ask is because my FastQC reports still have a lot of Kmer, overrepresented sequences, and bad GC content. I figured these can be corrected by removing PCR contamination. This is after trimming adapter and low quality (10) bases using BBDuk.

rna-seq
• 4,361 views • link
updated by Carlo Yague Carlo Yague • written by umn_bist • umn_bist

Depends on what data you have, but a slight bimodal distribution of GC content in whole exome data, seems to be the norm (I haven't figured out a reason why, but it appears to be commonplace)

• 0 views • link
written by andrew.j.skelton73 andrew.j.skelton73

1 answer

Overrepresented sequences / skewed GC content is expected in RNA-seq data. It usually comes from the most highly expressed transcripts (such as rRNA). However, it can also come from PCR duplicates and those can completely skew variant calling. For this reason, while people usually don't dedupe RNA-seq data for differential expression analysis, it is still recommended to do so for variant calling.

Some reference : http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0058815

• 0 views • link
written by Carlo Yague Carlo Yague

Log in to answer this question.

More posts like this