Hi everyone!
I have a lot of (short) RNAseq reads from a metagenome sample. I am only interested in the reads originating from one species in the microbial community, so i would like to filter my fastq files accordingly.
I have a reference genome of the species, so my naive approach would be to map all reads to this reference (let's say with bowtie2) and keep all reads that have mapped for further analysis.
However, this approach seems to introduce a bias that is much larger than i expected. I mapped my data twice: First to a bunch of reference genomes and then only to the single reference that i am interested in. In the second case, the amount of reads mapping to my genome of interest is much higher. Can you explain to me why this happens? Am i now also collecting all the reads of all related species?
Can i overcome the problem by adjusting parameters and making the mapping more strict? Or is the whole approach of mapping to only a single reference just not reasonable?
1 answer
I mapped my data twice: First to a bunch of reference genomes and then only to the single reference that i am interested in. In the second case, the amount of reads mapping to my genome of interest is much higher. Can you explain to me why this happens?
Aligners are generally designed to be greedy i.e. they will do their best to try and align a read. When you use a reduced representation of the reference (e.g. your single reference) it is possible that the aligner may shoe-horn a read in even though it may not have originated from that species. And then there is the genuine issue of multi-mapping where reads from similar genes from multiple genomes will map to that one reference. So yes you could be indeed be collecting reads from other species when you use a single reference.
Or is the whole approach of mapping to only a single reference just not reasonable?
Depends. Ultimately what you want to get out of this analysis may decide the answer.
The samples are from microbial communities, so there is no opportunity to put different tags on the individual species' DNA. Is the right course of action assembly everything in the sample, then mapping the reads, then extracting those which mapped contigs resembling my reference?
Log in to answer this question.
-
Tool to align metagenomic data to reference genome
written by Madde •
-
Subsampling reads
written by amy •
-
bbsplit not using all reads in library
written by ghd1998 •
-
What is an acceptable mapping rate on de novo transcriptomes for subsequent differential expression…
written by lebouquinaugustin •
-
Extracted mapped shotgun metagenomic reads to reference genome. SPAdes or metaSPAdes for de-novo as…
written by O.rka
-
Quick question: samtool idxstats and multiple mapping
written by Aurelie MLB
-
subsetting sets of scaffolds based on sequence similarity
written by second_exon
-
Better statergy for Gap Closing
written by second_exon
-
Should Rna-Seq Reads Map To One Strand Of Cdna Reference?
written by chjiao3456 •
-
Illumina Pe Reads Mapping To Reference Genome
written by rohita.sinha •
