Should I remove host RNA from host-pathogen metatranscriptomic data
1
1
Entering edit mode
9.9 years ago
biotech ▴ 570

Does removing host RNA before mapping reads to pathogen genome adds any benefits?

I'm worried about the false positive host RNA could give.

I'm working with pig and bacterial data.

Thanks!

metatranscriptomics • 3.3k views
ADD COMMENT
2
Entering edit mode
9.9 years ago

In my experience, one does not benefit that much from such an approach in this case (mammal vs bacterium). Besides, wouldn't you still have the false positive problem if you pre-mapped to pig? (i e wrongly assigning bacterial reads to pig) Anyway, I think there is enough divergence that you would not get many false positive hits.

I did such a project (mouse + a certain bacterium) and there the problem was rather that we had a bacterial genome that we mapped to and got a lot of hits to, but in the end it turned out that most of those were from other bacteria (mostly E coli) that were contaminating the sample. A useful way to spot such erroneous cross-bacterial mappings turned out to be looking at the "mismatch spectrum" I e how many of your reads have 0, 1, 2 ... mismatches in the alignment to the genome you are mapping against (this is easy to get from the BAM/SAM). It should be mostly zeroes and ones, but in our case the most common number of mismatches was 6 or so. So the bacteria were close enough that reads from one could be mapped to the other, but always with many mismatches. Perhaps that could be helpful for you as well.

ADD COMMENT
0
Entering edit mode

Thanks Mikael, some references would be also appreciated.

ADD REPLY
0
Entering edit mode

This is the study I mentioned:

http://journals.plos.org/plospathogens/article?id=10.1371/journal.ppat.1004600

I'm afraid I don't have anything else I can think of right now.

ADD REPLY
0
Entering edit mode

A great reference Mikael, although I'm surprised how you were able to publish this in vivo data. Extremely low abundance! Well, finding positive correlation among replicates seems to help!

ADD REPLY
0
Entering edit mode

I was surprised too :)

ADD REPLY
0
Entering edit mode

The paper you refer has a LOT of work, apart from the RNA-seq runs.

ADD REPLY

Login before adding your answer.

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