Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Tophat2 giving unreasonably long junctions

    I'm using Tophat2 to align my RNA-seq reads. Many of the junctions produced by Tophat have extremely long introns and span multiple genes in my dataset. I've set the max intron length to be slightly larger than the biggest intron in my GFF file, but I'm still observing this weird behavior. I've attached a few IGV screenshots of what I'm talking about. These very long "intron" mapped reads occur consistently throughout my dataset. Often, however, the junctions perfectly match with my gene model (also in an IGV screenshot), suggesting that something here is working.

    Has anyone else observed this behavior?

    EDIT:
    Discussion with a lab mate has solved the problem. I'm adding the solution here in case anyone runs into the same issue and searches for it. The GFF file I was using was incorrectly formatted, and Tophat was interpreting each exon as a separate transcript. Lesson learned - make sure your GFF is correctly formatted in the first place and if you're concerned, double check the reconstructed transcriptome file.
    Attached Files
    Last edited by allind; 08-28-2013, 10:41 AM. Reason: Solved problem

  • #2
    Your second screen shot looks especially weird because it looks like Tophat is assigning splice junctions between divergently transcribed genes!

    Comment

    Latest Articles

    Collapse

    • 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
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 08:47 AM
    0 responses
    12 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    60 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    59 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    54 views
    0 likes
    Last Post seqadmin  
    Working...
    X