Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Cuffdiff - results replicates

    Hi:

    I have 4 different cell types. A , B, C and D. Each class has 4 replicates.

    In Cuffdiff, I gave replicates with comma separated lines

    A1.bam, A2.bam, A3.bam,A4.bam B1.bam,b2.bam,b4.bam ... and so on.


    From gene_exp.diff, I wanted to get list of differentially expressed genes.
    I observe that cuffdiff also compared q1 q2 (which is a1.bam a2.bam).

    why cuffdiff comparing between two replicates. Could anyone provide their expert input.

    Also is this the only file to look for differentially expressed genes?

    Thanks
    Adrian

  • #2
    q1 and q2 will be the replicate groups, not the sample id's if you've specified your replicate groupings correctly.

    Use -L to specify labels instead

    genes_exp.diff is the gene level summaries of the differential expression testing between conditions.

    The other diff files are summaries of the differential expression at transcript leve, tss level, cds level etc.

    Comment


    • #3
      Hi :
      Thanks for your reply.

      A1 A2 A3 B1 B2 B3 C1 C2 C3 D1 D2 D3

      If q1, q2, Q3 and Q4 are replicate groups, I should not get anything above q4 since I have only 4 replicate groups. However, I am seeing in columns Sample_1 and Sample_2:

      q1 q4

      q1 q8

      q1 q9

      q1 q12

      etc.

      This lead me to think that cuffdiff is comparing every replicate sample to other samples in other groups and reporting the significant FDR (yes or no calls).

      Yes, I did not specify -L option.


      Is the following correct syntax for -L option?

      cuffdiff combined.gtf -L sample1, sample2, sample3 and sample4 a1.bam,a2.bam,a3.bam b1.bam,b2.bam,b3.bam c1.bam,c2.bam,c3.bam d1.bam,d2.bam,d3.bam


      Thanks a lot again.

      -Adrian

      Comment


      • #4
        cuffdiff combined.gtf -L sample1,sample2,sample3,sample4 a1.bam,a2.bam,a3.bam b1.bam,b2.bam,b3.bam c1.bam,c2.bam,c3.bam d1.bam,d2.bam,d3.bam

        This should be the right one.

        Comment


        • #5
          Thanks a lot for your reply.

          I have one more question:

          Is is okay if I do cuffdiff only on two replicate groups.

          Say cuffdiff between a1,a2,a3 and b1,b2,b3 only.

          since my cuffmerge gtf assembly file is made from all samples in the experiment, will I loose power by only comparing two replicate groups and not all.

          Thanks again.

          -Adrian

          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