Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • carmeyeii
    replied
    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?

    Leave a comment:


  • BAJ
    replied
    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:


  • CompBio
    replied
    Originally posted by Lee Sam View Post
    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.
    Actually, that's one of the approaches I'm comparing it to. Supersplat can use a lot of memory that appears to increase with the size of the query file as well as with parameter changes. If you're not careful to split up your files, your machine can spend a lot of time swapping.

    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:


  • Lee Sam
    replied
    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:


  • john_mu
    replied
    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:


  • CompBio
    replied
    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:


  • anecsulea
    replied
    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 Post
    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.

    Leave a comment:


  • IrisZhu
    replied
    Originally posted by CompBio View Post
    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?
    Could the reason be your "-g 1"? Is "1" too small.
    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:


  • CompBio
    started a topic TopHat Misses Splice Junctions

    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?

Latest Articles

Collapse

  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin


    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM
  • seqadmin
    Essential Discoveries and Tools in Epitranscriptomics
    by seqadmin




    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
    04-22-2024, 07:01 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
23 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
44 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-09-2024, 02:46 PM
0 responses
58 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-07-2024, 06:57 AM
0 responses
44 views
0 likes
Last Post seqadmin  
Working...
X