Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • aguffanti
    replied
    Cuffllinks variance, gene expression level and ifferential expression

    Hello ! which kind of sequencing was that - fragment or paired ends ? it changes a lot

    You have to look at the *count.tracking files and examine the count variance values

    There are instances (especially with fragment sequencing) for which the count variance is far greater than the counts - actually I would like to know from Cole wether this is somewhat expected

    Here is an example (see the values for q2), with the associated gene expression reported below - looks like for q2 a FPKM of 0 is calculated due to the enormous variance, am I right ?

    tracking_id q1_count q1_count_variance q1_count_uncertainty_var q1_count_dispersion_var q1_status q2_count q2_count_variance q2_count_uncertainty_var q2_count_dispersion_var q2_status

    GAPDH 308519 40753000 0 308519 OK 207559 10265500000 0 10265500000 OK


    test_id gene locus status MCF7 MCFS MCF7+MCFS log2(fold_change) test_stat p_value q_value significant

    NM_001256799 GAPDH chr12:6643584-6647537 OK 11.25 0.00 11.25 =-inf nan 0.00005 0.0142079 yes

    I am running software versions as follows on Color Space reads (so I have to stick to the old bowtie..)

    bowtie version 1.0.0
    TopHat v2.0.9
    cufflinks v2.1.1

    Thanks in advance for any feedback,

    Alessandro

    Leave a comment:


  • jp.
    replied
    Hi
    May please someone check my attached cuff_diff results. I have problem with high error bar and low FPKM value. Is there something wrong with sequencing or its a cufflinks problem?
    Second, there are only few genes in either samples (not in all) with 0 FPKM value; for example, control 3840.28 FPKM and treated 0 FPKM with Log2_fold_change=Inf. Attached File 0_FPKM.png Expression bar plot indicates high fold expression ? Can I ignore this one which seems significant ?

    Thank you very much

    Click image for larger version

Name:	1 csDensity.png
Views:	1
Size:	15.5 KB
ID:	304263

    Click image for larger version

Name:	csBoxplot_genes.png
Views:	1
Size:	11.6 KB
ID:	304264

    Click image for larger version

Name:	dispersionPlot_genes.png
Views:	1
Size:	22.5 KB
ID:	304265

    Click image for larger version

Name:	seq.png
Views:	3
Size:	13.3 KB
ID:	304266

    Click image for larger version

