Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Tuxedo pipeline: cufflinks with -G option vs cuffdiff with reference

    Dear all,

    I have been analysing RNAseq experiments using the default tuxedo pipeline for about a year, i.e., TopHat2 + cufflinks + cuffmerge + cuffdiff2 and one thing I have never understood is the value that running cufflinks with the -G/--GTF option adds to the analysis.

    Considering that cufflinks acts as a transcript assembler, with the -G option it does not search for novel transcripts, which means that the output GTF contains the same transcripts as the original reference. Therefore, cuffmerge will merge as many reference annotation files as available samples which will produce, again, the reference annotation. The only difference is that if cuffmerge is run including the fasta reference, it will add some attributes that allow cuffdiff to test for differential expression at the TSS and CDS level.

    Considering all this, my question is whether there is any big difference between running the whole pipeline without searching for novel isoforms (cufflinks -G) and running directly cuffdiff with the reference annotation directly. I must say that I have tried this in the past and that the gene level results where slightly different, but most importantly, the second option is much faster than the whole pipeline.

    Any input will be highly appreciated.

    Thanks,

    Asier

  • #2
    Just like you discovered some folks are not interested in finding new transcripts/splice isoforms and are happy to use the "known transcriptome" for the expression analysis. So doing a special one-time TopHat run with the -G option to create a transcriptome specific index subset of genome speeds the analysis process up.

    Comment


    • #3
      Thank you GenoMax,

      So would you say there are any differences between running only TopHat + CuffDiff and TopHat + Cufflinks (-G) + Cuffmerge + CuffDiff if you only want to stick to the annotated transcriptome?

      I am asking because I work with wheat, which has a 16 Gb genome (12 Gb latest reference) and skipping cufflinks and cuffmerge means that you can save up to one week depending on the experiment.

      Asier

      Comment


      • #4
        If you are not interested in novel isoforms, but you are interested in saving time, it's much faster to skip TopHat and map to the transcriptome only, using a different aligner. Or even better, an alignment-free tool such as Seal or Sailfish.

        Comment


        • #5
          Thanks for your input Brian,

          I accept that using the transcriptome as the reference is much quicker and the tools like Sailfish and Kallisto are amazing in that sense. However, mapping to the transcritome means that if a gene is not annotated properly, say it is a partial gene, a read could map beyond its boundaries in the genome but would not map in the transcriptome.

          Having saying this, I have used the whole Tuxedo pipeline searching for new isoforms, but cufflinks is not flawless. One of the problems it has is that it may merge genes thar are close.

          Considering all these, I still try to understand if running cufflinks with the -G option has any advantages over skipping it.

          Thanks again,

          Asier

          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