This is a test version of Biostars. For the public version, visit https://www.biostars.org.
rRNA detection (for contamination) in RNA-seq

Hi all,

I am working with human RNA-seq data.I want to know whats the best way to get information on rRNA reads in RNA-seq? Does rRNA gets passed on in the bam files after alignment(using human genome) or we have to check the raw fastq reads only for this?

Also I want to look for this irrespective of the library they came from (polyA or rRNA removal kits),as I just want to have an estimate of this rRNA reads as a QC measure.

let me know your thoughts.

Thanks,

Ron

rna-seq next-gen alignment star
• 13,245 views • link
updated by Charles Plessy Charles Plessy • written by Ron Ron

5 answers

I used RseQC for rRNA detection. They have a bed file too.

http://rseqc.sourceforge.net/

• 0 views • link
written by Ron Ron

The rRNA-cluster are distributed in repeats over the genome. Therefore, reads origin from rRNA-molecules map too many times and get discarded by most RNA-SEQ aligners.

I'd download the rRNA fasta files from NCBI and map all reads against them (e.g. BBmap, bowtie2); the mapping rate is equal to your rRNA contamination. Additionally, you can use this as a first step to sort these reads before continuing your pipeline.

• 0 views • link
written by michael.ante michael.ante

I asked a related question recently that might be useful: RNA-seq rRNA contamination

• 1 views • link
written by igor igor

A very simple approach would be to download rRNA bed coordinates (from Ensembl) and then count the total number of reads that fall into this ROI (using bedtools suite) and divide this by total number of reads in the bam (samtools idxstats). This would give you a rough estimate of your background rRNA levels.

• 0 views • link
written by nick.a.rouse • nick.a.rouse

You can count or filter out rRNA reads from FASTQ files (single- or paired-end) using TagDust (and do many other things).

• 0 views • link
written by Charles Plessy Charles Plessy

Log in to answer this question.

More posts like this