Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Is cufflinks still useful if you have a reference genome?

    Sorry, the actual heading for my post should be:

    If you have a high quality genome annotation, should you always still use cufflinks.

    Disclaimer: I suspect I'm missing something so this may be a silly question but I've been thinking about this for years and am still puzzled so here goes.

    I don't think I should use cufflinks for my expt because:

    a) I am working on species with high quality annotation (human).

    b) Since I'm working with short read fragments (few hundred bp), it is likely that using cufflinks to redefine transcripts would introduce as many/more errors in transcript assembly than it would true novel transcripts, simply because transcripts are best described by long (sanger) sequences and only so much information can be inferred about the which transcript a short read originated from. Sometimes it is impossible to determine which differentially exons at the 5' end are connected to others at the 3' end. This is no fault of the algorithm.

    c) Because cuffdiff calculates gene-level expression by summing isoform-level expression, any false/artifactual isoforms may reduce the certainty with which a read can be associated with any individual isoform. This comes down to the principle that the more isoforms there are, the more difficult it is for cuffdiff to accurately deconvolute which isoform the read originated from (as per cuffdiff manual?).

    d) aside, I'm inclined to think that if I don't intend to validate novel genes/isoforms then I should stick with published ones - global-level RNA-seq of known genes should be just that.

    So, in a nutshell this leads me to thinking the following:

    If you're not specifically interested in novel gene/isoform discovery and suspect that you would do more harm than good to your genome annotation by including transcripts which were not supported by Sanger-length sequencing, then this should be your pipeline:

    tophat -> bam -> + reference.gtf -> cuffdiff

    i.e. no cufflinks step.

    Yet every protocol seems to have a cufflinks step, and the -G parameter seems to be the used a lot in this kind of situation (despite the potential issues I mentioned above).

    Can anyone comment whether I have missed something/made a fundamental mistake or does my reasoning seem sound? Should I be running cufflinks for general-purpose human genome RNA-seq?

    P.S. This is assuming you wan't differential expression, of course. If you only have 1 sample, of course you run cufflinks.
    Last edited by KAP; 04-17-2013, 04:03 AM. Reason: typo

  • #2
    If I only want a list of DEG I go straight from Tophat to Cuffdiff (skipping de novo steps using Cufflinks and cuffmerge). This is what is written in the cuffdiff manual, so I assume it must be valid!

    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
    49 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