Question: ribosomal RNA contamination
1
gravatar for lkmklsmn
2.7 years ago by
lkmklsmn890
United States
lkmklsmn890 wrote:

Hi,
I am working on a RNAseq data set of about 20 human tissue samples. However, some of the samples have a very low alignment rate (~30-40%) while the rest has a a decent alignment rate around 80%. I aligned the reads using STAR with default settings. The log file tells me that for those samples with low alignment rate many reads did not map due to multimapping. My first guess is that there could have been an issue with the rRNA depletion for those samples. How can I show that the rRNA depletion failed?
Take unmapped reads and align to rRNA sequences? Has anybody had similar issues before?
Any feedback is greatly appreciated!

contamination rrna rnaseq • 1.3k views
ADD COMMENTlink modified 2.7 years ago by Devon Ryan91k • written 2.7 years ago by lkmklsmn890
2
gravatar for Devon Ryan
2.7 years ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k wrote:

Yes, you can align against the 45S sequence, noting that this will (in my experience) somewhat underestimate the actual amount of rRNA in the sample.

As an aside, I should note that rRNA depletion only seems to work well if you do it on fresh samples. If you freeze them beforehand then the % of rRNA remaining is more likely to be high (this is purely observational).

ADD COMMENTlink written 2.7 years ago by Devon Ryan91k

Thanks! Its true that some of the samples had low RINs. Is there any way to correct for this? And still use the samples. As of right now, PCA clustering does not reveal my treatment structure and I am worried these samples are unusable.

ADD REPLYlink written 2.7 years ago by lkmklsmn890

Low RINs tend to correlate with degradation, so you might give salmon a try with whatever the position bias correction option is. Perhaps that'd help.

ADD REPLYlink written 2.7 years ago by Devon Ryan91k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 918 users visited in the last hour