Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Cuffdiff normalization issue

    I think that CuffDiff (from cufflinks ) might have a normalization issue!

    I'm using the Tophat/cufflinks/cuffdiff pipeline to analyze RNA-seq data from a WT and knockout (KO) mouse. Both done in dublicates.

    But when I make the scatterplot of WT vs KO using the FPKM values from Cuffdiff there is a general shift towards WT being more expressed than KO as seen here:

    Click image for larger version

Name:	CuffDiifScatter.jpeg
Views:	1
Size:	108.0 KB
ID:	308021

    I've plotted each replicates combination separately to be consistent with the plots below. The red line is the 1:1 ratio (where the bulk of the values are expected to be) and the blue line is the dy/dx mean created with geom_smooth() from ggplot2 (Generalized additive models with integrated smoothness estimation).

    This pattern probably originates from the raw read counts as the shift is also seen in the raw data, as showed here:

    Click image for larger version

Name:	RawScatter.jpeg
Views:	1
Size:	132.9 KB
ID:	308022

    Even though the shift is not as pronounced in the raw data.

    I think it is a normalization problem because when I take the raw read counts and normalize them using EdgeR ( calcNormFactors() ) the shift is much smaller, as seen here:

    Click image for larger version

Name:	EdgeRscatter.jpeg
Views:	1
Size:	132.5 KB
ID:	308023

    I'm using Cufflinks v 2.0.2 and the iGenome as refrences.
    I get these results regardless of how I run Cuffdiff. I've tried using the following parameter combination:

    only mandatory parameters
    -N
    -N -M myGTF
    -N -M myGTF -b
    -N -M myGTF -b -u

    And I always get the same results. Furthermore I have had a similar problem with another mouse data-set created by a completely different experimental protocol.

    Does anyone have good suggestions of how to make CuffDiff work?
    Last edited by Kristoffer Knudsen; 02-11-2013, 06:10 AM.

Latest Articles

Collapse

  • seqadmin
    Genetic Variation in Immunogenetics and Antibody Diversity
    by seqadmin



    The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
    11-06-2024, 07:24 PM
  • seqadmin
    Choosing Between NGS and qPCR
    by seqadmin



    Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
    10-18-2024, 07:11 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 11-08-2024, 11:09 AM
0 responses
225 views
0 likes
Last Post seqadmin  
Started by seqadmin, 11-08-2024, 06:13 AM
0 responses
166 views
0 likes
Last Post seqadmin  
Started by seqadmin, 11-01-2024, 06:09 AM
0 responses
80 views
0 likes
Last Post seqadmin  
Started by seqadmin, 10-30-2024, 05:31 AM
0 responses
27 views
0 likes
Last Post seqadmin  
Working...
X