SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
question about affy microarray annotation on UCSC tujchl Bioinformatics 1 07-28-2013 01:23 AM
tophat with/without annotation,and cufflink with annotation? louis7781x Bioinformatics 19 04-05-2013 07:09 AM
Question about exon limits and annotation.. shyam_la Bioinformatics 16 07-26-2012 03:09 PM
question related to Annovar annotation kjaja Bioinformatics 4 03-27-2012 02:32 AM
genome annotation question slny Bioinformatics 3 06-01-2011 07:26 PM

Reply
 
Thread Tools
Old 04-27-2014, 11:36 AM   #1
bongbimit
Member
 
Location: Ha Noi

Join Date: Dec 2013
Posts: 29
Default a question about annotation data in tophat

Let me start by saying I am very new to this so I may be asking naive questions

i am align my rna-seq reads (human) to genome by tophat but i am confusing about some parameter of annotation data in tophat

-G/--GTF : it mean i also have genes.gtf and i understand i know the exon, transcript and also the splicing junction

But
--no-novel-juncs: i dont understand what 's it function because i said above i have the genes.gtf file

thanks you all.
bongbimit is offline   Reply With Quote
Old 04-27-2014, 01:35 PM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

A "novel junction" would be a new splice junction (i.e., a new connection between exons that doesn't exist in your GTF file). Finding novel junctions makes tophat take longer (it's slow enough as is, try STAR instead if you have enough memory) and may or may not be needed, depending on exactly what you're studying.
dpryan is offline   Reply With Quote
Old 04-27-2014, 06:38 PM   #3
bongbimit
Member
 
Location: Ha Noi

Join Date: Dec 2013
Posts: 29
Default

Quote:
Originally Posted by dpryan View Post
A "novel junction" would be a new splice junction (i.e., a new connection between exons that doesn't exist in your GTF file). Finding novel junctions makes tophat take longer (it's slow enough as is, try STAR instead if you have enough memory) and may or may not be needed, depending on exactly what you're studying.
but i think when i set parameter -G genes.gtf, tophat will map reads only into exons, other unmapping-reads will be slided smaller and map back into exon

so what is the different between from 2 parameter [-G --no-novel-junc] and only [-G]?
bongbimit is offline   Reply With Quote
Old 04-27-2014, 07:02 PM   #4
blancha
Senior Member
 
Location: Montreal

Join Date: May 2013
Posts: 367
Default

With the option --GTF or -G, TopHat will not only map the reads to the exons annotated in the GTF file. TopHat will first map the reads to the exons annotated in the GTF file, then TopHat will map the unmapped reads to the rest of the genome. You need to add the parameter -T/--transcriptome-only to map the reads only to the transcriptome defined in the GTF file.

The difference when you add the parameter --no-novel-juncs is that TopHat will not try to splice the reads at junctions other than those defined in the GTF file.

I generally set the parameter --no-novel-juncs since the runtime of TopHat is much longer without this parameter, as mentioned by dpryan, and there are more false positive mappings. However, if you are interested in discovering novel splicing junctions, not already discovered and annotated in the GTF file, you should leave this setting on.

Last edited by blancha; 04-27-2014 at 07:20 PM.
blancha is offline   Reply With Quote
Old 04-27-2014, 09:48 PM   #5
bongbimit
Member
 
Location: Ha Noi

Join Date: Dec 2013
Posts: 29
Default

Quote:
Originally Posted by blancha View Post
With the option --GTF or -G, TopHat will not only map the reads to the exons annotated in the GTF file. TopHat will first map the reads to the exons annotated in the GTF file, then TopHat will map the unmapped reads to the rest of the genome. You need to add the parameter -T/--transcriptome-only to map the reads only to the transcriptome defined in the GTF file.

The difference when you add the parameter --no-novel-juncs is that TopHat will not try to splice the reads at junctions other than those defined in the GTF file.

I generally set the parameter --no-novel-juncs since the runtime of TopHat is much longer without this parameter, as mentioned by dpryan, and there are more false positive mappings. However, if you are interested in discovering novel splicing junctions, not already discovered and annotated in the GTF file, you should leave this setting on.
Thanks you all, i got it
bongbimit is offline   Reply With Quote
Reply

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 07:18 PM.


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