Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Cufflinks2 analysis without significant genes

    Dear all,

    I'm currently analysing RNAseq data of two distinct cancer conditions using the tuxedo pipeline and I highly doubt my results so far. According to cuffdiff, there are absolutely NO genes found to be differentially expressed with a significant p-value which I cannot truly believe. When evaluating my results with cummeRbund, more than 60% of the transcripts are flagged with test status other than "OK".

    Therefore, I would like to ask whether anyone has encountered similar problems and how to circumvent them. Since I already heard that cufflinks seems to have serious flaws, I run scripture in parallel although I still have to find tools to further analyse the assembled transcripts's abundance.

    So here some details on my analysis:
    • paired end RNAseq reads with 101 bases length
    • mapping done via bowtie2/tophat2 against Hg19
    • using latest releases of bowtie2, tophat2 and cufflinks2

    here's the commands I use for tuxedo:
    • tophat --b2-sensitive -p 8 -o mapping1 bowtie2-2.0.0-beta6/hg19/hg19 R1_001.fastq R2_001.fastq
    • cufflinks --GTF UCSC_annotation_hg19.gtf -p 8 -o annotation1 mapping1/accepted_hits.bam
    • cuffmerge -g UCSC_annotation_hg19.gtf -s bowtie2-2.0.0-beta6/hg19/hg19.fa -p 8 assemblies.txt
    • cuffdiff -m 199 -s 38 -o diff_out -b bowtie2-2.0.0-beta6/hg19/hg19.fa -p 9 -L G1,G2 -u merged_asm/merged.gtf mapping1/accepted_hits.bam,mapping2/accepted_hits.bam,mapping3/accepted_hits.bam mapping4/accepted_hits.bam,mapping5/accepted_hits.bam,mapping6/accepted_hits.bam


    note: I currently use the --GTF flag for cufflinks, though I will change that to --GTF-guide as soon as the results make sense.

    and for scripture:
    • for i in {1..22};
      do java -jar scripture-beta2.jar -alignment merged_align.bam -out chr$i.scriptureESTest.segments -sizeFile Hg19/chr$i\_size.txt -chr chr$i -chrSequence Hg19/Homo_sapiens.GRCh37.67.dna_rm.chromosome.$i.fa;
      done


    Any help and suggestions are greatly appreciated!
    Best regards

  • #2
    you aren't alone in thinking the cufflinks/cuffdiff output doesn't make sense. i've never been happy with it either with novel isoform discovery and quantification or -G quantifaction only mode. if you're not after isoform level DE then you can go with raw read counts at the gene level and then use an R package like DESeq or EBSeq to analyze for differential expression. you'll need something to count your reads (you can try HTSeq, authored by the author of DESeq). HTSeq is a python library that comes with a script called htseq-count.

    even if you ARE after isoform level DE it might be helpful to at least start at the gene/locus level using this alternative approach.

    another possible pipeline is RSEM / EBseq. RSEM works by quantifying isoform level expression against a known annotation (any GTF). EBSeq is a new R package that takes isoform level uncertainty into account when calculating differential expression. I have not yet run all the way through this pipeline yet so I can't say how the output looks. of course with this pipeline you lose novel isoform discovery. i'd think you could take the output of Scripture and build a GTF which you could then use with RSEM for quantification.
    /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
    Salk Institute for Biological Studies, La Jolla, CA, USA */

    Comment


    • #3
      Yes, we have been experiencing similar problems (along with other problems) with the latest version of Cufflinks and are considering other methods for our DE analysis. I have some variance in my biological replicates which may be an issue for cuffdiff, but I have data for example where the RPKMs for Disease v Control are 250 vs 10 and this is not significant.

      Here are some of the issues we have been trying to deal with. I have also emailed the cufflinks helpdesk but have not had a response.

      http://seqanswers.com/forums/showthread.php?t=23962
      http://seqanswers.com/forums/showthread.php?t=23910
      http://seqanswers.com/forums/showthread.php?t=21708

      Comment


      • #4
        Nice to see that I am not the only one experiencing this (though I did already find a lot of related threads).

        I was hoping to get some replies by the cufflinks team as well, but so far it seems that they are busy preparing the next release version.

        Despite that, my current plan is to switch to easyRNASeq for read counting and edgeR as well as DEXSeq for downstream analysis of my samples.
        Thus far, transcript assembly is only a future prospective, but it will definitely end up on my table and scripture might be an option once I get it working properly. Unfortunately, the last time I tried to run it, scripture did not proof to be very user friendly.

        Comment

        Latest Articles

        Collapse

        • 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
        • seqadmin
          The Impact of AI in Genomic Medicine
          by seqadmin



          Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
          02-26-2024, 02:07 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 03-14-2024, 06:13 AM
        0 responses
        34 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-08-2024, 08:03 AM
        0 responses
        72 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-07-2024, 08:13 AM
        0 responses
        82 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-06-2024, 09:51 AM
        0 responses
        68 views
        0 likes
        Last Post seqadmin  
        Working...
        X