How to assess biases in ribosomally depleted RNAseq data?
0
1
Entering edit mode
4.8 years ago
Lina F ▴ 200

Hi all,

I am interested in evaluating four different ribosomal depletion wet lab protocols and I want to see if any of them are biased against transcripts other than rRNA. Ultimately, I want to find out which of the protocols is the least biased towards other transcripts.

I quantified all the samples with sailfish and now I have number of reads and TPM for each transcript. For each of the four different wet lab protocols I have between 6 and 30 samples.

Is it enough to just compare TPM counts for each transcript and maybe look at TPM variances, or should I do a regular DESeq differential expression analysis?

Also, unfortunately there are some samples that were only evaluated in one protocol (sorry, that's what the lab gave me...). Is it worth to even include these? If yes, should they be treated differently?

Thanks for any suggestions!

Edit: I ran SortMeRNA to remove ribosomal RNA reads before quantifying with sailfish. This is ecoli data, in case it matters.

RNA-Seq bias transcripts • 1.2k views
ADD COMMENT
0
Entering edit mode

Have you aligned the samples against rDNA repeat for your organism to estimate what fraction of reads are mapping to it?

ADD REPLY
0
Entering edit mode

Yes I have. I ran SortMeRNA to remove ribosomal RNA reads before quantifying with sailfish. This is ecoli data, in case it matters.

ADD REPLY
0
Entering edit mode

If you look at the number of rRNA reads that got removed from various samples do you see any correlation with the TPM counts you have? Ultimately getting the right DE result is the ultimate litmus test. Individual samples (not tested against the entire spectrum) won't be informative to make a decision across the tests, won't they?

ADD REPLY
0
Entering edit mode

Could you elaborate on how to correlate the data?

I plotted three samples below that were treated with three different rRNA depletion protocols. The third protocol removed the most rRNA (i.e. the percentage value on the 2nd y-axis is highest).

enter image description here

enter image description here

enter image description here

ADD REPLY

Login before adding your answer.

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