Hi!
This is a somewhat old thread, but I would like to know more about the biological constraints TopHat uses to call a splice junction... is there anyway to override this?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Hi,
I would be interested in how you actually compared the different algorithms?
You say you used an existing data set... Why don't you use a constructed one? and how do you know what is true or not?
Could you also please summarize your findings?
As I am planing to do similar things within the next few weeks I would be of course very interested in any advice you could have.
Best,
Bernd
Leave a comment:
-
Originally posted by Lee Sam View PostThere is also the consideration that (to my admittedly shallow understanding) TopHat only uses the set of canonical splice signal sequences. Unfortunately, I can't recommend anything better - I was playing with the supersplat aligner which does spliced alignment but the program was very very slow and never ran to completion.
As for TopHat, I'll have to dig into their documentation again to see if I can figure out how it uses the canonical sites. Another good idea -- thanks!
Leave a comment:
-
There is also the consideration that (to my admittedly shallow understanding) TopHat only uses the set of canonical splice signal sequences. Unfortunately, I can't recommend anything better - I was playing with the supersplat aligner which does spliced alignment but the program was very very slow and never ran to completion.
Leave a comment:
-
There is a new program called HMMsplicer that was designed for the Arabidopsis genome. Maybe you can give that a try.
I haven't personally used it, but they were benchmarking it against SpilceMap and TopHat before.
Leave a comment:
-
Good suggestion -- relaxing the mismatches requirement may not help much if I don't allow a few multi-reads. I'll apply them both and see what happens.
Leave a comment:
-
Hi,
I don't know very much about Arabidopsis biology, but in mammals there is one class of genes for which TopHat is unable to find the junctions: the genes that generate duplicates through retrotransposition. This is because these duplicate copies are intronless, and thus the reads that
map on the splice junctions of the parental genes also map perfectly on the genome in the duplicated region. Since TopHat only uses reads that were not mapped on the genome in order to detect junctions, for this class of genes the junctions will not be detected.
So it might be worth to check if these genes have (recent) retrocopies in Arabidposis.
Hope this helps !
Best,
Anamaria
Originally posted by CompBio View PostI'm currently running a comparison between TopHat and other methods for finding splicing patterns in RNA Seq data. I've got about 250 million 32-nt reads from Arabidopsis. For the genes I am interested in, using other methods I have been able to align reads across many known and novel splice junctions at 100% ID that are unique to these locations (i.e., no multi-hits), with at least 8nt on either side of each junction (anchor regions). However, I am unable to get TopHat to find these junctions.
Leave a comment:
-
Originally posted by CompBio View PostI'm currently running a comparison between TopHat and other methods for finding splicing patterns in RNA Seq data. I've got about 250 million 32-nt reads from Arabidopsis. For the genes I am interested in, using other methods I have been able to align reads across many known and novel splice junctions at 100% ID that are unique to these locations (i.e., no multi-hits), with at least 8nt on either side of each junction (anchor regions). However, I am unable to get TopHat to find these junctions.
These are the TopHat parameters I'm using:
-g 1 # report only unique hits
-F 0.01 # report even poorly-represented junctions
--segment-mismatches=0 # enforce 100% ID for all reads
--splice-mismatches=0 # enforce 100% ID across junctions
--min-coverage-intron=10 # minimum allowed intron size for Arabidopsis
--max-coverage-intron=11000 # maximum intron size found in Arabidopsis
-i 10 # (as above)
-I 11000 # (as above)
--min-segment-intron=10 # (as above)
--max-segment-intron=11000 # (as above)
-j TAIR9_GFF3_genes.juncs # pre-processed splice junctions from gene model
-a 8 # minimum overlap/anchor of 8nt
-p 4 # allowing up to 4 threads on 8-processor machines
I can relax the constraint on splice mismatches to see if it helps, but ultimately I would like TopHat to find the junctions at 100% ID that I've already seen using other methods. I would like to make my comparison as fair as possible.
Any ideas? Am I misinterpreting/misusing any of these parameters?
I am always a bit confused about this "-g", coz the default is 40 --- I can't imagine it will keep all the reads that mapped to less than 40 places in the genome.
Leave a comment:
-
TopHat Misses Splice Junctions
I'm currently running a comparison between TopHat and other methods for finding splicing patterns in RNA Seq data. I've got about 250 million 32-nt reads from Arabidopsis. For the genes I am interested in, using other methods I have been able to align reads across many known and novel splice junctions at 100% ID that are unique to these locations (i.e., no multi-hits), with at least 8nt on either side of each junction (anchor regions). However, I am unable to get TopHat to find these junctions.
These are the TopHat parameters I'm using:
-g 1 # report only unique hits
-F 0.01 # report even poorly-represented junctions
--segment-mismatches=0 # enforce 100% ID for all reads
--splice-mismatches=0 # enforce 100% ID across junctions
--min-coverage-intron=10 # minimum allowed intron size for Arabidopsis
--max-coverage-intron=11000 # maximum intron size found in Arabidopsis
-i 10 # (as above)
-I 11000 # (as above)
--min-segment-intron=10 # (as above)
--max-segment-intron=11000 # (as above)
-j TAIR9_GFF3_genes.juncs # pre-processed splice junctions from gene model
-a 8 # minimum overlap/anchor of 8nt
-p 4 # allowing up to 4 threads on 8-processor machines
I can relax the constraint on splice mismatches to see if it helps, but ultimately I would like TopHat to find the junctions at 100% ID that I've already seen using other methods. I would like to make my comparison as fair as possible.
Any ideas? Am I misinterpreting/misusing any of these parameters?Tags: None
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
159 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: