This is a test version of Biostars. For the public version, visit https://www.biostars.org.
rRNA removal needed for salmon / DESeq2 analysis

Hi all,

I am re-analysing some RNA-seq data that I'm told have high levels of rRNA contamination (not sure how high yet). My usual workflow is salmon -> tximport -> DESeq2

I've seen recommendations that any reads matching rRNA need to removed prior to analysis. Is this still the case with the salmon (alignment-free) protocol? if the rRNA are not in my gtf file, then they will not be quantified and appear in the analysis. Obviously I'll get smaller counts than expected, but should be Ok to proceed with the counts for all the other genes?

Thanks,

rrna contamination salmon deseq2
• 4,348 views • link
updated by ATpoint ATpoint • written by mark.dunning mark.dunning

So the transcript file I was planning to use from Gencode does include rRNA sequences that I will have to remove before calling salmon

• 0 views • link
written by mark.dunning mark.dunning

1 answer

Unless rRNA is present in your transcriptome fasta file there's no reason to remove it (likewise, there's no reason to remove it in other cases unless it's actually annotated). What you likely should do is quantify how variable rRNA levels are. Sometimes rRNA level ends up being a nice surrogate for other issues (e.g., RNA degradation) with a given sample and you might need to use it for reliably excluding samples.

• 0 views • link
written by Devon Ryan Devon Ryan

Thanks! That seems like sensible advice.

• 0 views • link
written by mark.dunning mark.dunning

Are there any papers that study these issues (rRNA level as surrogate for RNA degradation, etc.)?

• 0 views • link
written by grp2009 • grp2009

Probably, search pubmed for them.

• 0 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this