Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • zeam
    Member
    • Oct 2010
    • 43

    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!
  • taozuo
    Junior Member
    • Jul 2011
    • 8

    #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

    • zeam
      Member
      • Oct 2010
      • 43

      #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

      • taozuo
        Junior Member
        • Jul 2011
        • 8

        #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

        • zeam
          Member
          • Oct 2010
          • 43

          #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

          • SEQadmin2
            From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
            by SEQadmin2


            Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


            The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
            ...
            Yesterday, 10:05 AM
          • SEQadmin2
            Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends
            by SEQadmin2


            With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.


            Introduction

            Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
            05-22-2026, 06:42 AM
          • SEQadmin2
            Environmental Genomics in the Age of NGS: From Microbes to Conservation Strategies
            by SEQadmin2

            Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.


            Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...
            05-06-2026, 09:04 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, Yesterday, 12:03 PM
          0 responses
          19 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, Yesterday, 11:40 AM
          0 responses
          14 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 05-28-2026, 11:40 AM
          0 responses
          29 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 05-26-2026, 10:12 AM
          0 responses
          31 views
          0 reactions
          Last Post SEQadmin2  
          Working...