Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • cufflinks : analysis comparison with and without a gtf reference file

    Hi,
    I have many questions about cufflinks output. Here one of them :
    First I've used tophat to map my RNAseq (100pb) to obtain a accepted_hits.bam file.
    Then I've used cufflinks in two ways :
    - simply : cufflinks accepted_hits.bam
    - with a gtf file, that is the actually annotation of my genome (eucalyptus) :
    cufflinks -g annotGenome.gtf accepted_hits.bam
    Note that I've used the –g and not the –G option.

    One example result :
    - without reference gtf :
    one gene / one isoform : 12110-17714
    The first part of this isoform 12-16530 has the same structure intron/exon than isoforms formed with the reference.
    Then I have a last exon 16530-17714.

    -with reference gtf
    one gene/two isoforms
    + transcript 1: 12024-17350 = exact transcript from the reference
    full_read_support "no";
    The corresponding no reference last exon is now :
    16530 - 16561
    16597 - 17350
    That's my reference, but in my run this intron is mapped. There is no read that split in two parts. A few reads begin at position 16595. I've checked no read ending at 16561.
    I thing this RNA doesn't exist in my transcriptome.
    + transcript 2 : 12024-17714
    full_read_support "yes";
    last exon : 16530-17714, the same exon than the no reference version
    Why this transcript2 contains the 12024-12109 portion that is not mapped with RNAseq (instead the reference=transcript1 begin with this sequence) ?

    for the two isoforms, I have FPKM values (4 for transcript1 that doesn't seem to exist in my transcriptome and 13 for the transcript2). How cufflinks attributes those values ?

    With the version without gtf reference, I have a FPKM=36, that is the double comparing with the version with reference (13+4=17) while the mapping file is the same.

    At least, note that those transcripts are located on the forward strand of the genome and that there is nothing in gtf and cufflink results on the opposite strand at this location.

    Many thanks for your suggestion,

    S.

  • #2
    Hello sohnic,

    I have no answer to your question, but I would like to ask you one because you're posting about the very same issue I am wondering about.

    Can you just give me some hints about the difference of using (or not) a local reference to run cufflinks?
    I would like to understand what does cufflinks use to assembly the information contained in target.bam file when you don't provide it a local reference. Further, I'd also like to know whether using it will make cufflinks name the transcripts with the gene names instead of CUFF.65279 and the like.

    Thank you for your help.

    Comment


    • #3
      My experience with cufflinks

      Cufflinks will use the mapped reads and the reference genome to create a parsimonious assembly (a minimum spanning tree-esque structure) that explains the reads. The logic of the algorithm is based on the Burset and Guigo paper (Burset, Guigo, Evaluation of Gene Structure Prediction, 1996) and their observations regarding high false positive rates of gene finding programs that are guided by annotation. This is my understanding of the matter based on the past year I have spent working on similar issues.

      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
      8 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      8 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      49 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      66 views
      0 likes
      Last Post seqadmin  
      Working...
      X