Name:	0_FPKM.png
Views:	1
Size:	6.2 KB
ID:	304267
    Last edited by jp.; 07-24-2013, 10:49 PM. Reason: adding info

    Leave a comment:


  • jwaage
    replied
    Hi all; I too, have some similar problems - I'm running the newest cuffdiff (2.0.2) on a two-sample, two-replicates illumina run using

    cuffdiff --mask-file xxx.gtf -upper-quartile-norm

    genes.diff and isoforms.diff are fine (although not many are called DE), but splicing.diff and promoters.diff only have NOTEST, LOWDATA or FAIL.

    Any suggestions? I've tried running with -c 1, but to no avail.

    All the best,
    Johannes Waage
    Uni of. Copenhagen

    Leave a comment:


  • billstevens
    replied
    Originally posted by MeixiaZhao View Post
    Hello,
    I used the latest version cufflinks 2.0.2, comparing with my previous results from 1.3.0, the results are unbelievable:
    Gene level: I got 340 significant genes with cuffdiff2 and 2692 with cuffdiff1.3;
    Isoform level: I got 272 significant isoforms with cuffdiff2 and 11380 with cuffdiff1.3;
    My samples are single replicate.
    Any solutions to judge the results? Or which version should I use for later analysis?
    Thanks a lot!
    I've spoken to Cole about this extensively, and he says the new version is just much more accurate. It is an absolute shock going from one to the other.

    Leave a comment:


  • MeixiaZhao
    replied
    Hello,
    I used the latest version cufflinks 2.0.2, comparing with my previous results from 1.3.0, the results are unbelievable:
    Gene level: I got 340 significant genes with cuffdiff2 and 2692 with cuffdiff1.3;
    Isoform level: I got 272 significant isoforms with cuffdiff2 and 11380 with cuffdiff1.3;
    My samples are single replicate.
    Any solutions to judge the results? Or which version should I use for later analysis?
    Thanks a lot!

    Leave a comment:


  • billstevens
    replied
    Originally posted by Cole Trapnell View Post
    That's excellent - thanks to you and others in this thread for the helpful feedback. The official release of 2.0.1 is imminent.
    Hi Cole,

    Thanks for all the notes. However, I'm getting the same issue here. I have 46,000 reads that test, and 40,000 that have NOTEST. I am using the -b option, but not -c 0, and not the multi-read correct or frag bias. I'm using Cufflinks 2.0.1 and Tophat 1.4. Was the shortcoming in regards to -b addressed in the newest Cufflinks?

    Leave a comment:


  • sdodson
    replied
    Hi,

    I'm a novice in RNA-seq analysis, and I have a question regarding Cuffdiff...

    We generated our RNA-seq data using the Illumina HiSeq2000, and have been following the data analysis protocol for Tophat, Cufflinks, ect. from Nature Protocols. We are currently using Tophat v1.4.1 and Cufflinks v1.3.0. Our Cuffdiff output shows 10371 of 27080 genes are NOTEST. In an effect to decrease this value, can we re-run our data through the newest version of Cuffdiff without first re-running it through the newest version of Cufflinks? Or is something else we can change to minimize the number of NOTEST genes? (such as using -c 0 ?)

    Thanks!

    Leave a comment:


  • Cole Trapnell
    replied
    Originally posted by pinin4fjords View Post
    I also confirm that the numbers produced using the new version at #32 look more sensible (using -b and min-outlier-p).
    That's excellent - thanks to you and others in this thread for the helpful feedback. The official release of 2.0.1 is imminent.

    Leave a comment:


  • Cole Trapnell
    replied
    Originally posted by gesdy View Post
    Hi Cole,
    I used the cuffdiff version you posted, and now it works.
    I repeated the same analysis with the new version and the 1.3 (cuffdiff)
    I got 43 genes called significant with cuffdiff2 and 940 with cuffdiff1.3
    I have just one replicate for each conditions.
    Do you think it's normal?
    here my command line:

    cuffdiff -b genome.fa -p 10 -u genes.gtf

    any suggestions?
    thank you very much!
    Mat
    Hi Mat,

    That sounds fair to me - Cuffdiff 2 is far more conservative at the gene level than was 1.3 (as detailed on the site), and with a single replicate per sample, even more so. This isn't too surprising to me.

    Leave a comment:


  • pinin4fjords
    replied
    I also confirm that the numbers produced using the new version at #32 look more sensible (using -b and min-outlier-p).

    Leave a comment:


  • gesdy
    replied
    Hi Cole,
    I used the cuffdiff version you posted, and now it works.
    I repeated the same analysis with the new version and the 1.3 (cuffdiff)
    I got 43 genes called significant with cuffdiff2 and 940 with cuffdiff1.3
    I have just one replicate for each conditions.
    Do you think it's normal?
    here my command line:

    cuffdiff -b genome.fa -p 10 -u genes.gtf

    any suggestions?
    thank you very much!
    Mat

    Leave a comment:


  • Cole Trapnell
    replied
    Hmm, we see that sometimes when one or more of the libraries have extremely low sequencing yield, for example. Can you show me a csDendro plot with replicates=T? It's important to verify that the replicates for each sample really segregate together. It might also be helpful to see a dispersion plot for each condition. If you're worried about revealing sample names you can email it to [email protected].

    Leave a comment:


  • glados
    replied
    Cole, can you please take a look at this boxplot of one of the genes. I don't think the error bars are due to variances in my samples.
    Attached Files

    Leave a comment:


  • NicoBxl
    replied
    I've exactly the same problem with the -b param.

    Leave a comment:


  • pinin4fjords
    replied
    Originally posted by Cole Trapnell View Post
    We are testing a pre-release that fixes the issues reported so far, at least in our hands. Can you try:



    (No mac build or source yet - we were only able to reproduce this with our pre-compiled linux binary anyways)
    Okay, will do. Do you still recommend the min-outlier-p option? Have corrections been made to Cuffdiff only, or is it necessary to re-run the whole Cufflinks sequence?

    Leave a comment:

Latest Articles

Collapse

  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 04:11 PM
1 response
20 views
0 likes
Last Post kim897
by kim897
 
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
35 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
25 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-14-2024, 09:19 AM
0 responses
57 views
0 likes
Last Post seqadmin  
Working...
X