Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Cuffcompare not showing all query transcripts

    Hi,

    I have 3 reference annotation GTF files and 3 query GTF files which were generated by Tophat-Cufflinks combination.

    I am using cuffcompare to compare the 3 query GTF files against each reference annotation GTF file.

    While doing this I noticed that the resulting combined.gtf and .tracking file are not reporting all the transcripts from the query GTF files.

    When I checked the .stats file I noticed in 'Summary for dataset' section that only a subset of the transcripts are reported.

    Example:
    GTF_A - 265244 transcripts
    GTF_B - 267648 transcripts

    .stats file for GTF_A reports always the same number against all references
    #= Summary for dataset: GTF_A.gtf :
    # Query mRNAs : 216023 in 57092 loci (188694 multi-exon transcripts)

    .stats file for GTF_B also reports always the same number against all references
    #= Summary for dataset: GTF_B.gtf :
    # Query mRNAs : 218034 in 58267 loci (189378 multi-exon transcripts)

    And even though the reference annotation is different the list of missing transcripts from query GTF is exactly the same.
    So they must have been skipped when importing the transcripts from query GTF and it is independent of the reference GTF.

    When comparing GTF_B (Query) vs GTF_A (Reference).
    .stats file reports

    #= Summary for dataset: GTF_B.gtf :
    # Query mRNAs : 218034 in 58267 loci (189378 multi-exon transcripts)
    # (21371 multi-transcript loci, ~3.7 transcripts per locus)
    # Reference mRNAs : 238139 in 45259 loci (223605 multi-exon)

    One thing to notice above is that the statistics for GTF_A (Reference) are different when it is used a reference from earlier when it was used in query.

    I tried to find any trend for the missing transcripts but cannot identify anything.

    Can anyone help me understand the issue at hand.

    Regards,

    Veer

  • #2
    I suspect this may be related to cuffcompare excluding "contained" transcripts. Maybe using the -C flag will give less mysterious results.

    From the manual:
    "-C Enables the "contained" transcripts to be also written in the <outprefix>.combined.gtffile, with the attribute "contained_in" showing the first container transfrag found. By default, without this option, cuffcompare does not write in that file isoforms that were found to be fully contained/covered (with the same compatible intron structure) by other transfrags in the same locus."

    Comment


    • #3
      I am sorry I should have posted this earlier. I later found out the reason for this behavior, it was because there were predicted transcripts in the GTF which did not have any strand information and these transcripts without strand information were skipped.

      So what I did is... I assigned them as + strand :-p. I decided to check for the correct strand info later.

      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, 03-27-2024, 06:37 PM
      0 responses
      12 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-27-2024, 06:07 PM
      0 responses
      11 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      53 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      69 views
      0 likes
      Last Post seqadmin  
      Working...
      X