SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
one gene one RPKM cufflinks or DEseq frankyue50 Bioinformatics 5 10-30-2014 12:37 AM
some gene no RPKM in every library by Cufflinks stunlay RNA Sequencing 0 08-19-2011 09:06 AM
cufflinks problem: one gene have multiple rpkm values calliopsis Bioinformatics 1 05-25-2010 08:18 AM
tophat/cufflinks bam vs. RPKM mgogol Bioinformatics 5 04-26-2010 09:58 AM
Tophat : percentage of reads contributing to RPKM values sdriscoll Bioinformatics 19 10-16-2009 04:01 PM

Reply
 
Thread Tools
Old 10-23-2009, 08:43 AM   #1
mgogol
Senior Member
 
Location: Kansas City

Join Date: Mar 2008
Posts: 197
Question RPKM and multiple reads, tophat and cufflinks

My understanding of the RPKM calculation in Tophat is that it includes multi-reads that match < 40 times in the genome (by default).

It seems like cufflinks does something more complex involving some kind of allocation of multiple reads.

http://cufflinks.cbcb.umd.edu/manual.html

Cufflinks models the sequencing process by asking what the probability is of observing each read, given a set of transcripts and a set of abundances. The program then multiplies these probabilities to compute the overall likelihood that one would observe the reads in the experiment, given the proposed abundances on the transcripts. Because Cufflinks' statistical model is linear, the likelihood function has a unique maximum value, and Cufflinks finds it with a numerical optimization algorithm.


I know we're supposed to start using cufflinks' RPKM now, but I'd like to understand tophat's as well. Does anyone know if my description is correct?
mgogol is offline   Reply With Quote
Old 10-23-2009, 12:31 PM   #2
Cole Trapnell
Senior Member
 
Location: Boston, MA

Join Date: Nov 2008
Posts: 212
Default

Correct - both TopHat and Cufflinks include multireads in their RPKM calculations
Cole Trapnell is offline   Reply With Quote
Old 11-16-2009, 06:49 AM   #3
sandmann
Junior Member
 
Location: Germany

Join Date: May 2009
Posts: 3
Default

I would like to use Cufflinks to assemble SOLiD reads, so unfortunately I can't use bowtie. Are there any recommendations about the maximum number of mappings any read should not exceed to be included in Cufflink's analysis ? In other words, should I discard reads that map e.g. to >40 locations in the genome before running Cufflinks ?

Any advice is greatly appreciated.
sandmann is offline   Reply With Quote
Old 12-03-2009, 07:46 AM   #4
Cole Trapnell
Senior Member
 
Location: Boston, MA

Join Date: Nov 2008
Posts: 212
Default

You should try leaving them in, but make sure their mapping quality is set to 0. Cufflinks will use them to assemble transcripts, but they will contribute little to nothing to RPKM's of transcripts that contain them. Leaving them out could create voids in read coverage where there are transcribed repeats (as happens not uncommonly in UTRs, etc), which could break up the gene models.
Cole Trapnell is offline   Reply With Quote
Old 12-07-2009, 08:05 AM   #5
HTS
Member
 
Location: Toronto

Join Date: Nov 2009
Posts: 24
Default

Please allow me to ask two related questions:

1. If I do use TopHat and Bowtie to do the alignment, will the mapping qualities of multireads be automatically adjusted? From the Bowtie manual it seems that mapping quality is always set to 255 at the moment.

2. Is the default value of 40 a good one for mammalian genomes? I have been using 10 in my own calculation so far and I am a bit curious about this default value.

Thanks a lot!

-- Leo
HTS is offline   Reply With Quote
Old 12-07-2009, 09:16 AM   #6
Cole Trapnell
Senior Member
 
Location: Boston, MA

Join Date: Nov 2008
Posts: 212
Default

Quote:
Originally Posted by HTS View Post
Please allow me to ask two related questions:

1. If I do use TopHat and Bowtie to do the alignment, will the mapping qualities of multireads be automatically adjusted? From the Bowtie manual it seems that mapping quality is always set to 255 at the moment.
Yes, TopHat calculates its own mapping quality. After ranking the alignments for each read, TopHat outputs all the alignments in the top rank. If there are more than one, then the probability that the read is incorrect (i.e. the mapping quality) is (1 - 1/p), assuming one of the alignments in the top rank *is* the correct one.

Quote:
Originally Posted by HTS View Post
2. Is the default value of 40 a good one for mammalian genomes? I have been using 10 in my own calculation so far and I am a bit curious about this default value.
It should be fine - Reducing it might cut down on your running time a bit, but not drastically in my opinion.
Cole Trapnell is offline   Reply With Quote
Old 09-03-2012, 02:46 PM   #7
carmeyeii
Senior Member
 
Location: Mexico

Join Date: Mar 2011
Posts: 137
Default

Hi,

VEEERY basic question: When calculating FPKM, what is the "correct" M value to use: Total unique mapped reads (i.e., the awk-ed file of accepted hits), or ALL the hits, i.e, all the alignments that TopHat output? What does cufflinks take as the M in FPKM?

Thanks for your help.

Carmen
carmeyeii is offline   Reply With Quote
Reply

Tags
cufflinks, multireads, rpkm, solexa, 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 04:52 AM.


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