Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • StringTie - isoforms not overlapping

    Hi,
    I am using StringTie for transcriptome reconstruction and identification of new isoforms.
    While I was exploring the "[file]-transcripts.gtf" output file from stringtie I found something that intrigued me... in the example below I show three isoforms resulting from the same gene ("STRG.14686") and the last one was present in the reference annotation. However the start and end coordinates do not match. The first two isoforms end at 394180 bp and 389422 bp, respectively, while the third starts at 396257 bp...


    scaffold_96 StringTie transcript 383404 394180 1000 + .gene_id "STRG.14686"; transcript_id "STRG.14686.1"; cov "72.829010"; FPKM "6.506798"; TPM "8.058224";
    scaffold_96 StringTie transcript 383404 389422 1000 + .gene_id "STRG.14686"; transcript_id "STRG.14686.2"; cov "61.675678"; FPKM "5.510321"; TPM "6.824155";
    scaffold_96 StringTie transcript 396257 398001 1000 + .gene_id "STRG.14686"; transcript_id "STRG.14686.3"; reference_id "scaffold_96.g39603.t1"; ref_gene_id "scaffold_96.g39603"; cov "2963.938721"; FPKM "264.808624"; TPM "327.947357";
    Why is StringTie "clustering" these isoforms in the same gene?
    Last edited by pbarros; 04-26-2017, 03:22 AM.

  • #2
    first of all let me say that I agree with your interpretation. the third transcript does not overlap the first two and it has given all three the same 'gene_id' value.

    the only thing that comes to mind is that the assembly from stringtie, or cufflinks for that matter, is an attempted explanation, and often a simplification, of the alignment data. you may learn more about this by looking at the alignments in an alignment browser such as UCSC or IGV.

    while i doubt this, stringtie could be very, very smart and have assembled the third isoform from reads that multimapped between it and the first two transcripts which would imply they are all the same gene but one that is repeated in more than one position.
    /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
    Salk Institute for Biological Studies, La Jolla, CA, USA */

    Comment


    • #3
      thank you for the input sdriscoll ... maybe I was overthinking this

      cheers,
      pedro

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Strategies for Sequencing Challenging Samples
        by seqadmin


        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
        03-22-2024, 06:39 AM
      • seqadmin
        Techniques and Challenges in Conservation Genomics
        by seqadmin



        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

        Avian Conservation
        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
        03-08-2024, 10:41 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 06:37 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      9 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      50 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Working...
      X