Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • FPKM values differ for sample sample between comparisons?

    Hi I have three biological samples (lets call them Sample_A-C)that I want to find DEGs between. Firstly, I'm performing pairwise comparisons between samples using cuffdiff using alignments generated using tophat.

    The problem is when I compare the fpkm value of genes in Sample_A (in the gene_exp.diff file), the values differ depending on whether the comparison of Sample_A was against Sample_B or Sample_C.

    How is this difference generated? I thought when running cuffdiff the fpkm value would be static and would originate from the mapping data in the bam/sam file used as input for cuffdiff.

    Here is just the header and the first gene in the gene_exp.diff file.
    Sample_A vs Sample_B
    test_id gene_id gene locus sample_1 sample_2 status value_1 value_2 log2(fold_change) test_stat p_value q_value significant
    ACYPI000001 ACYPI000001 acyp2eg0000191 chr1:2365250-2374181 q1 q2 OK 21.6621 25.2228 0.219556 0.389799 0.71015 0.991694 no
    Sample_A vs Sample_C
    test_id gene_id gene locus sample_1 sample_2 status value_1 value_2 log2(fold_change) test_stat p_value q_value significant
    ACYPI000001 ACYPI000001 acyp2eg0000191 chr1:2365250-2374181 q1 q2 OK 26.0189 587.227 4.49629 4.16814 0.0011 0.422675 no

    So as you can see, when compared against Sample_B, the fpkm value of ACYPI000001 in Sample_A is 21.6621 and in the second comparison it's 26.0189

    There are even more pronounced differences in these datasets and I just don't get why these differences are arising.

    Can someone please shed some light on what may be underlying this?

    Thanks

    EDIT: I'm using cufflinks version 2.1.1
    cmd# cuffdiff -p 5 -o sample_A_vs_sample_B annotations.gff sample_A_accepted_hits.bam sample_B_accepted_hits.bam
    Last edited by beeman; 04-16-2014, 07:06 PM. Reason: adding more info

  • #2
    You might want to read about 'Library Normalization Methods' on the cufflinks manual:



    FPKMs and fragment counts are scaled via the median of the geometric means of fragment counts across all libraries, as described in Anders and Huber (Genome Biology, 2010). This policy identical to the one used by DESeq. (default for Cuffdiff)

    Comment


    • #3
      Great thanks for the info

      Comment


      • #4
        Hello,
        Did you try to compare the three conditions in the same process with cuffdiff ?
        As explained, Cuffdiff should support that with "cuffdiff sample1 sample2 sample3". The 3 conditions should be compared instead of doing several pairwise comparison. In my case, it doesn't work ! Do you know if some options should be used to do that.
        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
        31 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        33 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
        53 views
        0 likes
        Last Post seqadmin  
        Working...
        X