SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat RNASeq mapping - right reads map and left reads map 50% less airad22 Bioinformatics 2 08-14-2013 07:26 AM
TopHat: SE reads maria_mari RNA Sequencing 4 12-10-2012 07:55 AM
unexpected BWA output when enabling multihit jorge Bioinformatics 1 08-31-2011 01:03 AM
retain multihit reads for CNV analysis jorge Bioinformatics 0 08-30-2011 12:41 AM
Genomic reads and Junction reads-Tophat repinementer Bioinformatics 0 08-10-2010 11:05 PM

Reply
 
Thread Tools
Old 02-03-2015, 07:21 AM   #1
allerzbul
Junior Member
 
Location: KCMO

Join Date: Feb 2015
Posts: 1
Default Tophat Multihit Reads

Hi, I'm running the top hat aligner for a ribo-depleted human RNAseq dataset. I'm interested in looking at expression in snRNAs. These loci in particular are very repetitive and when I run Tophat using the default settings (allowing up to 20 hits for each read and reporting these secondary alignments), none of these loci are assigned any single mapping reads. I originally wanted to assign the other multi mappers using the single read count ratios (there is an option for this in Cufflinks) but since there are no single mappers, this probably won't work. I tried an alternative approach, setting Tophat -g 1, which based on my understanding assigns each read to it's best mapping location (based on the mapping score) or randomly, if all of the map scores are the same. This seems to have worked for my snRNAs. For the loci i've looked at, it seems that their ratios when compared to each other are about the same as when I allowed multi-mapping.

I want to know if using the -g 1 has any downfalls? It seems to work really well (based on my initial look) and it removes a lot of the complications that come with dealing with multi mapping reads in differential expression analyses. But if it does work so well, why is default to allow 20 multi mappers? Also, has anyone else come up with different ways to deal with multi mapping reads in differential expression analyses?

Thanks so much for your help!
allerzbul is offline   Reply With Quote
Old 02-24-2015, 05:36 AM   #2
swaraj
Member
 
Location: Naples, Italy

Join Date: Feb 2012
Posts: 50
Default

This issue is discussed in a publication (http://www.biomedcentral.com/1471-2164/12/516), which reports that accepting multi mapping reads (up to 10) is beneficial (Fig 1). It accounts for maximum
number of gene detection especially genes with low expression levels. The authors argue
that multihits allows the inclusion of homologous sites from the RNAseq data. Using a
reference GTF file with predicted or known transcript models can alleviate the problem. But it is not an ideal case for studies where transcription of repeats is the focal point.
Secondly this combination increases the false negative transcript population
drastically.
swaraj is offline   Reply With Quote
Reply

Tags
multi-mapping, tophat

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:38 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO