Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • d f
    Member
    • Feb 2010
    • 18

    #1

    how does cuffcompare choose which transcript to put in combined.gtf file?

    When cuffcompare is run with multiple input gtf files from different sequencing expts, how does it decide which transcript to be the "representative" transcript to list in the combined.gtf file? If multiple transcripts have the same intron structure but are not identical, which one does it choose to put in combined.gtf? It does not always choose the longest one as the "representative" transcript (which is what I thought it would do).

    For example, from these two input gtf files,

    Code:
    contig00177     Cufflinks       transcript      230     711     1000    .       .       gene_id "CUFF.68933"; transcript_id "CUFF.68933.1"; FPKM "107.2929288684"; frac "1.000000"; conf_lo "86.576469"; conf_hi "128.009389"; cov "2.269504";
    contig00177     Cufflinks       exon    230     711     1000    .       .       gene_id "CUFF.68933"; transcript_id "CUFF.68933.1"; exon_number "1"; FPKM "107.2929288684"; frac "1.000000"; conf_lo "86.576469"; conf_hi "128.009389"; cov "2.269504";
    and

    Code:
    contig00177     Cufflinks       transcript      230     1047    1000    .       .       gene_id "CUFF.71009"; transcript_id "CUFF.71009.1"; FPKM "86.1874620509"; frac "1.000000"; conf_lo "67.620022"; conf_hi "104.754903"; cov "2.055336";
    contig00177     Cufflinks       exon    230     1047    1000    .       .       gene_id "CUFF.71009"; transcript_id "CUFF.71009.1"; exon_number "1"; FPKM "86.1874620509"; frac "1.000000"; conf_lo "67.620022"; conf_hi "104.754903"; cov "2.055336";
    the transcript listed in the combined.gtf file is the first and shorter transcript:

    Code:
    contig00177     Cufflinks       exon    230     711     .       ^@      .       gene_id "XLOC_000001"; transcript_id "TCONS_00000001"; exon_number "1"; oId "CUFF.68933.1"; class_code ".";
    I am running version 0.9.2 of cufflinks/cuffcompare, and am running it without a "reference" annotation GTF.

    Incidentally, my goal here is to create a transcriptome for a mostly unannotated, novel genome. I have RNA-seq data from two different sequencing runs. I ran TopHat (with a couple of parameter sets) and Cufflinks on the reads to predict transcripts. I am now using cuffcompare to consolidate the results of the transcript predictions from the two sequencing/TopHat runs to create the most "comprehensive" transcriptome.

Latest Articles

Collapse

  • SEQadmin2
    Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
    by SEQadmin2



    CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

    Despite this, “CRISPR helped turn genome editing from a specialized technique into
    ...
    07-31-2026, 11:01 AM
  • SEQadmin2
    Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
    by SEQadmin2


    Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

    The systematic characterization of the human proteome has
    ...
    07-20-2026, 11:48 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 08-06-2026, 07:41 AM
0 responses
14 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-03-2026, 10:13 AM
0 responses
31 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 07-31-2026, 02:55 AM
0 responses
40 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 07-24-2026, 12:17 PM
0 responses
26 views
0 reactions
Last Post SEQadmin2  
Working...