Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Run cufflinks with or without annotation?

    Hi all
    From what I see in this forum, the common pipeline for Cufflinks is to run it without a reference annotation and supply the annotation only to cuffcompare.

    If I run cufflinks with annotation, it will not build novel transcripts, but are there any additional considerations?
    i.e. If I don't care about novel transcripts, is it better/recommended for me to run cufflinks with annotation?

    thanks in advance,
    Reut

  • #2
    Dear Reut,

    The latest version of Cufflinks allow you to supply an annotation file and with two options: consider the reads that only compatible with the annotation, or building novel transcripts in addition to the annotation.

    Check the online manual for details.

    Douglas

    Comment


    • #3
      I would think it should be better to run it with an annotation if you are not interested in novel transcripts.

      Comment


      • #4
        I run it with a reference. In fact, unless you have very good coverage, I think running it with a reference is a better way to go. Unless you have split or paired end reads linking every part of a gene together, running Cufflinks without a reference annotation could cause it to split a single gene into multiple parts. Avoiding an artifact like that is more important to me than using Cufflinks to identify novel transcripts. There are other tools to use for that.

        Comment


        • #5
          Thanks!

          Thank you everyone and thank you pbluescript, what you are describing is exactly what I see - when running Cufflinks without reference I see a single gene split into multiple parts.
          So I will run Cufflinks only with reference.

          Comment


          • #6
            Originally posted by reut View Post
            Thank you everyone and thank you pbluescript, what you are describing is exactly what I see - when running Cufflinks without reference I see a single gene split into multiple parts.
            So I will run Cufflinks only with reference.
            I have a question here.

            As mentioned in the mannul of cufflinks, "Output will include all reference transcripts as well as any novel genes and isoforms that are assembled." is for the "-g/--GTF-guide" option.

            My question is if the predicted transcripts don't include the
            reference transcripts
            , how could I differentiate the predicted transcripts from the output of cufflinks using the -g option? It seems I could only see the novel ones, compared to the reference.

            And how could you prove the "multiple parts" are artifacts and remove them?

            Cheers,

            Comment


            • #7
              My question is if the predicted transcripts don't include the reference transcripts, how could I differentiate the predicted transcripts from the output of cufflinks using the -g option?
              From the manual, it sounds like it will output all provided reference transcripts in addition to any novel isoforms. If I'm understanding your question correctly, it sounds like you'd be able to identify novel isoforms by subtracting the reference GTF file you provide from the Cufflinks output GTF.

              And how could you prove the "multiple parts" are artifacts and remove them?
              If you load the Cufflinks output GTF file into say, the UCSC genome browser, and look at your favorite gene, you'll see that a multi-exonic transcript appears as several, smaller "fragment" transcripts--some single exon, some just a few exons--but all incomplete as compared to the reference. In theory, you could strip all of these by comparison to the reference set of transcripts, but then you'd eliminate the novel isoform discovery function of Cufflinks, and might as well run with that function disabled. The conclusion reut and pbluescript reached is that, when run with the --GTF-guide option, most of these spurious transcripts will be eliminated while still allowing novel isoform discovery.
              Last edited by polyatail; 06-01-2011, 12:57 PM. Reason: with -> when

              Comment


              • #8
                Has anyone found any differences in Cufflinks output when providing a reference via -G at the TopHat step? With or without also providing a reference to Cufflinks?

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Current Approaches to Protein Sequencing
                  by seqadmin


                  Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                  04-04-2024, 04:25 PM
                • 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

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 04-11-2024, 12:08 PM
                0 responses
                30 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 04-10-2024, 10:19 PM
                0 responses
                32 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 04-10-2024, 09:21 AM
                0 responses
                28 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 04-04-2024, 09:00 AM
                0 responses
                52 views
                0 likes
                Last Post seqadmin  
                Working...
                X