Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • What's the best RNA Seq Tuxedo suite tutorial/streamlined intro protocol?

    Hi all, new to this community. I've been attempting to learn the Tuxedo suite using their protocol "Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks." (http://www.ncbi.nlm.nih.gov/pubmed/22383036) but keep dead ending with errors. Is there a better tutorial, more up to date and streamlined tutorial out there? Thanks

  • #2
    What errors are you getting? I had followed the protocol with newer versions of tuxedo suite a few months ago (at least to cuffdiff part) and recall that it had worked.

    Comment


    • #3
      One of the GFF ID is duplicated, but even our resident lab coder couldn't find the problem. A lot of things aren't explained all that well in the protocol, hence my troubles. I'd really like to find an additional resource to aid my understanding.

      Code:
      cuffmerge -g genes.gtf -s genome.fa -p 8 assemblies.txt
      
      [Thu Jan 22 11:26:25 2015] Beginning transcriptome assembly merge
      -------------------------------------------
      
      [Thu Jan 22 11:26:25 2015] Preparing output location ./merged_asm/
      [Thu Jan 22 11:26:29 2015] Converting GTF files to SAM
      [11:26:30] Loading reference annotation.
      [11:26:31] Loading reference annotation.
      [11:26:33] Loading reference annotation.
      [11:26:34] Loading reference annotation.
      [11:26:36] Loading reference annotation.
      [11:26:38] Loading reference annotation.
      [Thu Jan 22 11:26:40 2015] Quantitating transcripts
      Warning: Your version of Cufflinks is not up-to-date. It is recommended that you upgrade to Cufflinks v2.2.1 to benefit from the most recent features and bug fixes (http://cufflinks.cbcb.umd.edu).
      Command line:
      cufflinks -o ./merged_asm/ -F 0.05 -g genes.gtf -q --overhang-tolerance 200 --library-type=transfrags -A 0.0 --min-frags-per-transfrag 0 --no-5-extend -p 8 ./merged_asm/tmp/mergeSam_tmp.6.EZGXyZ 
      [bam_header_read] EOF marker is absent. The input is probably truncated.
      [bam_header_read] invalid BAM binary header (this is not a BAM file).
      File ./merged_asm/tmp/mergeSam_tmp.6.EZGXyZ doesn't appear to be a valid BAM file, trying SAM...
      [11:26:40] Loading reference annotation.
      [11:26:52] Inspecting reads and determining fragment length distribution.
      Processed 11167 loci.                       
      > Map Properties:
      >	Normalized Map Mass: 68011.00
      >	Raw Map Mass: 68011.00
      >	Fragment Length Distribution: Truncated Gaussian (default)
      >	              Default Mean: 200
      >	           Default Std Dev: 80
      [11:26:54] Assembling transcripts and estimating abundances.
      Processed 11167 loci.                       
      [Thu Jan 22 11:27:27 2015] Comparing against reference file genes.gtf
      Warning: Your version of Cufflinks is not up-to-date. It is recommended that you upgrade to Cufflinks v2.2.1 to benefit from the most recent features and bug fixes (http://cufflinks.cbcb.umd.edu).
      Error: duplicate GFF ID 'FBtr0079998' encountered!
      	[FAILED]
      Error: could not execute cuffcompare

      Comment


      • #4
        Just want to confirm that you are using the test data (and not your own). I will take another look later this evening at this.

        Comment


        • #5
          Originally posted by GenoMax View Post
          Just want to confirm that you are using the test data (and not your own). I will take another look later this evening at this.
          Yes I am.I downloaded GSE32038_simulated_fastq_files.tar.gz from http://www.ncbi.nlm.nih.gov/geo/quer...i?acc=GSE32038 and alligned it against the igenome http://support.illumina.com/sequenci...e/igenome.html BDGP5

          Comment


          • #6
            I tried cuffmerge (using cufflinks v.2.2.1). No problems in getting the command to complete.

            Comment


            • #7
              Originally posted by GenoMax View Post
              I tried cuffmerge (using cufflinks v.2.2.1). No problems in getting the command to complete.
              Thanks, I'll update it. Is there an easy command for updating?

              Comment


              • #8
                Are you administering your own machine or do you need to get someone else to do it on a central server/cluster?

                Comment


                • #9
                  Originally posted by GenoMax View Post
                  Are you administering your own machine or do you need to get someone else to do it on a central server/cluster?
                  It's our lab's own mini supercomputer, so I do have root access to it

                  Comment


                  • #10
                    You will have to download and compile it: https://github.com/cole-trapnell-lab/cufflinks Make necessary changes to your path so you start using the new version (if you are keeping the old version around).

                    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, 03-27-2024, 06:37 PM
                    0 responses
                    13 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-27-2024, 06:07 PM
                    0 responses
                    11 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-22-2024, 10:03 AM
                    0 responses
                    53 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-21-2024, 07:32 AM
                    0 responses
                    69 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X