Where does my RNASeq contamination fits on the Tree of Life?
I have RNASeq of algal cultures so my samples are not really axenic. Instead, the presence of some bacteria or fungus is to be expected. I assembled my reads with Trinity and now I would like to estimate the origin of each individual contig. Ideally, I would like to get visualization of where the contamination is coming from on the Tree of Life as a:
- quality metrics that the origin of my contamination makes sense (and I see what I expect to see for algal cultures)
- to remove contaminants and "clean" the assembly
Is there any tool that could do this for me?
I started by automatically outputting the "best" blast hit for each contig, but I am getting large variety of the hits and I am not sure how to summarize them or properly assign them phylogenetically.
Thanks for help.
• 2,848 views
•
link
written
by
Biomonika (Noolean)
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
De novo transcriptome assembly (Trinity) with heavy fungal contamination in fruit/flower samples
written by soundslikedaisy •
-
How to handle contamination in 16s amplicon data
written by 1769mkc
-
Remove microbes(cobionts) contamination: before or after insect de novo genome assembly?
written by hankool •
-
How does BWA handle duplications?
written by Nikolas •
-
Could adaptor contamination and primer dimers lead to low number of cells detected during a single …
written by charlesgwellem •
-
Unequally pooling libraries for RNAseq?
written by jingjin2203 •
-
Trinity transcriptome quality assessment
written by cwbenson1993 •
-
Bowtie2 orientation mapping question
written by Andrés Ribone •
-
Puzzle - Strange Contig Assembly In Newbler - Trimmed Ends
written by Biomonika (Noolean)
-
Retain Capitalization For Features
written by Lee Katz
Why not start with filtering out reads that can be mapped to known bacterial/fungal species?
Can you point me to such list/database?
I'm a big fan of Kraken for screening against contamination, the program assigns a taxid to each read, with a little leg work you could filter off of that. If you use the kraken-translate tool you should be able to get the whole taxonomy for each read and filter there using keywords. E.g. get a list of reads with the word "bacteria" in their kraken-translate entry, then toss all of those reads from your reads.
https://ccb.jhu.edu/software/kraken/MANUAL.html#output-format
I am honestly not sure which is better: clean before assembly or after assembly.
NCBI has a new
ref_prok_rep(representative prokaryotic genomes) pre-made blast database available. Since you have assembled sequences you could do a quick blast against that to see if you can find any low hanging fruits in terms of identification.