Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Cuffdiff with more than two samples

    Recently,I'm using Cufflinks to process my RNA-seq data,and I have 4 tissues(more than two).In the mannual of cuffdiff,you said,if more than two are provided, Cuffdiff tests for differential expression and regulation between all pairs of samples. But when I used the cuffdiff,it only report one pairs of samples.There was only one gene_exp.diff in the output directory.Command was just as below,
    ======================================
    cuffdiff -o DEG_cuffdiff -p 20 cuffcmp.combined.gtf ./tophat_sd/accepted_hits.bam ./tophat_em/accepted_hits.bam tophat_en1/accepted_hits.bam ./tophat_en2/accepted_hits.bam
    ======================================
    What's my problem in using cuffdiff?
    Also I want to compare gene expression of duplicate gene pairs, and is it neccesery to remove the reads that mapped to multiple location in the BAM file generated by Tophat as duplicate gene pairs have high sequence similarity ?
    Thanks!

  • #2
    For Question 1:

    Yes, Cuffdiff would only generate one output even though you have more than one comparisons. However, all other comparisons are also included in that file. Check the output file more carefully.

    Comment


    • #3
      I have checked the outpu directory.

      In fact, I have checked the outpu directory, and there are only 11 files as follows:cds.diff , cds_exp.diff,genes.fpkm_tracking, isoforms.fpkm_tracking ,splicing.diff, tss_groups.fpkm_tracking , cds.fpkm_tracking , gene_exp.diff isoform_exp.diff , promoters.diff , tss_group_exp.diff. And I did not find any other comparison files between genes.
      I have posted this problem on the other thread ,and his met the similar issue.
      Application of sequencing to RNA analysis (RNA-Seq, whole transcriptome, SAGE, expression analysis, novel organism mining, splice variants)

      Comment


      • #4
        I am not quite sure right now for your case. But I did this before and I can get all the comparisons. For example, if I open the gene_exp.diff file, I will have the q1 vs q2 first. And when I continue check the result by clicking pagedown, I will find q1 vs q3 and etc.

        Comment


        • #5
          Thanks for your help

          Now I know all of results of the comparisons betweem pairs are saved in one file.
          Thanks!

          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
          18 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 10:19 PM
          0 responses
          22 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 09:21 AM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-04-2024, 09:00 AM
          0 responses
          47 views
          0 likes
          Last Post seqadmin  
          Working...
          X